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A  theoretical  development  of  small  amplitude  wave  disturbances 
on  oceanic  fronts  is  presented.  The  prototype  front  is  the  inshore 
boundary  of  the  Gulf  Stream,  although  the  model  encompasses  a  wider 
range  of  applicability.  This  work  is  an  extension  of  recent  research 
by  Garvine  which  showed  the  importance  of  dissipation  near  the  surface 
front,  and  earlier  work  by  Duxbury  on  large  scale  time  varying  inviscid 
flow  regimes.  A  two- layer  model  is  considered  in  which  the  lower  layer 
is  much  deeper  than  the  upper,  lighter  layer.  The  upper  layer,  includ¬ 
ing  the  frontal  zone,  is  then  dynamically  uncoupled  from  the  lower 
layer.  The  frontal  zone  that  forms  the  horizontal  boundary  of  the  upper 
layer  is  divided  into  two  regions,  a  smaller,  inner  region  in  which  the 
flow  is  dissipative  and  depends  upon  Interfacial  mass  entrainment  and 
turbulent  friction,  and  a  larger,  outer  region  in  which  the  motion  is 
Inviscid.  The  boundary  between  these  two  regions  is  placed  where  the 
internal  Froude  number  is  of  order  one. 

% 

Both  geostrophlc  and  ageostrophlc  basic  state  (non  time  varying) 
flows  are  considered  under  the  f-plane  approximation.  The  geostrophlc 
case  is  a  modification  of  Duxbury 's  two-layer  model  to  include  the 
inner  dissipative  zone  described  by  Garvine.  The  ageostrophlc  basic 
state  is  similar  to  one  originally  suggested  by  Stommel  as  a  more 
realistic  model  of  the  geometry  of  the  Gulf  Stream  front.  Nonzero 
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cross-stream  flow  is  allowed  in  each  basic  state  which  is  con¬ 
sistent  with  mass  entrainment  within  the  inner  region. 

Small  perturbations  to  the  basic  state  are  then  assumed  with 
time-harmonic  and  space-harmonic  dependence  in  the  direction  parallel 
to  the  front.  A  system  of  coupled,  linear  ordinary  differential 
equations  is  derived  in  the  perturbation  variables.  Their  solution 
is  developed  using  a  small  parameter  expansion,  the  small  parameter 
being  the  ratio  of  the  cross-stream  to  the  downstream  basic  state 
velocities.  Application  of  appropriate  physical  boundary  conditions 
leads  to  zeroth  and  first  order  dispersion  relations.  To  lowest  order, 
for  both  the  geostrophic  and  ageostrophic  basic  state  flows,  the  dis¬ 
persion  relation  is  linear.  All  waves  are  therefore  stable  to  lowest 
order,  and  a  continuum  of  wave  modes  is  allowed.  To  first  order,  that 
is,  when  finite  cross-stream  flow  occurs,  the  dispersion  characteris¬ 
tics  are  nonlinear  and  admit  complex  wave  frequencies  that  lead  to 
instabilities.  The  stability  characteristics  are  found  to  be  inde¬ 
pendent  of  the  shear  velocity  between  the  frontal  zone  and  ambient 
fluid  and  depend  only  upon  the  frontal  zone  interfacial  depth  at 
the  origin.  Comparison  is  made  between  Duxbury’s  analytical  model 
and  observational  data  for  Gulf  Stream  meanders.  Disturbance 
periods  generally  in  the  range  3-20  days  are  predicted  which  de¬ 
crease  with  increasing  latitude.  The  results  compare  well  with 
recent  satellite  data  for  the  Gulf  Stream. 
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CHAPTER  1 


INTRODUCTION 

1.1  OBSERVATIONS  OF  OCEANIC  FRONTS 

Oceanic  fronts  have  undoubtedly  been  observed  for  centuries  by 
those  who  sailed  the  sea.  Perhaps  they  have  even  been  understood  on 
an  Intuitive  level,  as,  for  example,  by  the  Japanese  fishermen  who 
regularly  exploited  their  biologically  rich  waters.  However,  their 
scientific  study  as  prominent  and  important  geophysical  features  is 
of  a  far  more  recent  origin.  One  of  the  earliest  published  accounts 
describing  what  we  recognize  today  as  an  oceanic  front  appeared  in 
William  Beebe's  book  The  Arcturus  Adventure  (1926),  in  which  he  devel¬ 
oped  the  metaphor  of  an  invisible  wall  separating  two  very  different 
water  masses.  Steaming  in  the  mid-Pacific  just  north  of  the  equator, 
Arcturus  encountered  the  filamentary  surface  manifestation  of  Beebe's 
invisible  wall,  a  region  he  described  as  "...a  world,  not  of  two,  but 
to  all  intents  and  purposes,  of  a  single  plane  -  length."  He  continued: 
"From  first  to  last  we  followed  its  course  along  a  hundred  miles ,  and 
yet  ten  yards  on  either  side  of  the  central  line  of  foam,  the  water 
was  almost  barren  of  life.  The  thread-like  artery  of  the  currents' 
juncture  seethed  with  organisms  -  literally  billions  of  living  creatures, 
clinging  to  its  erratic  angles  as  though  magnetized."  Beebe's  quali¬ 
tative  description  included  many  of  the  important  features  character¬ 
istic  of  oceanic  fronts,  features  which  have  been  observed  often  in 
subsequent  investigations. 
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Fronts  occur  at  the  boundary  between  two  water  masses  with 
different  physical  properties.  Thus,  they  generally  exhibit  a  long, 
thin  geometry.  Their  presence  is  frequently  revealed  at  the  sea- 
surface  by  a  foam  line  (as  observed  by  Beebe) ,  debris  accumulation, 
or  an  abrupt  change  in  some  visible  surface  property,  such  as  sea 
state,  water  color  or  transparency.  Although  none  of  these  surface 
manifestations  is  necessary  for  the  existence  of  a  front,  one  or  more 
of  them  often  accompanies  an  oceanic  front,  as  observed  by  Beebe  and 
others  (Amos,  et  al.,  1972;  Cromwell  and  Reid,  1956;  Knauss,  1957; 
Garvlne  and  Monk,  1974;  Stommel,  1976;  Zaneveld,  et  al. ,  1969). 

The  near-surface  zone  associated  with  an  oceanic  front  is  also 
frequently  characterized  by  considerably  enhanced  biological  activity 
compared  to  the  surrounding  waters.  Beebe  (1926)  considered  this 
discovery  one  of  the  most  remarkable  features  of  the  "current  rip" 
described  in  his  book,  and  attributed  the  large  number  and  variety 
of  living  organisms  to  the  shelter  provided  by  the  debris  accumulated 
along  the  foam  line.  There  is  not  necessarily  any  correlation  between 
biological  activity  and  debris  accumulation,  however,  since  fronts  have 
been  observed  in  which  the  former  was  present  but  not  the  latter 
(Knauss,  1957).  The  front  studied  by  Knauss  in  the  eastern  equatorial 
Pacific  was,  nevertheless,  easily  observed  by  eye.  The  colder  water  on 
one  side  was  darker  in  color  with  shorter-crested  waves  than  the  water 
on  the  warm  side,  where  most  of  the  biological  activity  was  concentrat¬ 
ed.  In  this  case,  the  biological  activity  was  apparently  promoted  by 


3 


higher  levels  of  suspended  and  dissolved  nutrients  near  the  surface 
front  than  in  the  surrounding  waters. 

Observational  data  also  show  that  fronts  tend  to  be  regions 
of  surface  convergence,  that  is,  surface  water  moves  toward  the  front 
from  both  sides,  subsequently  sinking  (Garvine  and  Monk,  1974;  Ingram, 
1976;  Voorhis,  1969;  Wooster,  1969).  No  observation  of  surface  diver¬ 
gence  has  apparently  been  made.  In  addition,  fronts  are  not  fixed  in 
space,  but  frequently  wander  about  the  ocean  (Garvine  and  Monk,  1974; 
Katz,  1969;  Wooster,  1969).  They  occur  on  a  variety  of  length  and 
time  scales,  from  tens  of  meters  to  thousands  of  kilometers  and  from 
hours  to  years  (Beebe,  1926;  Garvine  and  Monk,  1974;  Ingram,  1976; 
Knauss,  1957;  Pak  and  Zaneveld,  1974). 

In  view  of  the  number  of  characteristics  that  may  be  associated 
with  an  oceanic  front,  Cromwell  and  Reid  (1956)  provided  the  first 
precise  definition  of  the  term  in  an  attempt  to  introduce  a  standard¬ 
ized  nomenclature  and  to  isolate  the  essential  physical  processes. 

They  defined  an  oceanic  front  as  a  band  along  the  sea  surface  across 
which  abrupt  changes  in  the  water  density  occurred.  In  their  study, 
they  correlated  abrupt  temperature  changes  with  abrupt  density  changes, 
since  their  measurements  were  made  where  appreciable  salinity  contrast 
across  the  front  was  absent.  More  recent  studies,  however,  have  been 
made  of  fronts  with  nearly  uniform  density  as  a  result  of  offsetting 
temperature  and  salinity  effects  (Maul  and  Hansen,  1972),  and  of  fronts 
whose  density  contrast  was  a  consequence,  instead,  of  salinity  differ¬ 
ences  alone  (Amos,  ejt  al . ,  1972).  It  is  thus  apparent  that  Cromwell 
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and  Reid's  definition  should  now  be  broadened  to  include  these  other 
possibilities.  The  essential  notion  conveyed  by  the  term  "oceanic 
front"  is  that  the  front  constitutes  a  well-defined  boundary  between 
two  water  masses  with  different  physical  properties  (density,  temper¬ 
ature  or  salinity),  across  which  sharp  changes  occur  with  little  mixing 
between  the  two  regions  (Katz,  1969).  Changes  in  the  physical  proper¬ 
ties  of  the  water  masses  are,  moreover,  often  accompanied  by  sharp 
gradients  in  other  properties,  for  example,  current  velocity,  concen¬ 
tration  of  suspended  and  dissolved  materials,  water  color,  transparency, 
and  surface  wave  action. 

In  describing  the  Arc turns  observations,  Beebe  (1926)  referred 
to  the  surface  expression  of  the  water  mass  boundary  as  a  "current  rip", 
since  he  observed  a  strong  velocity  shear  parallel  to  the  foam  line, 
as  have  other  investigators  (Cromwell  and  Reid,  1956;  Zaneveld,  £t  al. , 
1969).  But  Cromwell  and  Reid  recognized  the  similarity  between  oceanic 
and  atmospheric  fronts  as  boundaries  between  fluid  masses  of  different 
density,  and  therefore  adopted  the  meteorologist’s  terminology.  However, 
their  definition  was  too  limited,  since,  in  addition  to  allowing  only 
density  differences,  it  considered  only  the  surface  expression  of  the 
boundary.  Instea',  the  surface  separating  the  two  water  masses  gener¬ 
ally  slopes  down  and  away  from  the  surface  front  into  the  ambient  water. 
In  this  paper,  therefore,  the  term  "front"  refers  to  the  entire  discon¬ 
tinuity  surface  that  separates  two  water  masses  with  different  proper¬ 


ties. 
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Three  observational  papers  that  provide  some  quantitative  insight 
into  typical  frontal  length  and  time  scales,  density,  temperature  and 
salinity  contrasts,  and  circulation  patterns  are  Katz  (1969),  Voorhis 
(1969)  and  Knauss  (1957).  These  papers  are  particularly  relevant,  since 
they  examine  typical  large  scale  upper  ocean  fronts,  and  the  theoreti¬ 
cal  work  presented  here  addresses  the  dynamics  of  such  large  scale 

O  O 

structures.  Katz'  investigation  involved  a  front  near  30  N,  70  W 

whose  linear  extent  was  at  least  200  km  oriented  approximately  north- 

west-southeast.  The  entire  frontal  zone  was  observed  to  translate  in 

a  west-southwest  direction  with  a  maximum  average  velocity  of  15  cm/s. 

In  addition,  a  definite  meander  (or  "wiggling")  of  the  frontal  boundary 

was  observed,  giving  rise  to  a  roughly  sinusoidal  surface  shape  of 

small  amplitude  and  about  180  km  wavelength.  Water  to  the  north  of 

the  frontal  boundary  was  found  to  be  characteristically  1°C  colder  than 

water  on  the  southern  side,  with  a  corresponding  salinity  contrast  of 

approximately  0.1°/oo.  The  variation  of  a  across  the  front  was  on  the 

order  of  0.25.  Over  a  surface  range  of  20  km,  the  depth  of  the  boundary 

between  the  warm  southern  water  and  the  colder  northern  water  increased 

from  zero  to  about  150  m,  corresponding  to  an  Interfacial  slope  on  the 
-2 

order  of  10  ,  and  thereafter  more  gradually  increased  to  about  175  m 

depth.  Measurements  were  made  with  a  depth  resolution  of  approximately 
2  m,  and  on  that  scale  no  mixing  was  observed  between  the  two  water 
masses.  The  frontal  boundary  was  therefore  very  sharp  and  well-defined. 
Evidence  was  found  for  a  circulation  pattern  that  included  convergence 
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of  che  surface  waters  followed  by  sinking  along  the  Interface  between 
the  two  water  masses.  The  frontal  zone  Rossby  number  was  estimated  to 
be  0.1,  indicating  nearly  geostrophic  flow.  Estimates  of  surface  veloc¬ 
ity  were  then  based  upon  geostrophy.  A  surface  velocity  shear  of 
approximately  80  cm/s  was  computed  and  found  to  agree  quite  favorably 
with  limited  direct  velocity  measurements. 

The  Sargasso  Sea  front  studied  by  Voorhis  in  1969  exhibited 
features  similar  to  those  described  above.  The  surface  expression  of 
the  front  was  tracked  over  1000  km  extent  and  was  marked  by  the  21.5*C 
and  23.5°C  isotherms,  with  the  2aC  temperature  contrast  being  confined 
to  a  region  sometimes  less  than  1  km  in  lateral  extent.  The  correspond¬ 
ing  salinity  and  ot  contrasts  were  found  to  be  approximately  0.1°/oo 
and  0.4,  respectively.  A  high  velocity  jet  was  observed  on  the  warm 
water  side  with  central  currents  as  high  as  80  cm/s.  The  oppositely 
directed  cold  water  current,  which  provided  the  characteristic  frontal 
zone  shear,  was  much  lower  at  10  to  30  cm/s.  The  warm-cold  water 
boundary  gradually  sloped  down  from  the  surface  to  a  depth  of  50  m  over 

_3 

a  horizontal  range  of  10  km  (slope  of  order  10  ) ,  with  the  less  dense 

warm  water  wedge  floating  on  the  more  dense  colder  ambient  water.  A 
drogue  <*tudy  provided  evidence  for  surface  convergence  along  the  frontal 
interface.  In  addition,  the  surface  front  was  marked  by  the  accumulation 
of  Sargassum  weed,  a  further  indication  of  surface  convergence.  Air¬ 
craft  infrared  studies  were  performed  during  a  five  month  period  in 
order  to  determine  frontal  zone  persistence  and  motion.  A  well-defined 
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front  was  observed  to  preserve  its  identity  throughout  this  period, 
but  with  marked  changes  in  its  configuration. 

The  front  studied  by  Knauss  in  the  Pacific  was  close  to  the 
equator  (2°46'N,  120®39'W)  corresponding  to  small  values  of  the 
Coriolis  parameter,  unlike  the  Sargasso  Sea  fronts.  Bathythermograph 
measurements  in  the  frontal  zone,  whose  surface  expression  was  confined 
to  a  band  some  100  yards  vide,  revealed  temperature,  salinity  and  o 
contrasts  of  about  3*C,  0.08°/oo  and  1.0,  respectively.  The  front  was 
observed  to  translate  normal  to  itself  at  a  speed  of  about  50  cm/s.  In 
addition,  surface  convergence  and  a  strong  velocity  shear  of  nearly 
2  m/s  were  observed. 

A  very  important  feature  of  oceanic  fronts,  and  the  one  addressed 
theoretically  in  this  paper,  is  their  tendency  to  develop  wavelike  dis¬ 
turbances  which  propagate  along  the  frontal  interface.  On  fronts 
marked  by  some  visible  surface  expression,  such  as  a  foam  or  debris 
line,  these  wavelike  distortions  of  the  boundary  between  two  water 
masses  are  frequently  manifested  as  roughly  sinusoidal  lateral  displace¬ 
ments  of  the  surface  foam  or  debris  line  (Zaneveld,  et  al.,  1969). 
Frontal  waves  have  been  observed  on  several  large  scale  upper  ocean 
frontal  systems,  for  example,  the  Kuroshio  Current  (Uda,  1964),  the 
California  Current  (Bernstein,  et  al. ,  1977),  Sargasso  Sea  fronts 
(Katz,  1969),  the  Norwegian  Current  (Mysak  and  Schott,  1977),  and  the 
Gulf  Stream  (Hansen,  1970;  Maul,  et  al.,  1978;  Stommel,  1976;  Webster, 
1961).  The  Gulf  Stream  Front  separates  cold  Atlantic  slope  water  from 


Che  warmer  Sargasso  Sea  water.  The  waves,  or  meanders  as  they  are 
frequently  called,  sometimes  grow  to  such  an  amplitude  that  vortex 
rings,  or  eddies,  actually  detach  from  the  front. 

In  his  study  of  Gulf  Stream  meanders  off  Onslow  Bay,  North  Carolina, 
Webster  (1961)  performed  a  spectral  analysis  that  revealed  dominant 
periods  of  6.9  and  3.9  days,  with  amplitudes  of  about  10  km.  Periods 
on  the  order  of  1-2  weeks  seem  to  be  characteristic  of  wave  distur¬ 
bances  on  the  Gulf  Stream  Front.  Hansen  (1970)  examined  wave  distur¬ 
bances  on  the  Gulf  Stream  east  of  Cape  Hatteras.  He  found  an  eastward 
propagating  wave  pattern  out  to  about  60°W  with  typical  wavelengths  in 
the  range  200-400  km  and  phase  speeds  on  the  order  of  5  cm/s.  He 
also  observed  a  generally  increasing  amplitude  with  the  wave's  eastward 
progression,  with  initially  small  disturbances  growing  to  some  200  km 
amplitude,  indicating  a  possible  instability  in  the  system.  Recent 
satellite  infrared  studies  of  the  Gulf  Stream  (Maul,  et  al. ,  1978) 
revealed  dominant  wave  periods  of  30  and  6  days  off  Onslow  Bay  and  45 
and  5  days  off  New  England.  Wavelengths  of  500-900  km  have  been 
observed  for  meanders  of  the  Kuroshio  system  (Uda,  1964);  300-500  km 
in  the  California  Current  (Bernstein,  et  al. ,  1977);  50-125  km  in 
the  Norwegian  Current  (Mysak  and  Schott,  1977);  and  approximately  100  km 
in  a  typical  Sargasso  Sea  front  (Katz,  1969). 

Atmospheric  frontal  systems  also  support  wavelike  disturbances 
which  are,  in  many  respects,  quite  similar  to  those  observed  on  oceanic 
fronts.  Oceanic  and  atmospheric  fronts  are  not  isomorphic,  however. 
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The  latter  may  be  conveniently  divided  into  two  categories,  lower  and 
upper  tropospheric  fronts. 

Lower  tropospheric  fronts,  the  first  category,  are  those  where 
a  thin  wedge  of  cold  air  at  the  earth's  surface  is  overlain  by  a  warm 
air  mass.  These  are  least  similar  to  upper  ocean  fronts.  For  the 
lower  tropospheric  front,  the  important  external  boundary  is  the 
earth's  surface.  There  the  no-slip  condition  must  be  satisfied  with 
the  consequent  development  of  a  frictional  boundary  layer.  The  hori¬ 
zontal  convergence  of  this  layer  at  the  surface  front,  in  turn,  pro¬ 
duces  a  vertical  component  of  flow  there.  In  the  oceanic  case,  on  the 
other  hand,  the  relevant  external  boundary  is  the  ocean's  free  surface, 
along  which  a  specified  pressure  distribution  must  be  maintained.  In 
the  absence  of  wind  shear,  frictional  surface  boundary  layers  will  not 
be  present,  so  that  some  other  physical  mechanism  is  required  to 
explain  the  observed  sinking  along  oceanic  frontal  interfaces. 

Upper  tropospheric  fronts  in  the  atmosphere,  the  jet  stream  or 
polar  front,  for  example,  are  far  more  similar  to  upper  ocean  fronts. 
Newton  (1978)  has  made  a  detailed  comparison  of  the  Gulf  Stream  front 
and  the  atmospheric  jet  stream  front  under  the  assumption  that  similar¬ 
ity  of  structure  implies  similarity  of  physical  processes.  The  Gulf 
Stream  and  jet  stream  fronts  are  geometrically  very  similar,  each 
having  a  high  velocity  filament  of  fluid  ("jet")  that  floats  buoyantly 
on  a  relatively  deep  layer  of  denser,  ambient  fluid  (the  "troposphere"). 
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The  fully  developed  frontal  zone  in  each  case  has  a  depth  which  is 

typically  10-20X  of  the  total  ambient  depth,  that  is,  on  the  order  of 

1.2  km  and  12  km  for  the  Gulf  Stream  and  jet  stream,  respectively.  The 

frontal  zone  width  is  likewise  in  the  ratio  1/10  for  the  oceanic  vs. 

the  atmospheric  front,  while  the  interfacial  slopes  are  of  the  same 
_2 

order  (10  )  for  each.  The  velocity  fields,  however,  are  considerably 

different,  with  the  maximum  velocity  in  the  jet  stream  being  some  30 
times  larger  than  in  the  Gulf  Stream.  Typical  wavelengths  for  dis¬ 
turbances  on  the  Gulf  Stream  are  on  the  order  of  200  km,  while  the 
corresponding  value  for  the  jet  stream  is  on  the  order  of  4000  km,  a 
ratio  of  1/20.  The  same  ratio  is  preserved  for  the  small  amplitude 
phase  velocities,  which,  according  to  Newton  (1978),  are  typically 
0.4  m/s  and  8  m/s.  The  amplitude-to-wavelength  ratio  in  both  cases  is 
approximately  1/10.  Small  amplitude  disturbances,  which  are  the  subject 
of  this  investigation,  thus  apparently  exhibit  considerably  larger  phase 
velocities  than  the  large  amplitude  waves  into  which  some  of  them 
eventually  evolve  (only  about  5  cm/s  for  the  Gulf  Stream  from  Hansen, 
1970). 

The  limited  observations  that  have  been  made  of  oceanic  fronts 
emphasize  the  variety  of  phenomena  associated  with  them  and  the  wide 
range  of  spatial  and  temporal  scales  on  which  they  exist.  Their  complex 
structure  and  dynamical  behavior,  however,  remain  largely  unexplained 
by  available  theory.  The  work  presented  in  this  paper  is  a  theoretical 


model  of  baroclinic  instabilities  on  large-scale  upper  ocean  density 
fronts.  Such  instabilities  produce  the  meanders  and  eddy  detachment 
observed  on  many  oceanic  frontal  systems.  In  view  of  the  nature  of 
this  work,  a  brief  review  of  earlier  theoretical  models  for  atmospheric 
and  oceanic  fronts  is  presented  in  the  next  two  sections. 

1.2  THEORETICAL  STUDIES  OF  ATMOSPHERIC  FRONTS 

In  a  theoretical  consideration  of  oceanic  fronts,  the  similarity 
to  their  atmospheric  counterparts  leads  naturally  to  a  review  of 
attempts  at  understanding  the  dynamical  processes  in  atmospheric  fronts. 
One  of  the  most  complete  discussions  of  the  theory  of  atmospheric  fronts 
was  provided  by  Stoker  (1953).  He  made  no  attempt  to  solve  the  frontal 
problem,  but  pursued  the  more  limited  objective  of  rigorously  formu¬ 
lating  the  full  nonlinear  problem.  Starting  with  the  fluid  dynamic 
equations  for  momentum  and  continuity,  as  well  as  certain  simplifying 
assumptions.  Stoker  developed  a  sequence  of  four  hydrodynamic  models 
of  decreasing  complexity,  each,  however,  being  more  restrictive  than 
its  predecessor.  In  each  of  these  problems,  the  flow  was  taken  to  be 
inviscid  and  incompressible  (no  thermodynamic  processes),  the  earth's 
sphericity  was  ignored  (f-plane  approximation) ,  and  the  initial  steady 
state  was  taken  to  be  geostrophic,  that  is,  a  balance  between  horizontal 
pressure  gradient  and  Coriolis  force.  The  front  was  defined  as  the 
discontinuity  surface  separating  a  thin  wedge  of  cold  air  adjacent  to 
the  earth's  surface  from  an  overlying  region  of  warm  air  of  lower 


density.  Problem  I  consisted  of  the  momentum  and  continuity  equations 
subject  to  the  restrictions  listed  above,  in  which  the  velocity  and 
pressure  fields  were  treated  as  unknowns,  along  with  appropriate 
boundary  and  initial  conditions  as  follows:  the  pressure  field  and 
normal  velocity  were  required  to  be  continuous  across  the  front,  the 
no-slip  condition  was  imposed  at  the  earth's  surface,  and  the  shape 
of  the  discontinuity  surface  was  presumed  known  initially. 

Stoker's  Problem  II  added  the  assumption  of  vertical  hydrostatic 
balance,  that  is,  vertical  accelerations  were  considered  negligible, 
as  a  consequence  of  which  the  horizontal  velocity  field  was  vertically 
uniform.  In  Problem  III,  two  further  assumptions  were  added:  fluid 
elements  initially  on  the  discontinuity  surface  remained  there  for  all 
time,  and  the  motion  of  the  warm  air  region  was  completely  unaffected 
by  the  motion  of  the  cold  air  region.  This  latter  assumption  is 
equivalent  to  requiring  that  the  cold  air  (lower)  layer  is  vanishingly 
thin  compared  to  the  upper  layer.  In  Problem  IV,  Stoker  introduced  the 
displacement  distance  of  the  front  normal  to  its  stationary  position 
as  an  additional  dependent  variable  and  assumed  that  the  horizontal 
velocity  perpendicular  to  the  Initial  front  vanished  far  from  the 
front.  However,  even  Problem  IV  could  not  be  solved  analytically  and 
presented  difficulties  that  made  numerical  solution  impractical  at 
that  time.  In  a  companion  paper  to  Stoker's,  Whitham  (1953)  pro¬ 
vided  a  heuristic  development  of  an  essentially  graphical  solution  to 
Problem  IV  using  the  method  of  characteristics  which  exhibited  the 
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nonlinear  nature  of  wave  propagation  along  frontal  surfaces  and  also 
demonstrated  the  tendency  towards  occlusion  and,  ultimately,  cyclo¬ 
genesis.  Unfortunately,  Problem  IV  was  burdened  by  so  many  restrictions 
that  it  constituted  a  rather  poor  model  of  atmospheric  frontal  zones, 
and.,  as  a  consequence,  it  was  of  limited  utility. 

Stoker's  frontal  models  led  to  mathematical  problems  that  could 
not  be  solved  analytically.  At  the  time  his  work  was  completed,  even 
numerical  solutions  could  not  be  implemented,  since  adequate  computers 
were  not  available.  Recognizing  the  futility  of  attempting  a  solution 
for  even  the  simplest  of  his  frontal  models.  Stoker  went  no  further. 

The  problem  of  nonlinear  frontal  dynamics  lay  dormant  for  more  than 
ten  years,  when,  with  the  advent  of  sufficiently  large,  high-speed 
computers,  Kasahara,  et  al. ,  (1965)  provided  the  first  numerical  solu¬ 
tion.  Kasahara  built  directly  upon  Stoker's  previous  modelling  efforts, 
and  solved  Problem  III  described  above  using  a  finite  difference  numeri¬ 
cal  method.  His  results  showed  clearly  that  initially  sinusoidal 
frontal  disturbances  propagate  from  west  to  east  and  have  a  tendency 
towards  occlusion  and,  eventually,  cyclogenesis. 

Another  numerical  study  of  essentially  the  same  nonlinear  problem 
was  performed  by  Grammeltvedt  (1970) .  His  results  also  clearly  demon¬ 
strated  the  tendency  for  frontal  disturbances  to  develop  the  charac¬ 
teristic  asymmetry  that  ultimately  leads  to  occlusion.  To  obtain 
numerical  solutions,  however,  he  found  it  necessary  to  introduce  an 
artificial  dissipative  mechanism  in  order  to  maintain  stability  in  the 
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numerical  integration.  It  appears,  in  fact,  that  some  physical 
dissipative  mechanism,  turbulent  friction,  for  example,  is  an  essential 
ingredient  in  any  realistic  frontal  model.  The  instabilities  encounter¬ 
ed  in  performing  a  numerical  Integration  of  the  nonlinear  dynamical 
equations  are  apparently  not  simply  numerical  instabilities  associated 
with  a  particular  integration  algorithm,  but  rather  important  conse¬ 
quences  of  the  inviscid  dynamics  assumed. 

Even  though  no  known  method  of  exact  solution  exists  for  the  com¬ 
plicated  system  of  equations  that  describes  frontal  motion,  certain 
theoretical  results  are  nonetheless  available  using  the  approximations 
of  classical  stability  and  small  perturbation  analysis.  The  numerical 
models  that  demonstrated  frontal  wave  instability  provided  impetus  for 
renewed  analysis,  since  even  approximate  solutions  to  the  dynamical 
equations  can  provide  considerably  more  insight  into  the  physical  pro¬ 
cesses  at  work  than  numerical  results.  Since  wavelike  disturbances  of 
the  frontal  zone  are  observed,  it  is  reasonable  to  assume  a^  priori  that 
there  exist  wavelike  solutions  to  the  dynamical  equations;  in  particular, 
solutions  which  contain  an  explicit  time-harmonic  dependence  and  which 
permit  complex-valued  circular  frequencies.  This  is  the  approach 
adopted  in  classical  stability  analysis.  Solutions  exhibiting  purely 
real  circular  frequencies  are  bounded  for  all  time  and  therefore  termed 
stable.  Solutions  whose  frequency  contains  a  nonzero  imaginary  part, 
however,  correspond  to  a  time-dependent  amplitude  that  either  decays  or 
grows  exponentially,  depending  upon  the  sign  of  the  imaginary  part. 

Waves  of  the  first  type  are  termed  evanescent,  while  those  of  the 
second  type  are  termed  unstable. 
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The  simplification  afforded  by  assuming  wavelike  solutions, 
however,  is  not  alone  sufficient  to  render  the  dynamical  equations 
tractable,  since  they  remain  inherently  nonlinear.  Linearization 
of  the  fundamental  equations  is  accomplished  by  using  classical  small 
perturbation  theory,  in  which  it  is  assumed  that  wave  disturbances  to 
the  system  constitute  small  departures  from  a  known,  time  invariant 
equilibrium  state.  Each  dependent  variable  in  the  problem  is  thus 
decomposed  into  the  sum  of  an  equilibrium  state  component  and  a  harmonic 
perturbation  component  whose  amplitude  is  small  by  comparison.  This 
scheme  results  in  a  linear  system  of  equations  for  the  perturbation 
variables  which  may  then  be  solved  using  standard  analytical  techniques. 
The  principal  objective  of  this  synthesis  of  stability  and  perturbation 
analysis  is  the  derivation  of  a  dispersion  relation,  that  is,  a  function 
which  relates  the  disturbance  wavelength  to  the  frequency.  This  rela¬ 
tion  then  can  be  used  to  Identify  unstable  modes  in  the  system. 

Using  such  a  combined  stability/perturbation  approach,  Orlanski 
(1968)  examined  the  propagation  of  atmospheric  frontal  waves  and  achieved 
some  measure  of  success.  The  system  he  modelled  was  essentially  the 
same  as  that  used  by  Kasahara  and  Stoker,  although  he  considered  the 
coupling  of  the  dynamics  in  each  of  the  two  frontal  regions  to  first 
order  in  a  small  parameter,  which  was  not  a  feature  of  the  Kasahara/ 
Stoker  model.  Orlanski's  efforts  were  subsequently  extended  by  Kasahara 
and  Rao  (1972) ,  who  employed  the  same  frontal  model  and  stability/ 
perturbation  treatment  in  order  to  extend  the  analysis  to  regions  of  the 
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problem's  parameter  space  that  had  not  been  previously  considered. 

Their  results  were  very  similar  to  Orlanskl's.  However,  like  his, 
theirs  did  not  include  any  dissipative  effects,  even  though 
Grammeltvedt's  numerical  work  provided  indirect  evidence  of  the  impor¬ 
tance  of  dissipative  mechanisms  in  modelling  frontal  dynamics. 

1.3  THEORETICAL  STUDIES  OF  OCEANIC  FRONTS 

Garvine  (1974)  published  a  model  of  quasisteady,  small  scale 
oceanic  frontal  dynamics  in  which  frictional  dissipation  and  mass 
entrainment  figured  prominently.  His  model  showed  clearly  that  a 
steady  state  frontal  structure  can  exist  only  if  friction  and/or  mass 
entrainment  are  present  to  balance  the  net  horizontal  pressure  gradient. 
Otherwise,  the  sloping  interface  between  two  different  water  masses 
cannot  be  maintained.  In  contrast,  the  models  of  large  scale  atmos¬ 
pheric  fronts  described  in  Section  1.2  are  all  inviscld  and  invoke  a 
basic  balance  between  horizontal  pressure  gradient  and  Coriolis  force. 
Only  in  large  scale  oceanic  fronts  where  Coriolis  effects  will  be  signi¬ 
ficant  can  dissipative  effects  be  expected  to  play  a  lesser  role. 

Garvine  (1979a,  1979b)  has  expanded  his  small  scale  model,  in 
which  Coriolis  effects  were  neglected,  to  include  quasisteady  large 
scale  oceanic  fronts  whose  dynamics  involve  friction,  mass  entrainment, 
wind  stress  and  Coriolis  force.  In  view  of  the  inherently  turbulent 
nature  of  geophysical  flows  and  the  attendant  analytical  difficulties 
in  treating  the  turbulence  structure  in  detail,  Garvine  parameterized 
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the  effects  of  friction  and  mass  entrainment  by  introducing  friction 
and  entrainment  coefficients  whose  values  reflected  the  bulk  influence 
of  these  effects.  He  also  wrote  the  momentum  and  continuity  equations 
in  vertically  integrated  form.  His  model  showed  that  the  steady  state 
frontal  zone  is  characterized  by  two  length  scales,  the  baroclinic 
Rossby  radius  of  deformation  and  a  turbulent  transport  length  scale. 

Their  ratio  he  called  the  "rotation  parameter",  which  was  the  major 
parameter  characterizing  the  frontal  zone  structure.  For  nearly  non- 
rotating  systems  in  which  Coriolis  effects  are  small,  the  frontal  zone 
scale  is  the  turbulent  transport  scale,  and  the  frontal  zone  is  dissi¬ 
pative  throughout  its  domain.  This  limit  corresponds  to  Garvine's 
1974  model  which  considered  only  nonrotating  systems.  Small  scale 
frontal  structures,  such  as  a  river  plume  front  (Garvine  and  Monk,  1974), 
exhibit  small  values  of  the  rotation  parameter. 

In  contrast,  large  scale  fronts,  such  as  the  Gulf  Stream,  are 
strongly  affected  by  earth  rotation,  but  only  weakly  by  dissipation. 

Their  corresponding  values  of  the  rotation  parameter  are  large,  and 
the  scale  of  the  front  is  the  Rossby  radius.  The  weak  dissipative  pro¬ 
cesses  of  friction  and  entrainment  are  active  only  in  an  inner  zone 
bordering  the  surface  front  where  the  gradient  Richardson  number  falls 
below  a  value  of  order  one  due  to  the  large  vertical  shear  there. 
Typically,  the  lateral  extent  of  this  region  is  10-20%  of  the  Rossby 
radius.  Hence,  in  this  zone,  inviscid  dynamics  are  not  valid.  Beyond 
this  zone,  the  flow  is  nearly  inviscid;  but  even  here  the  motion  is 
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indirectly  affected  by  the  dissipative  effects  of  the  inner  zone. 

Any  vertical  mass  entrainment  across  the  frontal  interface  occuring  in 
the  inner  zone  must  be  balanced  by  horizontal,  cross-stream  mass  flux 
above  the  interface;  i.e. ,  flow  toward  or  away  from  the  surface  front. 
This  mass  flux,  in  turn,  must  be  continuous  across  the  boundary  between 
the  inner  region  and  the  outer  inviscid  one.  Hence,  the  inviscid  region 
responds  Indirectly  to  inner  zone  mass  entrainment.  If  the  entrainment 
there  is  downward,  mass  is  lost  to  the  ambient  water  below,  and  the  in¬ 
viscid  outer  zone  must  supply  balancing  cross-stream  mass  flux  (flow 
toward  the  surface  front) .  The  opposite  occurs  for  upward  mass  entrain¬ 
ment. 

In  this  paper,  we  will  consider  only  upper  ocean  density  fronts 
with  large  values  of  the  rotation  parameter,  so  that  dissipative  effects 
are  subordinated  to  earth  rotation  effects.  Typical  of  the  large-scale 
oceanic  fronts  to  which  this  model  applies  are  the  Gulf  Stream  Front  and 
the  Sargasso  Sea  fronts,  for  which  Garvine  (1979b)  estimates  rotation 
parameter  values  of  9  and  4,  respectively.  For  the  Gulf  Stream,  the 
baroclinlc  Rossby  radius  is  on  the  order  of  100  km,  so  that  dissipation 
is  confined  to  a  zone  approximately  10-20  km  wide.  It  will  be  seen  be¬ 
low  that  there  does  exist  a  third  length  scale  in  the  frontal  problem 
which  is  imbedded  in  the  dissipative  zone.  On  this  scale,  which  is  an 
advectlve  scale  length,  the  cross-flow  Rossby  number  becomes  order  one. 

The  work  presented  in  this  paper  is  an  analytical  study  of  the 
stability  characteristics  of  large-scale  upper  ocean  density  fronts. 

The  model  is  developed  by  assuming  wavelike  disturbances  to  the  system. 


which  are  then  treated  as  small  amplitude  perturbations  to  a  stationary 
equilibrium  flow.  As  outlined  in  Section  1.2,  the  perturbation  variable 
decomposition  is  a  standard  technique  in  the  atmospheric  literature, 
and  also  in  prior  modelling  efforts  on  time-dependent,  inviscid  oceanic 
fronts  (Duxbury,  1963;  Mysak  and  Schott,  1977;  Orlanski,  1969;  Stommel, 
1953;  Stommel,  1976).  Time-dependent  numerical  models  have  also  been 
developed  of  both  small  scale  fronts  (Kao,  al. ,  1977)  and  large 
scale  fronts  (Kao,  et  al.,  1978),  but  not  with  the  objective  of  study¬ 
ing  wave  disturbances  on  established  frontal  structures.  Both  of  these 
papers  instead  address  the  problem  of  frontogenesis  with  the  ultimate 
development  of  a  steady  state  frontal  structure. 

The  present  model  first  extends  Duxbury' s  (1963)  treatment  of  an 
inviscid  density  front  with  a  planar  discontinuity  surface.  It  is 
significantly  different,  however,  in  that  it  subsequently  examines  ah 
exponential  discontinuity  surface.  This  form  was  suggested  by  Stommel 
(1976)  for  the  Gulf  Stream  front.  The  model  also  includes  the  important 
effect  of  a  cross-stream  velocity  component  (flow  perpendicular  to  the 
surface  front)  in  the  equilibrium  flow,  Duxbury' s  system  included  a 
constant  velocity  shear  across  the  boundary  between  the  frontal  zone 
and  the  ambient  ocean,  but  had  no  cross-stream  velocity.  It  was  found 
to  be  dynamically  stable  for  all  values  of  frequency  and  wave  number. 
Orlanski 's  (1969)  work  showed  that  systems  in  which  the  dynamics  of  the 
lower  layer  are  neglected  in  a  two-layer  frontal  model  are  inherently 
stable.  However,  he  too  did  not  admit  any  equilibrium  cross-stream 
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flow.  In  this  paper,  we  see  that,  even  neglecting  the  dynamics  of  the 
lower  layer,  a  frontal  system  subject  to  harmonic  wave  disturbances  may 
be  unstable  If  a  nonzero  cross-stream  velocity  Is  present  In  the  equi¬ 
librium  state.  The  source  of  this  cross-stream  velocity,  in  turn,  is 
mass  entrainment  in  the  inner  dissipative  region. 


CHAPTER  2 


THE  EQUILIBRIUM  STATE  AND  PERTURBATION  EXPANSION 

2.1  INTRODUCTION 

As  discussed  in  Chapter  1,  the  term  "oceanic  front"  in  this  paper 
refers  to  the  discontinuity  surface  that  separates  two  water  masses  of 
different  density.  The  "frontal  zone"  is  a  wedge  shaped  body  of  fluid 
that  floats  on  an  ambient  ocean  of  slightly  greater  density,  as  illus¬ 
trated  in  Figure  2.1.  Its  thickness  is  small  compared  to  the  total 
depth,  and,  as  a  result,  the  ambient  horizontal  pressure  gradient  field 
is  unaltered  by  the  presence  of  the  frontal  zone.  The  dynamical  be¬ 
havior  of  the  frontal  zone  and  that  of  the  underlying  ambient  fluid  are 
therefore  uncoupled,  which  is  equivalent  to  assuming  that  the  ambient 
fluid  is  effectively  infinitely  deep.  Recall  that  this  restriction  was 
also  implicit  in  Stoker's  (1953)  model  of  a  lower  tropospheric  front  in 
the  atmosphere.  The  frontal  zone  is  free  to  move  about  and  propagates 
into  the  ambient  fluid  with  a  velocity  *  (U^,  V^,  0)  relative  to  an 
earthfixed  Cartesian  coordinate  system  labeled  x' ,  y'  and  z'  in  Figure 
2.1.  The  location  of  the  coordinate  system  origin  is  arbitrary,  and 
the  z'  axis  coincides  with  the  local  vertical.  The  y'  axis  is  parallel 
to  the  surface  front,  with  the  x'  axis  being  perpendicular  to  it  into 
the  frontal  zone. 

In  nature,  oceanic  frontal  zones,  the  Gulf  Stream  front,  for 
example,  are  long  and  thin  with  a  characteristic  velocity  shear  across 
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the  front  in  the  direction  of  its  maximum  extent.  These  features 
motivate  the  "long-stream"  and  "cross-stream"  terminology  introduced  in 
the  figure.  The  long-stream  direction  coincides  with  the  direction  of 
maximum  extent  of  the  filamentary  frontal  zone,  while  the  cross-stream 
direction  is  perpendicular  to  it  and  taken  positive  toward  the  less 
dense  fluid.  The  frontal  zone's  long-stream  length  scale  is  typically 
an  order  of  magnitude  or  more  larger  than  its  cross-stream  length  scale. 
The  velocity  field  in  the  frontal  zone  exhibits  similar  behavior,  with 
the  long-stream  component  being  much  larger  than  the  cross-stream  com¬ 
ponent,  but  not  necessarily  in  the  same  ratio  as  the  length  scales. 

The  spatial  rate  of  change  of  the  velocity  and  pressure  fields,  however, 
behaves  in  the  opposite  manner,  with  derivatives  in  the  long-stream 
direction  being  generally  much  smaller  than  derivatives  in  the  cross¬ 
stream  direction. 

The  discontinuity  surface  which  defines  an  oceanic  front  thus 
divides  the  fluid  into  two  quite  distinct  regions,  the  (uncoupled) 
dynamics  of  which  are  taken  to  be  inviscid  in  this  model,  except  for 
the  innermost  portion  of  the  upper  layer.  The  model  assumes  that  the 
fluid  is  incompressible,  and  ignores  surface  wind  stress  as  a  driving 
force.  We  will  later  introduce  wavelike  disturbances  of  the  frontal 
zone,  such  as  would  be  initiated  by  some  driving  agent  (such  as  surface 
wind  stress)  which  is  no  longer  active.  In  typical  upper  ocean  frontal 
zones,  the  lateral  extent  is  sufficiently  small  compared  to  the  plane¬ 
tary  scale  that  Coriolis  acceleration  is  adequately  approximated  on  the 
f-plane,  which  is  the  assumption  made  here.  In  addition,  the  vertical 
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dynamical  balance  is  assumed  to  be  everywhere  hydrostatic,  in  view  of 
which  the  horizontal  velocity  field  is  vertically  uniform.  The  hori¬ 
zontal  dynamical  balance  in  the  ambient  fluid  is  taken  to  be  geo- 
strophic,  that  is,  a  balance  between  Coriolis  force  and  horizontal 
pressure  gradient,  and,  moreover,  is  unaltered  by  the  presence  of  the 
frontal  zone. 

The  model  is  hydrodynamic  in  nature,  not  thermodynamic,  since  the 
fluid  densities,  p  and  p  -  dp  in  the  ambient  fluid  and  in  the  frontal 

00  CO 

zone,  respectively,  are  specified  a  priori.  Both  p^  and  Ap  are  positive 

constants  with  <<  1.  Specification  of  the  density  field  in  this 
^00 

manner  leaves  only  the  velocity  and  pressure  fields  as  unknown  quanti¬ 
ties,  thereby  eliminating  the  need  for  an  energy  equation  that  couples 
the  density  field  and  frontal  zone  dynamics. 

The  model  is  restricted  to  upper  ocean  density  fronts  which,  in 
Garvine's  (1979a,  1979b)  model,  are  characterized  by  large  values  of  the 
rotation  parameter.  The  frontal  zone  dynamics  are  therefore  rotation- 
ally  dominated,  and  the  zone's  cross-stream  length  scale  is  the  baro- 
cllnic  Rossby  radius  of  deformation  X.  Even  though  dissipative  pro¬ 
cesses  are  unimportant  over  most  of  the  scale  X,  a  narrow  dissipative 
region,  whose  lateral  extent  is  approximately  0.1X  -0.2X,  exists  near 
the  surface  front  as  shown  in  Figure  2.1.  Garvine  (1979a)  has  developed 
the  quasi-steady  dynamics  for  this  case. 

In  the  present  work,  wavelike  disturbances  of  the  frontal  zone, 
including  those  leading  to  instability,  will  be  investigated  in  detail 
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for  Che  outer  Invlscid  scale.  The  model  will  assume  that,  as  Che  wave 
disturbances  develop  horizontal  distortions  of  the  surface  front  and 
adjacent  dissipative  region,  the  response  there  is  simply  to  displace 
this  region  without  changing  its  steady  state  structure.  In  other 
words,  at  a  fixed  long-stream  position,  wave  disturbances  will  simply 
advect  the  dissipative  zone  structure  in  the  cross-stream  direction. 

The  invlscid  zone  dynamics,  on  the  other  hand,  is  treated  as  explic¬ 
itly  time  dependent,  but  at  the  inner  boundary  of  this  zone,  boundary 
conditions  are  imposed  so  as  to  match  the  corresponding  values  at  the 
outer  edge  of  the  oscillating  dissipative  zone. 

2.2  THE  EQUILIBRIUM  STATE 

Figure  2.2  provides  a  cross-sectional  view  of  the  idealized 
geometry  that  is  assumed  for  the  equilibrium-state  frontal  zone.  In 
the  equilibrium,  or  "basic",  state,  the  frontal  zone  geometry  shown  in 
Figure  2.2  is  quasi-steady.  Variations  in  the  frontal  zone  structure 
are  thus  assumed  to  occur  on  a  time  scale  that  is  much  larger  than  that 
on  which  the  wave  disturbances  considered  in  this  paper  occur.  Obser¬ 
vations  show  that  wave  disturbances  on  large-scale  upper  ocean  fronts 
typically  exhibit  periods  measured  in  days  or  weeks  while  the  frontal 
zone  itself  persists  in  a  stable  configuration  for  a  much  longer  inter¬ 
val,  typically  months  or  years.  The  quasi-steady  assumption  is  there¬ 
fore  consistent  with  the  observed  characteristics  of  large-scale  frontal 
zones . 

The  dissipative  region,  whose  cross-stream  extent  is  some  small 

fraction  of  the  baroclinic  Rossby  radius,  A, lies  between  the  surface 
front  and  the  outer  invlscid  region.  The  ambient  ocean  is  labeled  Region 

II,  while  the  frontal  zone  is  labeled  Region  I.  An  earth-fixed 
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CROSS-SECTIONAL  VIEW  OF  THE 
IDEALIZED  FRONTAL-ZONE 
GEOMETRY  IN  EARTH-FIXED 
COORDINATES 


FIGURE  2.2 
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(primed)  right-handed  Cartesian  coordinate  system  is  established  in 

which  the  spatial  axes  are  labeled  x' ,  y*  and  z',  with  t'  designating 

the  time.  The  z'  axis  is  taken  to  be  along  the  local  vertical,  and  the 

y'  axis  is  oriented  parallel  to  the  surface  front.  The  boundary  between 

the  dissipative  and  inviscid  domains  in  the  frontal  zone  is  therefore 

parallel  to  the  plane  x'  *0.  The  sea  surface  elevation  above  the 

reference  plane  z'  •  0  is  designated  h^  and  h^  in  Regions  I  and  II, 

respectively,  while  the  depth  of  the  discontinuity  surface  in  Region  I, 

that  is,  the  "interfacial  depth",  is  denoted  D'.  The  fluid  densities  in 

Regions  I  and  II  are  p  -dp  and  p  ,  respectively.  Note  that  the  dissi- 

pative  region  is  part  of  the  frontal  zone  and  contains  fluid  of  density 

P  -dp.  Note  also  that  the  x'  velocity  at  the  boundary  between  the  dissi- 
00 

pative  and  inviscid  domains  in  the  frontal  zone  will  not,  in  general, 
be  zero,  which  reflects  mass  entrainment  through  the  discontinuity  sur¬ 
face  in  the  dissipative  region. 

Under  the  assumptions  made  in  Section  2.1,  the  general  time- 
varying  x',  y’  and  z'  momentum  equations,  and  the  continuity  equation  for 
the  outer,  inviscid  region  in  the  primed  coordinates  are,  respectively: 


au'/at'+u’Ou'/ax^+v'Ou'/ay'J+w'Ou'/az'J-fv' 

—  Op73x')/p  (2.1a) 

av’/at'+u^av'/ax^+v'Ov'/ay^+w'  (3v73z')+fu' 

— (3p'/3y')/p  (2.1b) 

0*pg+(3p,/3z')  (2.1c) 

au'/ax'+av'/dy'+aw'/sz'-o  (2. id) 
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In  (2.1),  u' ,  v'  and  w'  are  the  x',  y'  and  z'  velocities,  respec¬ 
tively,  p'  the  pressure,  p  the  density,  g  the  gravitational  acceler¬ 
ation,  and  f  the  Coriolis  parameter.  All  variables  are  dimensional. 

(2.1)  constitutes  a  system  of  four  coupled  nonlinear  partial  differen¬ 
tial  equations  (PDE's)  for  the  four  dependent  variables  u',  v' ,  w' 
and  p',  each  of  which  is  a  function  of  position  and  time  (x',  y',  z',  t'). 
The  horizontal  velocity  field  in  the  ambient  ocean  is  assumed 

to  be  invariant  in  space  and  time,  and  is  denoted  <$  *  (U  ,V  ,0), 

a  a  a 

where  U  and  V  are  the  ambient  x'  and  y'  velocity  components  as  shown 
a  a 

in  Figure  2.2.  Horizontal  pressure  gradients  in  the  ambient  fluid  are 
induced  by  the  sloping  free  surface  in  Region  II,  whose  elevation  above 
the  plane  z*  ■  0  is  related  to  U  and  V  by 

cL  <1 

fv  *g3h ' /3x 1  (2.2a) 

a  i 

-fu  *g3hl/3y '  (2.2b) 

a  £ 

as  a  consequence  of  the  assumed  geostrophic  balance.  The  ambient  sea 
surface  height  is  therefore  a  time-invariant  plane  in  the  primed  coord¬ 
inate  system  whose  uniform  slope  in  the  x'  and  y'  directions  drives  a 
uniform  geostrophic  flow  in  the  orthogonal  coordinate  direction. 

The  thin  wedge  of  lighter  fluid  that  constitutes  the  frontal 
zone  propagates  into  the  ambient  fluid  with  a  constant  velocity  in  the 
primed  system,  in  the  x*  direction  and  in  the  y'  direction  as 
shown  in  Figure  2.2.  The  free  surface  elevation  in  Region  I,  the 
frontal  zone,  is  therefore  time-dependent  in  earth-fixed  coordinates, 
h[-h[(x' ,y' ,t') ,  as  is  the  interfacial  depth  D'»D' (x' ,y' ,t’) .  The 
explicit  time  dependence  in  h j  and  D'  is  analytically  cumbersome,  and 


I 
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it  is  conveniently  removed  from  the  equilibrium  state  by  transforming 
into  a  new  Cartesian  coordinate  system  whose  origin  translates  rela- 
tive  to  the  earth  at  the  frontal  zone  velocity  Qfa(Uf ,Vf ,0) .  In  the 
new  coordinate  system,  the  same  symbols  are  retained  for  the  dependent 
and  independent  variables,  but  without  the  prime.  The  transformation 
is  defined  by 

x  -  x'-Uft 
y  -  y'-Vft 
z  =  z' 


T 

4 

T 

t  -  t' 

u  *  u'-U^ 

(2.3) 

i 

i 

v  *  v'-Vf 

w  *  w' 

in  view  of  which  (2.1)  becomes 

i 

v 

u  +uu  +vu  +wu  -f(v+V_)  “  -p  /p 
t  x  y  z  f  x 

(2.4a) 

I 

1 

v  +uv  +w  +wv  +f(u+U,)  =  -p  /P 
t  x  y  z  f  y 

(2.4b) 

1 

o  ■  g+p  /p 

Z 

(2.4c) 

u  +v  +w 
x  y  z 


0 


(2.4d) 
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CROSS-SECTIONAL  VIEW  IN  A  PLANE 
PARALLEL  TO  y=0 


GEOMETRY  OF  THE  EQUILIBRIUM 
FRONTAL  ZONE  IN  TRANSLATING 
COORDINATES 


FIGURE  2.3 
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(2.4)  has  been  simplified  notationally  by  marking  partial  der¬ 
ivatives  with  a  subscript  following  accepted  convention.  This  notat¬ 
ion  is  used  hereafter  whenever  possible,  with  the  significance  of 
subscripted  variables  that  are  not  derivatives  being  clear  from  con¬ 
text.  Figure  2.3  provides  a  cross-sectional  view  of  the  frontal  zone 
geometry  in  the  translating  coordinate  system.  As  before,  the  sea 
surface  elevations  in  Regions  I  and  II  are  denoted  h^  and  respec¬ 
tively,  while  D  denotes  the  interfacial  depth.  The  free  surface 
height  is  now  explicitly  time-dependent,  as  was  h|  previously;  but 
h^  and  D  in  the  equilibrium  state  are  not.  It  will  be  seen  in  what 
follows,  however,  that  only  the  uniform,  constant  velocity  in  the 
ambient  ocean  enters  the  equilibrium  state  equations  for  the  frontal 
zone  velocity  field  and  interfacial  depth,  so  that  the  time  dependence 
in  h2(x,y,t)  is  of  no  consequence.  The  ambient  fluid  velocity  in  the 
translating  coordinate  system  is  labeled  u  along  the  x-axis  and  v 

CO  CO 

along  the  y-axis,  where  u  and  v  are  constants  related  to  the  earth- 

00  00 

fixed  ambient  and  frontal  zone  velocities  by 
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The  pressure  field  Is  calculated  by  Integrating  the  hydrostatic 
equation,  (2.4c),  using  (2.6),  with  the  result 

P  *  g(P0#-Ap)(h1-z)  ,  -D  <  z  <  h1  (2.7a) 

P  ”  g(Pa8-4p)h1-g(Pooz+4pD)  ,  z  <  -D  (2.7b) 

Neglecting  terms  of  order  (Ap/p^)  compared  to  unity,  the  hori¬ 
zontal  pressure  gradients  in  the  half  space  x  >  0  are  computed  from 
(2.7)  as 

V  p/p  -  gV  h  -D  <  z  <  h  (2.8a) 

n  hi  -  -  1 

^hp/p  -  g  Vh(h  -ApD/pJ  z  <  -D  ,  (2.8b) 

A  A 

where  *(3/3x)i+(3/3y) j  is  the  horizontal  gradient  operator.  Note 
that  the  horizontal  gradients  of  the  pressure  field  are  discontinuous 
across  the  frontal  interface  marked  by  z  *  -D,  whereas  the  pressure 
field  is  continuous,  as  it  must  be. 

Below  the  frontal  interface  in  Region  I,  the  fluid  is  assumed 
to  be  in  isostatic  balance,  that  is,  the  ambient  horizontal  pressure 
gradient  field  is  unchanged  by  the  presence  of  the  frontal  zone  above. 
This  restriction  is  equivalent  to  assuming  that  the  ambient  fluid  is 
infinitely  deep,  and,  therefore,  that  its  (geostrophic)  dynamics  are 
unaltered  as  the  frontal  zone  propagates  into  the  ambient  fluid.  The 
velocity  and  pressure  graduent  fields  in  the  ambient  fluid  are  then 
related  by 

A  A 

V/P  "  f{<VVf)i-(U-+Uf)j} 


(2.9) 


(2.9)  follows  from  (2.4a)  and  (2.4b)  and  is  analogous  to  (2.2)  in 
earth-fixed  coordinates  in  which  the  pressure  gradients  have  been 
written  in  terms  of  the  free  surface  slope. 

In  view  of  (2.8),  the  horizontal  gradients  of  the  sea  surface 
elevation  and  interfacial  depth,  and  the  velocity  field  are  related 
as  follows: 

Vh(hrApD/pJ  =  (f/g){(V-+V  )i-(u<o+Uf)j}  .  (2.10) 

Integrating  (2.10)  subject  to  the  boundary  conditions  h^(0,0)=0 
and  D(0 ,0)=D^ (0) ,  where  D^(0)  is  some  known  constant,  yields  the 
following  relationship  between  h^  and  D: 

h(x,y)  =  (Ap/p  ){D(x,y)-D.(0)}+(f/g){(v  +V,)x-(u  +U,)y}  (2.11) 

X  £  OO  £  00  £ 

In  the  inviscid  zone,  mass  entrainment  is  not  permitted  through 
the  discontinuity  surface.  However,  entrainment  through  the  discontin¬ 
uity  surface  is  allowed  in  the  dissipative  zone  shown  in  Figure  2.3, 
which  produces  a  nonzero  cross-stream  volume  flux  in  the  plane  x=0. 

This,  in  turn,  leads  to  a  velocity  boundary  condition  at  x=0,  since 
we  demand  that  the  cross-stream  velocity  component  be  continuous  across 
the  boundary  between  the  inner  dissipative  region  and  the  outer  in¬ 
viscid  region. 

The  requirement  that  fluid  elements  initiall.  on  the  : lviscid 
zone  discontinuity  surface  remain  there  for  all  time  may  be  stated 
rigorously  as 

lim  Dt^hl“Z^  =  0  (2.12a) 

z-*h^ 

lim  ?r(IH-z)  -  0  (2.12b) 

z-*— D 

where  D/Dt  ■  9/3t+u(3/3x)+v(3/3y)+w(3/3z)  is  the  convective  derivative. 

Denoting  the  frontal  zone  vertical 
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2.3  PERTURBATION  ANALYSIS 

In  this  section,  the  general  nonlinear  time  dependent  frontal 
zone  momentum  and  continuity  equations  will  be  linearized  by  assuming 
a  perturbation  expansion  for  each  of  the  dependent  variables  u,  v,  and 
D.  Such  a  linearization  scheme  is  applicable  only  if  the  perturbation 
amplitude  is  small  compared  to  the  time  invariant  equilibrium  state 
value  of  a  given  variable,  and  this  restriction  is  implicit  in  all  of 
what  follows.  For  large  deviations  from  the  equilibrium  state,  the 
frontal  zone  dynamics  are  essentially  nonlinear  and  must  be  treated 
accordingly.  The  purpose  of  this  paper,  however,  is  to  address  the 
question  of  stability  of  the  frontal  zone  flow,  and  instabilities  that 
grow  to  large  amplitude  in  time  must  be  small  initially.  Large  ampli¬ 
tude  frontal  disturbances  are  therefore  not  within  the  purview  of  this 
analysis,  but  the  conditions  under  which  they  develop  are. 

In  analogy  co  (2.16),  the  general  time  varying  momentum  and  con¬ 
tinuity  equations  are 


u  +uu  +vu  -fv  *  -g'D  -fv 
t  x  y  x  00 

(2.17a) 

v  +uv  +w  +fu  =  -g'D  +fu 
t  x  y  y  00 

(2.17b) 

u  +v  +w  =0, 
x  y  z 

(2.17c) 

in  which  u,  v  and  D  are  functions  of  (x,  y,  t)  and  w  =  w(x,  y,  z,  t). 

g'  is  the  "reduced  gravity"  defined  as  g'  H  g  (— )  .  (2.17)  makes  no 

^00 

assumption  about  the  frontal  zone  geometry.  As  before,  all  variables 
are  dimensional.  In  (2.17a)  and  (2.17b),  the  horizontal  pressure 
gradients  have  been  written  in  terms  of  the  interfacial  slopes,  3D/3x 
and  3D/3y,  In  order  to  eliminate  the  vertical  velocity,  w,  in  (2.17c), 
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a  vertically  integrated  form  of  the  continuity  equation  is  introduced 
by  following  the  same  procedure  employed  in  the  previous  section,  but 
with  h^  and  D  now  functions  of  time  as  well  as  space.  Neglecting  terms 
of  order  h^/D  compared  to  unity,  the  counterpart  of  (2.15)  is  then 

(uD)  +(vD)  +D  ■  0,  (2.18) 

x  y  z 

which  replaces  (2.17c)  and  is  used  exclusively  hereafter. 

We  now  assume  that  the  dependent  variables  may  be  expanded  as 

u  *  u  +  u 

v  -  v  +  v  (2.19) 

D  -  5  +  D 

where  u,  v  and  D  are  the  equilibrium  state  velocities  and  interfacial 
depth,  respectively,  and  where  the  perturbation  (caret)  variables  are 
small  by  comparison.  This  decomposition  and  the  subsequent  lineariza¬ 
tion  of  (2.17)  are  therefore  valid  only  if 

|u/u|«l,  |  v/v  |  <<1 ,  ID/D  |  <<1.  (2.20) 

The  basic  state  variables  in  (2.19)  are  solutions  to  the  basic  state 
equations  given  in  (2.16).  The  planar  front  solution  to  this  system 
is  developed  in  the  next  chapter,  tfhlle  the  exponential  front  solution 
is  derived  in  Chapter  5. 


Substituting  (2.19)  into  (2.17a),  (2.17b)  and  (2.18),  expanding, 
and  neglecting  second  order  terms  (products  of  perturbation  variables) 
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leads  to  the  following  system  of  linear  partial  differential  equations 
(PDE's)  in  the  new  unknowns  u,  v  and  D: 


u  +uu  +vu  +u  d-fv+g'D  «  0 
t  x  y  x  x 

(2.21a) 

v  +uv  +vv  +v  u+fu+e 'D  =  0 
t  x  y  x  y 

(2.21b) 

(Gr>)  +(vD)  +(Du)  +(Dv)  +D  =  0  . 
x  y  x  y  t 

(2.21c) 

The  linearized  system  in  (2.21),  written  in  dimensional  variables, 
is  general  and  therefore  applies  to  both  the  planar  and  exponential 
fronts.  In  the  planar  front  case,  which  is  developed  in  detail  in  the 
next  two  chapters,  (2.21)  simplifies  considerably,  since  the  equilibrium 
velocity  is  spatially  invariant  and  its  spatial  derivatives  consequently 
vanish.  In  Chapter  5,  which  introduces  the  exponential  front,  an  ageo- 
strophic  equilibrium  state  is  considered.  Terms  in  (2.21)  involving 
spatial  derivatives  of  the  basic  state  flow  must  then  be  retained, 
leading  to  a  more  complicated  system  of  perturbation  equations.  Note 
that  the  time  independent  equilibrium  state  variables  appear  in  (2.21) 
as  non-constant  coefficients,  but  not  as  forcing  terms.  The  perturba¬ 
tion  system  is  thus  homogeneous.  Perturbations  to  the  basic  state  flow 
are  assumed  to  have  developed  from  an  external  forcing,  surface  wind 
stress,  for  example,  which  is  no  longer  active.  Moreover,  in  view  of 
the  homogeneity  of  (2.21),  the  perturbations  are  not  forced  by  the 
equilibrium  state  variables. 


CHAPTER  3 


THE  PLANAR  FRONT 


3.1  INTRODUCTION 

In  this  chapter*  the  equilibrium  state  solution  to  (2.16)  is  pre¬ 
sented  for  the  planar  front,  wavelilte  solutions  are  then  assumed  for 
the  perturbation  equations  in  (2.21),  and  the  general  equation  for  the 
amplitude  of  the  depth  perturbation  is  derived  and  solved.  Dimension¬ 
al  variables  are  used  throughout  in  order  to  compare  the  results,  de¬ 
veloped  here  with  those  of  Duxbury  (1963) .  It  is  shown  that  the  general 
ordinary  differential  equation  (ODE)  for  the  depth  perturbation  ampli¬ 
tude  (DPA)  reduces  to  Duxbury's  second  order  equation  for  the  case  of 
zero  cross-stream  flow.  In  its  general  form,  however,  which  allows 
nonzero  cross-stream  flow,  the  DPA  equation  is  a  seventh  order  ODE,  the 
implications  of  which  are  discussed  in  this  and  the  next  chapter. 

3.2  THE  PLANAR  FRONT  EQUILIBRIUM  STATE 

For  the  planar  front,  the  solutions  to  the  time  independent  basic 
state  equations  (2.16)  are 


-  D.(OH(fB/g')(-v  x+u  y) 

X  00  oo 

(3.1) 

u  -  (l-8)u 

00 

(3.2) 

v  -  (1-6) v 

OO 

(3.3) 

in  which  6  is  an  arbitrary  constant.  (3. 1)— (3. 3)  completely  specify 
the  equilibrium  state  velocity  field  and  interfacial  depth  for  the 
planar  front,  which  then  appear  in  the  perturbation  system  (2.21)  as 


38 


39 


''Mtm 


I 

I 

I 

I 

I 

I 

l 

I. 

• ». 


fT 


variable  coefficients.  The  discontinuity  surface  that  separates  the 
frontal  zone  from  the  ambient  ocean  is,  in  view  of  (3.1),  a  plane; 
hence,  the  term  planar  front.  The  free-surface  elevation,  which  is 
related  to  the  interfacial  depth  by  (2.11),  is  therefore  also  a  plane. 
The  horizontal  velocity  field  is  spatially  invariant,  and  includes  a 
nonzero  cross-stream  velocity.  The  resulting  nonzero  cross-stream 
volume  flux  at  x  ■  0  is  driven  by  mass  entrainment  through  the  discon¬ 
tinuity  surface  in  the  dissipative  region.  The  bulk  effect  on  the  outer 
(lnvlscid)  region  of  dissipation  near  the  surface  front  is  therefore 
parameterized  here  by  the  single  free  constant  S. 

3.3  HARMONIC  SOLUTIONS  TO  THE  PLANAR  FRONT  PERTURBATION  EQUATIONS 
*  + 

Even  with  the  level  of  simplification  achieved  by  introducing  the 

perturbation  expansion  in  Section  2.3,  the  analytical  solution  of  the 

system  (2.21)  constitutes  a  formidable,  if  not  impossible,  task.  Since 

we  anticipate  wavelike  disturbances  of  the  frontal  zone  flow  field  and 

Interfacial  depth,  a  further  simplification  is  introduced  by  assuming 

time  and  space-harmonic  solutions  to  (2.21)  of  the  form 

u(x,y,t)  -  X(x)exp(i(ky-Gt) }  (3.4a) 

$(x,y,t)  ■  Y (x) exp { i (ky-Qt ) }  (3.4b) 

D(x,y,t)  -  Z(x)exp{i(ky-Ot) }  (3.5) 

1/2 

where  i  ■  (-1)  ,  0  is  the  disturbance  circular  frequency,  and  k  the 

wavenumber.  X(x),  Y(x)  and  Z(x)  are  the  (complex-valued)  x-dependent 
amplitudes  of  the  perturbations  to  the  x  and  y  velocities  and  inter- 
facial  depth,  respectively.  In  addition  to  specifying  the  y  and  t 
dependence  explicitly,  solutions  harmonic  in  these  variables  provide 
the  additional  advantage  of  not  requiring  boundary  conditions  in  y  or 
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c.  The  x-direction  boundary  conditions  are  formulated  rigorously  in 
Chapter  7. 

Substitution  of  (3.4)  and  (3.5)  into  (2.21)  along  with  the  basic 
state  equations  (3. 1)— (3. 3)  reduces  the  system  of  perturbation  PDE's 
for  0,  V  and  D  to  a  system  of  ODE's  for  the  planar  front  perturbation 
amplitudes  X(x),  Y(x)  and  Z(x)  as  follows: 


imX-fY+uX'+g'Z'  -  0 

(3.6a) 

iwY+f X+uY+ikg ' Z  -  0 

(3.6b) 

io)Z-Htf'+(D  +ikD)Y+D  X+DX'  -  0, 

y  ^ 

(3.6c) 

where  u  ■  kv  -  !J  is  the  Doppler-shifted  disturbance  frequency.  Since 
(3.6)  contains  only  derivatives  with  respect  to  x,  it  has  been  simpli¬ 
fied  notationally  by  marking  derivatives  with  a  prime;  for  example, 

Z’  -  dZ/dx. 

3.4  DERIVATION  OF  THE  PLANAR  FRONT  DPA  EQUATION 

A  single  ODE  for  the  DPA  (depth  perturbation  amplitude),  Z(x), 
may  be  developed  by  eliminating  X(x)  and  Y(x)  in  (3.6),  which,  in  turn, 
is  most  easily  accomplished  by  writing  the  system  (3.6)  in  operator 
matrix  notation.  The  interfacial  slopes  in  the  x  and  y  directions, 
respectively,  will  be  denoted  a^  and  a^,  so  that  the  equilibrium  inter¬ 
facial  depth  (3.1)  becomes 

D(x,y)  -  Df (0)+a1xfa2y . 


(3.7) 
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Replacing  the  total  derivative  operator  by  the  symbol  p;  i.e.. 


p  ■  d/dx,  (3.6)  thus  becomes 


»j+Dp  a^+ikD 


(3.8) 


After  considerable  manipulation,  this  matrix  may  now  be  reduced  to  the 
triangular  form: 


g’p  X 


2  2 
f  +L1 


(3.9) 


2  2 

(f  +LpL5 
+Mp)L, 


where  the  operators  and  A(p)  are  defined  by 


*  i<*H-up  (3.10a) 
L2  -  a^up+f^  ^4*^)4!.*  (a2+ik5)  (3.10b) 
L3  -  a1g'up2+g'L1(a1+6p)p-L3  (3.10c) 
La  -  -fg'p+ikg'L^  (3.10d) 
L5  -  L2L4+(f2+L2)L3  (3.10e) 


A(p)/(2a  u)  ■  ku>2-uj2(f+3iTtu)p+u)u(2if-3ku)p2 

1  .  2  3  (3. lOf) 

4-u  (f-HLku)p 


The  final  form  of  the  ODE  satisfied  by  Z(x)  is  therefore 

{(f2+L2)L5+A(p)L4>Z(x)  -  0.  (3.11) 

An  examination  of  the  operators  in  (3.11)  shows  that  the  order  of 
the  highest  derivative  of  Z(x)  is  seven;  that  is,  Z  satisfies  a  seventh 
order,  linear,  homogeneous  ODE.  The  explicit  form  of  the  equation  for 
the  depth  perturbation  amplitude  is  obtained  from  (3.11)  by  expanding 
each  of  the  operators  using  the  appropriate  definitions.  Although 
straight  forward,  the  expansion  for  the  general  case  is  tedious,  and  is 
therefore  not  reproduced  here. 
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The  special  case  of  zero  equilibrium  cross-stream  flow  (u  =  0) , 
however,  is  readily  derived  and  is  of  considerable  interest,  since 
it  corresponds  to  Duxbury's  (1963)  model  of  the  planar  front.  In  (3.11) 
we  set  u  -  uo  ■  *  0,  and  upon  simplification  arrive  at  the 

following  equation  for  Z(x): 


{(u>2-f2)/g’+fa1k/a)-k2D)Z(x)+a1(dZ/dx)+5(d2Z/dx2)  -  0.  (3.12) 


Rewriting  (3.12)  in  Duxbury's  notation,  which  is  discussed  in  Appendix 
2,  recovers  precisely  his  equation  (49),  thereby  verifying  (3.11)  in 
this  limiting  case. 

For  the  general  case  (u  ^  0) ,  the  explicit  form  of  the  equation 
satisfied  by  Z(x)  is 
7 

y;  (VnX+Yn)(drZ/dxn)  =  0,  (3.13) 

n“0 

where  the  vn  and  Yn  are  complex  constants  given  in  Appendix  1,  eq. 
(A.I.).  Note  that  the  may  contain  a  parametric  dependence  on  y. 

3.5  SOLUTION  OF  THE  DPA  EQUATION  FOR  THE  PLANAR  FRONT 

The  ODE  for  the  depth  perturbation  amplitude  Z(x),  (3.13),  has  an 
Irregular  singular  point  at  infinity.  Its  nature  there  is  ascertained 
by  transforming  to  the  new  independent  variable  s  *  1/x  and  considering 
the  behavior  of  the  solution  Z(s)  for  s  -*■  0.  The  details  are  not  in¬ 
cluded  here.  However,  one  finds  that  the  point  at  infinity  is  an 
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irregular  singularity,  so  that  no  power  series  expansion  exists  for  the 
solution  Z(x)  as  x  +  «.  Thus,  the  method  of  Frobenlus  (Hildebrand, 

1962,  Ch.  4)  is  not  useful  in  developing  the  solution  for  large  x.  An 
asymptotic  normal  form  solution  does  exist,  however,  as  x  grows  arbi¬ 
trarily  large  (Ince,  1956,  Ch.  XVII),  its  form  being 

Z(x)  ‘v*  xreaXW(x)  for  x-*»  (3.14) 

where 

W(x)  -  l+(Wj/x)+(w  /x2)+...  (3.15) 

a,  r,  and  the  w  are  all  constants  to  be  determined.  Both  the  expo- 
n 

nentlal  factor  in  (3.14)  and  the  series  W(x)  In  (3.15)  are  non-dimen¬ 
sional.  It  is  therefore  understood  that  the  constant  1  with  dimensions 
of  length  raised  to  the  (1-r)  power  multiplies  the  right-hand  side  of 
(3.14),  rendering  Z(x)  dimensionally  correct.  The  constants  a,  r  and 
wn  are  determined  by  substituting  (3.14)  with  (3.15)  into  (3.13), 
expanding  in  powers  of  x,  and  equating  coefficients  of  like  powers. 

The  significance  of  (3.14)  is  that  it  sets  the  cross-stream  decay 
scale  through  the  parameter  a.  This  scale  length  turns  out  upon  ex¬ 
pansion  to  be  just  1/k.  Determining  it  explicitly  is  also  crucial  to 
the  normalization  scheme  Introduced  in  Chapter  6.  In  addition,  the 
parameter  a  serves  to  connect  the  asymptotic  form  of  the  solution  for 
x  •  to  an  integral  representation  of  the  solution  for  finite  x. 
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For  the  case  of  no  cross-stream  flow,  (3.14)  becomes 

Z(z)  -v  z_Ctd "e”2^2*  (l-(ctj/ z)+. . .  } ,  (3.16) 

a 

where 

ad  *  %tl+(f2-u>2)  /  a^g  'k-f /go}  (3.17) 

and  z  -  (2k/a1)(Df(0)+a1x}  (3.18) 

are  the  same  as  Duxbury's  (1963)  quantities  a  and  z.  (3.16),  which  is 
written  in  Duxbury's  notation,  recovers  his  expression  for  the  asymp¬ 
totic  behavior  of  the  DPA  which  appears  in  Appendix  A  of  his  paper. 

Note  that  z  above  should  not  be  confused  with  the  vertical  coordinate, 
nor  with  a  in  (3.14).  This  notation  is  employed  only  in  (3.17), 
(3.18)  and  (3.21)  in  order  to  compare  these  equations  with  Duxbury's 
expressions.  From  (3.16)  and  (3.18),  we  see  clearly  that  the  natural 
scale  in  the  cross-stream  direction  is  the  decay  scale  for  the  DPA,  1/k. 

Since  only  the  point  at  infinity  is  an  irregular  singularity  of 
the  DPA  ODE,  (3.13),  a  valid  power  series  expansion  of  Z(x)  exists  in 
any  finite  and  bounded  interval  in  x.  The  series  may  be  developed  by 
employing  the  method  of  Frobenius,  a  technique  which  is  straightforward 
in  principle  and  well  suited  to  numerical  computation.  This  procedure 
is  not  particularly  illuminating  from  an  analytical  viewpoint,  however, 
since  it  does  not  provide  the  relationship  between  the  series  expansion 
for  finite  x  and  the  asymptotic  form  of  Z(x)  as  x  -*■  ®  in  (3.14).  This 
limitation  is  overcome  by  Laplace  transforming  (3.13)  and  using  the 
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the  property  that  its  coefficients  are  at  most  linear  in  x,  whereupon 
the  solution  for  Z(x)  may  be  written 

r  «_ 

1  •  exp{t(x+Y7/v7) }  dt  .  (3.19) 

distinct,  are  the  roots  of  the 

-  0  ,  (3.20) 

and  the  integration  is  carried  out  along  a  line  parallel  to  the  real 
axis  in  the  complex  t-space.  The  A^  and  Bq  in  (3.19)  are  constants. 

For  the  case  of  zero  cross-stream  flow,  (3.19)  reduces  to  an  integral 
representation  of  the  confluent  hypergeometric  function  U  (Abramowitz 
and  Stegun,  1970,  Ch.  13).  Z(x)  then  becomes 

-z/2 

Z(x)  -  Be  Z/  -U(ad,l;z)  ,  (3.21) 

which  is  precisely  Duxbury's  equation  (66)  in  which  the  solution  that 
is  singular  at  infinity  has  been  removed.  The  a's  in  (3.19)  are  exactly 
the  same  a's  appearing  in  the  asymptotic  form  (3.14),  thereby  estab¬ 
lishing  the  connection  between  the  solution  for  finite  x  and  its  be¬ 
havior  at  Infinity. 

In  principle,  the  problem  of  the  planar  front  with  cross-stream 
flow  has  been  solved,  since  the  determination  of  the  DPA  Z(x)  has  been 
reduced  to  quadratures,  and  the  associated  perturbations  to  the  flow 
field  may  then  be  computed  from  (3.6).  Note  that  once  Z(x),  and 
consequently  its  first  derivative,  are  known  explicitly,  (3.6)  reduces 
to  an  algebraic  system  for  the  three  unknown  functions  X(x),  X'(x) 
and  Y (x) . 


CHAPTER  4 


NON-DIMENS IONALIZATION  OF  THE  PLANAR  FRONT  ODE  FOR  Z(x) 

4.1  INTRODUCTION 

In  Chapter  3,  a  single  ODE  was  developed  for  the  depth  perturba¬ 
tion  amplitude  in  the  planar  front  case.  The  analysis  was  presented 
entirely  in  dimensional  variables,  as  was  Duxbury's  (1963)  tre;  ment 
of  the  special  case  of  zero  cross-stream  flow.  It  led  to  an  asymptotic 
normal  form  solution  for  large  x,  an  integral  representation  of  the 
solution  for  finite  x,  and  the  connection  between  the  two  via  the 
parameters  a^.  In  principle,  then,  the  formal  solution  of  the  planar 
front  problem  with  allowance  for  cross-stream  flow  has  been  obtained, 
at  least  in  the  sense  that  it  has  been  reduced  to  quadratures.  In 
Duxbury's  model,  the  cross-stream  flow  was  taken  to  be  zero  at  the 
outset,  and  the  resulting  ODE  for  the  depth  perturbation  amplitude  was 
second  order.  In  this  model,  Duxbury's  second  order  ODE  and  all  of 
his  subsequent  results  are  recovered  if  the  cross-stream  velocity  com¬ 
ponent  vanishes;  however,  the  ODE  for  the  depth  perturbation  amplitude 
is  inherently  order  seven,  not  two,  as  long  as  any  cross-stream  flow 
exists,  however  small.  This  observation  is  somewhat  disquieting,  and 
is  reminiscent  of  the  singular  perturbation  problems  that  frequently 
occur  in  fluid  dynamics  in  which  the  order  of  a  system  changes  abruptly 
as  some  parameter  of  the  problem  becomes  zero.  This  peculiarity 
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demands  closer  examination,  and  motivates  the  non-dimens ionalization 

of  the  DPA  equation  that  is  presented  in  this  chapter. 

It  was  pointed  out  earlier  that  the  frontal  zone  is  characterized 

by  the  baroclinic  Rossby  radius  length  scale,  X,  on  which  the  dynamics 

are  invlscid  except  in  a  narrow  dissipative  region  near  the  surface 

front.  A  second  natural  length  scale  also  exists  for  the  problem, 

however,  the  ratio  of  the  basic  state  cross-stream  velocity  to  the 

Coriolis  parameter.  This  is  an  advectlve  scale  and  is  much  smaller 

than  the  Rossby  radius  scale.  It  is  the  scale  on  which  the  cross-stream 

Rossby  number  becomes  order  one;  i.e. ,  u/f.  In  a  typical  example  with 
-  -4 

u  ■  10  cm/s  and  f  *  10  /s,  this  scale  would  be  only  1  km.  Thus,  it 

would  generally  lie  well  within  the  inner  dissipative  zone  so  that  it 
would  not  require  separate  treatment  in  this  paper.  In  what  follows, 
the  ODE  for  the  planar  front  depth  perturbation  amplitude  will  be  non- 
dimensionalized  by  Introducing  coordinates  scaled  to  the  outer 
(inviscid)  and  advectlve  cross-stream  scale  lengths. 

4.2  INVISCID  SCALE  VARIABLES 

For  time  varying  wave  disturbances  on  the  planar  front,  the 

natural  long-stream  length  scale  is  clearly  the  disturbance  wavelength, 

or,  equivalently,  the  reciprocal  wavenumber  1/k.  From  the  asymptotic 

form  of  the  solution  to  the  DPA  in  Section  3.5,  in  which  an  explicit 
-lex 

e  dependence  appears,  it  is  clear  that  the  decay  scale  for  the  DPA 
in  the  cross-stream  direction  is  also  1/k,  and,  consequently,  that  the 
x-  and  y-  scales  for  wavelike  disturbances  of  the  planar  front  are  the 
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same.  A  careful  scaling  analysis  (cf.  Section  5.2  for  details),  more¬ 
over,  reveals  that  the  only  nontrivial  perturbation  of  interest  cor¬ 
responds  to  k  ^  1/A,  where  A  is  the  baroclinic  Rossby  radius  defined 
below.  Dimensionless  x  and  y  coordinates,  marked  by  a  tilde  (,v) ,  are 
therefore  introduced  as  follows: 

x  =  x/A  (4 . la) 

y  =  ky  (4. lb) 

The  natural  time  scale  is  the  inertial  period  1/f,  and  the  appropriate 
dimensionless  circular  frequency  is  thus 

o  =  oj/f  (4.  lc) 

The  advective  scale  is  characterized  by  the  x-  and  y-  direction 
lengths  u/ f  and  v/f,  respectively,  where  u  and  v  are  the  equilibrium 
state  velocities  introduced  in  Chapter  2.  The  ratios  of  the  advective 
to  inviscid  length  scales  for  the  x  and  y  directions  form  two  non- 
dimensional  parameters  of  the  problem,  which  are  denoted 

e  *  u/f  (4. Id) 

and 

e’*  kv/f  (4 . le) 

In  general,  we  anticipate  e  <<  1  while  e’^  1,  since  |u/v|  <<  1  for 
large  scale  fronts  with  large  values  of  the  rotation  parameter. 

Two  internal  Froude  numbers  exist  which  measure  the  ratio  of  the 
fluid  velocity  to  the  internal  wave  phase  speed.  The  latter  is  defined 

c2  =  g'Df(0)  (4. If) 

The  two  Froude  numbers  are  then, 

Fx  *  u/c  (4 . lg) 

for  the  x-direction,  and 

Fy  -  v/c  (4. lh) 

for  the  y-direction.  We  anticipate  Fx  <<  1  while  Fy  'v  1  for  the  same 
reasons  as  above  for  c  and  e'. 


(4.1i) 
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The  baroclinic  Rossby  radius  Is  defined  as 

A  -  c/f. 

The  free  parameter  6  appearing  In  the  basic  state  velocity  field 
developed  in  Chapter  3  for  the  planar  front  is  absorbed  into  a  new 
dimensionless  transport  parameter  defined  by 

*  *  1-8“ \  (4.1j) 

which  is  order  one. 

With  these  definitions,  the  non-dimensional  differential  equation 

for  the  planar  front  DPA  in  inviscid  scale  variables  becomes 

7 

Z  f\j  ^  %  n  'V,  'W* 

(v  x+Y  ) (d  Z(x)/dx)  -  0,  (4.2) 

n  n 

n**0 

% 

where  the  dimensionless  complex  constants  v  and  y  are  defined  in 

n  n 

Appendix  2B. 

We  now  assume  that  Z(x)  in  (4.2)  may  be  developed  in  a  small 
parameter  expansion  as  follows 


Z(x)  -  Z(0)+F  Z(1)+F2Z(2)+...  , 
x  x 


(4.3) 


in  view  of  which,  at  order  F^,  (4.2)  becomes 


(x+x { l-o2- (kA) 2 } /kAF-l/<J } Z (0) (x)  -  dZ(0)/dx 


+{kAT/Fy-x}(d2Z(0)/dx2)  -  0. 


(4.4) 


The  case  of  zero  cross-stream  flow  corresponds  to  F^  «  0  in  (4.2) 
and  (4.3),  which  again  recovers  Duxbury’s  (1963)  second  order  ODE  for 
Z(x),  as  it  must.  The  details  are  not  reproduced  here.  Even  for 
nonzero  F  ,  however,  we  see  that  at  order  F°  the  equation  for  Z(x)  is 
still  second  order,  not  seventh  order  as  in  the  general,  dimensional 
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formulation.  On  the  outer  (inviscid)  scale,  therefore,  which  is  the 
only  region  of  applicability  of  this  model,  the  lowest  order  ODE  for 
the  depth  perturbation  amplitude  is  still  second  order.  An  extension 

*\j 

to  first  order  in  F^  leads  to  a  similar  result:  the  equation  for  Z(x) 
is  still  second  order,  but  with  the  addition  of  lowest  order  forcing 
terms.  The  first  order  dynamics  in  the  small  parameter  F^  are  thus  no 
longer  homogeneous,  but  instead  are  driven  by  the  zeroth  order  solu¬ 
tions.  On  the  outer  scale,  the  third  and  higher  derivative  terms  in  Z 

2 

can'  enter  the  dynamics  at  most  at  order  F^  in  the  small  parameter  ex¬ 
pansion,  if  at  all.  Thus,  the  inviscid  scale  dynamics  of  the  planar 
front  are  fundamentally  second  order,  even  with  nonzero  cross-stream 
flow,  a  conclusion  which  is  consistent  with  Duxbury's  second  order 
model.  This  conclusion,  however,  has  been  reached  only  by  having 
introduced  appropriate  scalings  for  the  problem,  not  by  developing  the 
general  dimensional  form  presented  in  Chapter  3  or  in  Duxbury's  (1963) 
work. 

The  effect  of  nonzero  cross-stream  flow  on  the  outer  scale  is 
simply  to  modify  the  constants  in  the  coefficients  of  the  ODE  for  Z(x) 

by  introducing  terms  which  vanish  with  zero  cross-stream  flow;  the 

v 

functional  dependence  of  the  coefficients  on  x  remains  unchanged  even 
for  nonzero  cross-stream  flow.  The  formalism  developed  in  Chapter  3 
for  the  solution  of  the  depth  perturbation  amplitude  ODE  is  therefore 
directly  applicable  to  the  DPA  equation  at  lowest  order  in  the  para¬ 
meter  F  ,  but  the  results  thus  obtained  are  only  a  minor  extension  of 
Duxbury's  original  problem. 
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It  is  interesting  to  note  that  normalizing  the  cross-stream 
coordinate  to  the  advectlve  length  scale,  u/f,  and  redoing  the  above 

*\j  0 

analysis  leads  to  a  fifth  order  ODE  for  Z(x)  at  order  F^.  The  details 
are  not  reproduced  here,  but  this  observation  is  important,  since  it 
shows  that  the  higher  order  derivative  terms  in  the  DPA  equation  are 
significant  at  order  F^  only  in  the  very  narrow  advectlve  zone.  How¬ 
ever,  since  this  region  is  imbedded  in  the  dissipative  zone,  the  kind 
of  inviscld  dynamics  which  would  develop  on  this  scale  are  never  real¬ 
ized. 

The  Insight  provided  by  this  analysis  will  be  incorporated  in  what 
follows  in  order  to  develop  a  more  realistic  frontal  model  that  in¬ 
cludes  nonzero  cross-stream  flow  and  more  accurately  reflects  observed 
frontal  zone  geometry. 


CHAPTER  5 


THE  BASIC  STATE-  FOR  THE  EXPONENTIAL  FRONT 

5.1  INTRODUCTION 

The  planar  front  model  studied  in  Chapters  2  through  4  has  several 
physical  shortcomings.  The  interfacial  depth  cannot  increase  indefi¬ 
nitely  away  from  the  surface  front  in  any  realistic  frontal  model,  as 
it  does  in  the  planar  front,  unless  it  intercepts  the  bottom.  In 
nature,  the  interfacial  depth  gradually  approaches  some  fairly  constant 
value  far  from  the  surface  front  in  the  cross-stream  direction.  A 
model  possessing  this  feature  was  first  proposed  by  Stommel  (1976,  Ch. 

8)  for  the  Gulf  Stream  and  will  be  referred  to  here  as  the  "exponential 
front".  The  frontal  zone  geometry  in  this  case  is  shown  in  Figure  5.1, 
which  is  a  cross-sectional  view  in  the  plane  y  =  0.  All  of  the  restric¬ 
tions  introduced  earlier  are  also  assumed  here,  and  all  subsequent 
equations  are  written  in  the  moving  coordinate  system  first  introduced 
in  Chapter  2. 

Referring  to  Figure  5.1,  the  interfacial  depth,  which  is  denoted 
D^(0)  at  the  origin,  as  before,  gradually  increases  in  the  cross-stream 
(x)  direction  and  asymptotically  approaches  the  value  far  from  the 
dissipative  zone.  All  dissipation  through  mass  entrainment,  mixing, 
and  turbulent  friction  is  again  confined  to  the  narrow  dissipative 
region.  The  interfacial  depth  increases  from  D^(0)  to  in  the 
inviscid  (outer)  region  on  its  natural  length  scale.  A,  the  baroclinic 
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CROSS-SECTIONAL  VIEW  OF  THE 
EQUILIBRIUM  STATE  EXPONENTIAL 
FRONTAL  ZONE  GEOMETRY 


FIGURE  5.1 
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Rossby  radius.  With  the  spatial  coordinates  scaled  by  X  we  anticipate 
that  spatial  derivatives  are  order  one  in  the  x-direction  and  very 
small  in  the  y-direction. 


5.2  DIFFERENTIAL  EQUATIONS  FOR  THE  EXPONENTIAL  BASIC  STATE 

The  general  time  dependent  momentum  and  continuity  equations  for 
the  inviscid  (outer)  region  appear  in  (2.17)  and  (2.18)  and  are  repro¬ 
duced  here  for  convenient  reference: 

u  +uu  +vu  -fv  =  -g'D  -fv  (5.1a) 

txy  x  00 

v  +uv  +w  +fu  =  -g'D  +fu  (5.1b) 

txy  y  00 

D  +(uD)  +(vD)  -  0  (5.1c) 

c  x  y 


in  which  all  variables  are  dimensional.  (5.1)  makes  no  assumption  about 
the  frontal  zone  geometry. 

Differentiating  the  x-momentum  equation  with  respect  to  y  and  the 
y-momentum  equation  with  respect  to  x,  subtracting  and  using  (5.1c) 
leads  to  the  result 


£_  (_f±A_) 

Dt  '  D  ' 


(5.2) 


where  A  =  v  -  u  and  where  D/Dt  is  the  convective  derivative.  The 
x  y 

quantity  (f  +  A)/D  is  recognized  as  the  potential  vorticity,  and  (5.2) 
is  a  statement  of  conservation  of  potential  vorticity  following  a  fluid 
element.  This  result  follows  from  the  assumption  of  inviscid  dynamics 
which  permit  no  mechanism  for  altering  the  angular  momentum  of  a  fluid 
column. 

The  equilibrium  state  equations  for  the  exponential  front  are 
formulated  by  setting  —  -  0  and  uy  ■  vy  -  (vD)y  ■  0  in  (5.1),  reflec¬ 
ting  the  fact  that,  in  the  basic  state,  y-variations  are  assumed  very 
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small  compared  to  x- variations.  Although  the  long  stream  rate  of 
change  of  the  y-direction  volume  flux,  vD,  is  negligible  compared  to 
the  cross-stream  rate  of  change  of  the  x-direction  volume  flux,  uD, 
the  derivative  in  the  y-momentum  equation  is  retained.  This  is 
necessary  since  must  be  present  for  large  x  to  provide  the  y-pressure 
gradient  which  is  compatible  with  the  cross-stream  flow  u^  there.  The 
The  simplified  basic  state  continuity  and  momentum  equations  thus  be- 


(uD)  -  0 
x 

uu  -fv  «  -g'D  -fv 

X  X  00 

uv  +fu  ■  -g'D  +fu  , 
x  6  y  »  * 

while  the  potential  vorticity  equation  is  now 


(5.3a) 

(5.3b) 

(5.3c) 


(3/3x) { (f+v  )/D)  -  0. 
x 


(5.3d) 


Each  of  the  dependent  variables,  u,  v,  and  D  in  (5.3)  is  a  function 
of  x  and  y,  but  the  y-dependence  is  weak.  In  addition,  since  (5.3d) 
was  derived  by  combining  (5.3b)  and  (5.3c),  only  two  of  these  three 
equations  are  independent.  The  two  independent  equations  are  here 
taken  to  be  the  x-momentum  and  potential  vorticity  equations. 

(5.3a)  and  (5.3d)  may  be  integrated  immediately  to  obtain 


uD  *  u,  D, 


(5.4) 


f+v  *  f 5  , 

x 


(5.5) 


where  u^  is  the  (constant)  x-directlon  velocity  in  the  frontal  zone  far 
from  the  dissipative  region,  and  6  is  the  non-dimensional  interfacial 


depth  defined  by 


5  -  D/Dfe  . 


(5.6) 


Differentiating  the  x-momentum  equation  (5.3b)  with  respect  to  x 
and  substituting  for  vx  and  u  from  (5.4)  and  (5.5)  produces  the 
following  equation  for  the  normalized  basic  state  interfacial  depth: 


(1-F^/63)6cc+(3F^/64)6^+1-6  =  0  , 

(5.7) 

where  l* 

c  -  (g’Dbr 

(5.8a) 

Fb  ■  “b/c 

(5.8b) 

X  =  c/f 

(5.8c) 

5  =  x/X 

(5. 8d) 

In  (5.8),  c  is  the  frontal  zone  internal  wave 

phase  speed  (note 

that  it  is  defined  somewhat  differently  here  than  in  Chapter  4),  is 
the  x-direction  Froude  number  and  measures  the  strength  of  the  frontal 
zone  cross-stream  flow,  X  is  the  baroclinic  Rossby  radius  of  deforma¬ 
tion,  and  5  is  the  normalized  x-coordinate.  Under  the  linearization 
scheme  that  is  introduced  in  Chapter  6,  F^  constitutes  the  fundamental 
parameter  for  the  frontal  zone. 

In  general,  we  seek  solutions  for  which  the  disturbance  frequency 

is  small  compared  to  the  inertial  frequency;  i.e.,  we  are  interested  in 

low  subinertial  frequency  behavior.  Similarly,  we  are  interested  in 

fronts  for  which  the  inviscid  region  scale  length  is  of  order  X.  k 

careful  scaling  analysis  applied  to  the  general  equations  in  (5.1) 

2 

then  reveals  that  the  only  conditions  of  interest  are  that  kX  'v  F  ^  1, 
where  F  ,  is  the  y-direction  internal  Froude  number  (ratio  of  the 

y 

basic  state  long  stream  velocity  to  the  internal  wave  phase  speed) .  The 
reciprocal  wavenumber,  1/k,  is  therefore  of  the  order  of  the  Rossby 
radius,  X,  and  either  of  these  length  scales  may  be  used  to  normalize 
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the  inviscid  scale  spatial  coordinates.  (5.8d)  is  thus  consistent 
with  the  earlier  observation  for  the  planar  front  that  the  inviscid 
scale  length  is  order  1/k.  For  large  scale  upper  ocean  fronts,  ob¬ 
servations  tend  to  show  |F^|  <<  1,  while  for  smaller  scale  fronts 
|Fj  may  be  order  1  (Garvine,  1979b).  In  addition,  consistent  ob¬ 
servation  of  the  high  velocity  long-stream  shear  characteristics  of 
large  scale  fronts  is  reflected  in  the  condition  |F^|  v  l.  For 
analytical  convenience,  however,  any  explicit  F^  dependence  will  be 
avoided  in  the  sequel  by  choosing  F^  =  1  identically.  This  is  tanta¬ 
mount  to  scaling  the  long-stream  velocity  by  c. 

(5.7)  provides  now  a  second  order,  nonlinear,  inhomogeneous  ODE 
for  the  non-dimensional  interfacial  depth,  6.  Recall  that  5  is  only 
weakly  y-dependent,  and  that  its  y-dependence  may  be  computed  from  the 
y-momentum  equation  once  the  remaining  three  equations  in  the  system 
(5.3)  have  been  solved.  The  solution  to  (5.7)  will  be  developed  in 
the  next  section. 

5.3  SOLUTION  OF  THE  BASIC  STATE  INTERFACIAL  DEPTH  EQUATION 
Subject  to  the  boundary  conditions 

Lim  5(0  =  1  (5.9a) 

O* 

5(0)  -  5o  -  Df(0)/Db  , 


- 


. . . . — — ■ 


(5.9b) 


the  exact  solution  for  6(5)  in  (5.7)  is  given  implicitly  by 


5  -  ('9/ 6q)  -  (4>/6)  +  ln{ (6  +f) /(*+«)} 

+F,  (sin'V. /«)  -  ;vin_1  (F  /6  )} 

D  D  DO 

-  (l-F^) 3/2  In  |  («-1){6o-fJ+(1-F^)JjT}/ 


(5.10a) 


where 

*  ■  <5. 10b) 

*  -  (5.10c) 


Principal  values  are  to  be  understood  foi  he  arcsine  terms,  and  we 
2 

require  F^  <<  6  <1.  Note  also  that  6  ?  order  one. 

bo  o 

Although  the  exact  solution  is  very  u*ful  for  numerical  compu¬ 
tation,  it  is  very  cumbersome  for  analytic*  purposes,  since  the  de¬ 
pendent  variable  is  represented  implicitly.  ,\n  approximate  ,  but 
explicit  representation  may  be  developed,  he  ever,  by  introducting  a 
small  parameter  expansion  for  6(5)  in  terms  the  Froude  number  F^. 
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We  assume  that  6(0  may  be  written 

6(C)  -  6(0)  (O  +  F26(2)(£)  +  F46(4)(f-)  +  ...  ,  (5.11) 

2 

in  which  only  even  powers  of  Ffe  appear,  since  (5.7)  contains  only  Fj\ 
Each  of  the  functions  6^(0  is  order  one,  and  the  error  introduced 
by  truncating  the  series  at  6 is  order  f£.  Substituting  (5.11) 
into  (5.7)  yields  the  following  equations  for  6^  and  6^: 


«<°>  -  ««»  +  1  -  0 

««>  -  s(2>  -  -wW)!  1/(S(0))2 


The  approximate  solution  will  not  be  carried  beyond  second  order  in 
F^.  (5.12)  and  (5.13)  are  both  second  order,  linear  and  inhomogeneous. 

The  inhomogeneity  in  (5.12)  is  a  result  of  the  inhomogeneity  in  the 
original  equation  for  6(0 ,  while  the  Inhomogeneity  in  (5.13)  is 
characteristic  of  the  small  parameter  expansion  formulation,  in  which 
each  successive  solution  beyond  zeroth  order  is  dependent  on  the  lower 
order  solutions. 

The  boundary  conditions  (5.9)  will  be  satisfied  exactly  to  zeroth 
order  in  F^,  that  is. 


6(0) 


(0)  -  6q 


(5.14a) 


Lim  6  (0)  (O  «  1  . 


(5.14b) 


61 


In  view  of  (5.9)  and  (5.14),  all  higher  order  terms  in  F,  in  the 

b 

expansion  for  6(C)  must  satisfy  zero  boundary  conditions,  that  is, 


6(n)(0)  =  0 


Lim  6(n)(c)  =  0 


n  >  2 


(5.15a) 

(5.15b) 


The  lowest  order  solution  to  (5.12)  subject  to  the  boundary 
conditions  (5.14)  is 

6 (0) (C)  *  1  -  (1-6  )e‘C  ,  (5.16) 

o 


which  is,  of  course,  also  the  solution  to  the  original  equation  for 
6(C)  with  F^  =  0.  This  is  the  form  of  the  solution  found  by  Stommel 
in  his  Gulf  Stream  model  and  it  motivates  the  "exponential  front" 
nomenclature.  The  second  order  solution  to  (5.13)  subject  to  the 
boundary  conditions  (5.15)  is 


6(2)(C)/{6(0)-l}  =  In 


(1-5 (0)) 


1  (l-6o)6 


(0) 


(6(-0)-6  )(6(0)+6  +2) 


(5.17) 


The  set  of  equations  (5.11),  (5.16)  and  (5.17)  provides  an  expli¬ 
cit,  although  approximate,  representation  for  the  time  invariant  basic 
state  interfacial  depth  to  second  order  in  the  small  parameter  F,  .  The 

D 

difference  between  the  exact  and  approximate  solutions  is  expected  to 
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2 

be  order  F.  .  Having  developed  an  °.xplicit  form  for  <5  (5),  the  basic 

D 

state  x  and  y  velocities  in  the  fronts’ 1  zone  may  now  be  computed  from 
(5.4)  and  (5.3b)  if  desired. 

Figures  5.2a  through  5. 2e  compare  the  exact  equilibrium  inter¬ 
facial  depth  computed  from  (5.10)  (curve  A)  with  the  approximate  repre¬ 
sentation  computed  to  order  F^  using  (5.11)  with  (5.16)  and  (5.17) 
(curve  B).  The  value  of  5  is  fixed  at  0.4,  and  F.  ranges  over  0.1 

O  D 

to  0.3  in  steps  of  0.05.  For  F,  »  0.1  and  0.15,  the  approximate 

D 

solution  is  greater  than  the  exact  solution,  while  at  F^  *  0.25  and 
0.3,  this  is  reversed.  For  F^  *  0.2,  the  two  curves  are  essentially 
indistinguishable.  The  absolute  error  in  the  approximate  solution 

*\j 

therefore  passes  through  a  minimum  near  F^  =  0.2.  Agreement  between 

2 

the  two  solutions  at  order  F^  is  excellent  over  the  range  of  F^  that 

2 

is  considered,  and  the  error  is  generally  even  less  than  F^.  Although 
several  other  cases  were  computed,  the  results  are  not  included  here, 
since  Figure  5.2  serves  well  to  illustrate  the  behavior  of  the  approxi¬ 
mate  and  exact  solutions. 
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I 


< 


Figure  5.2(b) 


CHAPTER  6 


LINEARIZATION  OF  THE  TIME-VARYING  EQUATIONS  FOR  THE 
EXPONENTIAL  AND  PLANAR  FRONTS  IN  DIMENSIONLESS  VARIABLES 

6.1  INTRODUCTION 

In  this  chapter,  the  dependent  variables  in  the  time  varying 
momentum  and  continuity  equations  are  decomposed  into  basic  state  and 
perturbation  components,  as  was  done  in  Chapter  2,  but  with  the  intro¬ 
duction  of  dimensionless  variables.  The  resulting  perturbation  equa¬ 
tions  are  then  rewritten  using  a  small  parameter  expansion.  Both  the 
exponential  front,  whose  basic  state  was  formulated  in  the  previous 
chapter,  and  the  planar  front,  which  was  discussed  at  length  in  Chapters 
3  and  4,  are  considered.  The  perturbation  variable  decomposition  is 
the  same  as  that  introduced  in  Chapter  2,  and  the  restriction  that  the 
perturbation  amplitude  be  small  compared  to  the  basic  state  value  of  a 
given  variable  is  again  implicit  in  all  of  what  follows.  Since  the 
equilibrium  states  for  the  exponential  and  planar  fronts  are  different, 
the  equations  for  the  perturbation  variables  in  these  two  cases  are 
also  different,  and  therefore  examined  separately.  However,  essential¬ 
ly  the  same  non-dimensionalizatiom  scheme  is  employed  for  each  case, 
since  this  depends  only  on  the  frontal  zone  velocity,  length  and  time 
scales,  which,  in  turn,  are  only  weakly  dependent  upon  the  model  geo¬ 
metry.  The  perturbation  variable  decomposition  requires  an  explicit 
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form  for  the  basic  state  variables  u,  v  and  D,  which  then  appear  as 
coefficients  in  the  perturbation  equations.  This  is  accomplished  by 
introducing  the  same  small  parameter  expansion  in  terms  of  F^  that  was 
used  in  the  previous  chapter  for  both  the  exponential  and  planar  front 
cases.  This  provides  the  additional  advantage  of  facilitating  compari¬ 
son  between  the  two. 

! 

< 

6.2  PERTURBATION  VARIABLES  AND  NON-DIMENSIONALIZATION 

The  time  varying  momentum  and  continuity  equations  are  (2.17a), 

(2.17b),  and  (2.18)  for  both  the  exponential  and  planar  fronts,  in 
which  u,  v  and  D  are  each  functions  of  (x,  y,  t) .  These  equations  are 
general,  and  apply  to  the  basic  state  as  well  as  the  time  varying  flow. 

As  before,  each  of  the  dependent  variables  is  decomposed  into  a  basic 
state  component  and  a  perturbation  component  which  is  small  by  compari¬ 
son,  that  is,  the  decomposition  (2.19)  is  again  assumed.  For  the  ex¬ 
ponential  front,  the  resulting  perturbation  equations  are  different 
from  those  for  the  planar  front,  since  the  spatial  variation  of  the 
equilibrium  velocity  field  introduces  derivatives  of  the  basic  state 
flow. 

The  variables  for  the  planar  and  exponential  front  models  are 
scaled  as  follows:  x  is  scaled  by  X,  y  by  1/k,  the  interfacial  depth 
by  D^(0)  or  D^,  u  by  the  characteristic  x  velocity  u^,  v  by  c,  and  t 
by  the  inertial  time  scale  1/f.  As  before,  k  is  the  (dimensional) 
y-direction  wavenumber  for  harmonic  frontal  disturbances,  and  X  is  the 
(dimensional)  baroclinic  Rossby  radius.  The  product  kX  is,  therefore, 
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dimensionless.  The  long-stream  equilibrium  state  velocity  scale 
reflects  the  characteristic  frontal  zone  jet  predicted  by  other 
authors  (Charney,  1955;  Garvine,  1979a}  Morgan,  1956).  Since  per¬ 
turbations  to  the  system  are  homogeneous,  as  we  have  seen  in  section 
2.3,  the  normalization  for  the  perturbation  velocities  is  arbitrary, 
and  therefore  set  equal  to  c  for  convenience.  The  dimensionless  x 
and  y  coordinates  are  £  and  t,  respectively,  while  6  is  the  dimension¬ 
less  interfacial  depth.  F^  is  the  x-direction  internal  Froude  number 
introduced  in  Chapter  5,  while  its  y-direction  counterpart,  F^  in  the 
previous  chapter,  is  set  identically  equal  to  1.  Other  non-dimensional 
variables  are  marked  by  an  asterisk.  Thus,  we  define 


x  -  X?  ,  y  »  c/k  , 


_* 

u  =  u,  u 


-*  A  * 
v  -  CV  ,  U  ■  CU  ,  V  »  cv 


u«  *  v«  ’  v»  ■  cv»  *  Fb  "  ub/c 


t  *  t  /f  ,  X  *  c/f  (Rossby  radius) 


5  -  D.  6  I 

b  I 

A  A  I 

D  -  D.  6  | 

c’-k  I 

5  -  Df (0)6 
D  -  Df (0)6 
c2  -  g'Df(0) 


Exponential  Front 


Planar  Front 


(6.1) 
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In  view  of  (6.1),  the  fundamental  system  of  perturbation  equations 
for  both  the  planar  and  exponential  fronts  in  non-dimensional  variables 
is: 

fl  +  F,  uu_  +  kA  vur  +  F,  u  u  -  v  +  5  =  0  (6.2a) 

t  b  4  4  b  4  4 

v  +  F  uv  +  kAvv  +  v  u  +  u  +  kA6  =  0  (6.2b) 

t  b  4  4  4  4 

l  +  F  (uS)  _  +  kAv6  +  (6u)  _  +  kA(6<»  =  0  (6.2c) 

C  D  4  4  4  4 


The  planar  front  equilibrium  (overbar)  velocities  are  then 

u  =  (1-6)0^  (6.3a) 

v  =  (l-6)vra  (6.3b) 


where  6  is  the  dimensionless  transport  parameter  introduced  in  (3.1). 
Note  that  the  asterisks  have  now  been  dropped  for  notational  conveni¬ 
ence,  and  it  is  understood  that  all  quantities,  except  k  and  A,  which 
always  appear  as  a  product,  are  non-dimensional.  The  planar  front 
basic  state  interfacial  depth  is 


6(5,4)  =  6(0)  +  Fb5(1)  ,  (6.4a) 

where 


<5(0)  -  l-6v  4 

00 

(6.4b) 

6(1)  =  (BF^/kA)?  . 

(6.4c) 

Note  that  the  cross-stream  flow  enters  the  scaled  basic  state  only  at 
order  F^. 

D 
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For  the  exponential  front,  the  scaled  basic  state  (overbar)  vari¬ 
ables,  which  were  developed  explicitly  in  the  previous  chapter,  are 


given  by 


u  *  1/5 


(6.5a) 


v  -  6  -  F*5c/53  +  v 

£  b  £  <* 


(6.5b) 


6(0  -  6(0)  +  F26(2) 

D 


(6.5c) 


where  6^°^  and  5^  appear  in  (5.16)  and  (5.17),  respectively.  The 
basic  state  flow  is  derived  from  (5.3b)  and  (5.4)  using  (5.11). 


6.3  SMALL-PARAMETER  EXPANSIONS 

Each  of  the  dependent  variables  in  (6.2)  is  now  expanded  in 
terms  of  the  small  parameter  F^  as  follows: 


u  -  u(0)  +  ?bu(1)  +  F2u(2)  +  ... 

0  -  v(0)  +  F  0(1)  +  F2v(2)  +  . 

6  -  6<°>  +  F/1*  +  F ^(2)  +  ... 

5-6<°>+F26<2>+... 

3  -  6^  +  F  5^^ 

b 


(6.6a) 
(6.6b) 
(6.6c) 
(6 . 6d) 
(6 . 6e) 


where  (6.6d)  and  (6.6e)  apply  to  the  exponential  and  planar  fronts, 
respectively. 

In  (6.6),  each  of  the  functions  {i^°\  0^2\  etc.,  is  presumably  order 
one,  as  is  the  product  kX.  The  exponential  front  expansion  in  (6.6d) 
is  the  same  as  (5.11)  in  which  the  overbar  had  not  yet  been  used,  and 
only  even  powers  of  F^  appear  for  the  reason  discussed  in  section  5.3. 
For  the  planar  front,  the  basic  state  expansion  for  the  interfacial 
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depth  terminates  at  order  F^,  even  though  the  expression  is  exact;  and 
none  is  introduced  for  the  velocity  field,  since  it  is  constant. 

Substituting  (6.6)  into  (6.2)  leads  to  the  following  zeroth  and 
first  order  non-dimensional  perturbation  equations  for  the  planar  front: 


u<°>  +  kXiif)  -  v<0)  +  ;<°>  -  0 


«"»  +  +  a<0)  +  kxi(0)  -  0 

t  ■  '  c  £ 


i(0)  +  kx;i<°>  -  Bv  u(0)  +  (1-Bv  «a<°> 
t  C  00  °°  £ 

+  kA (l-8v  c)vf0)  =  0 
00  £ 

(PLANAR  FRONT,  ORDER  F°) 

b 


(6.7a) 

(6.7b) 

(6.7c) 


To  order  F°  for  the  planar  front,  the  basic  state  enters  as  nonconstant 

D 

coefficients  in  the  resulting  perturbation  equations.  The  lowest  order 
system,  however,  is  homogeneous,  whereas  the  first  order  system  is  not, 
as  a  result  of  forcing  by  the  lowest  order  terms.  The  lowest  order 
solutions  also  appear  as  nonconstant  coefficients  in  the  first  order 


system 


The  corresponding  systems  of  zeroth  and  first  order  equations  for 


the  exponential  front,  which  are  derived  by  using  (6.6)  in  (6.2)  and 


then  eliminating  the  basic  state  velocities  with  (6.5),  are  as  follows 


(EXPONENTIAL  FRONT,  ORDER  FH 


+  kX(l+v  -<5V  ')u 


+  kX(l+v  -<5 


+  kX(l+v  -<5V  ')<5 


(EXPONENTIAL  FRONT,  ORDER  F*) 


■ 
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As  with  the  planar  front,  the  lowest  order  system  is  homogeneous, 
while  the  first  order  system  is  not.  Note  that  the  coefficients  in 
(6.10)  have  been  simplified  somewhat  by  using  the  basic  state  relations 


and 


(6.11a) 


(6.11b) 


6.4  ZEROTH  ORDER  HARMONIC  SOLUTIONS 

In  order  to  simplify  the  notation,  the  caret  will  now  be  dropped 
from  the  perturbation  variables.  In  addition,  the  lowest  order  ex¬ 
ponential  basic  state  interfacial  depth,  6^°\  will  be  denoted  by  t;  i.e. 

T(£)  =  6(0)(O  ,  (6.12) 

where  6^  ^  is  defined  in  (5.16).  We  now  assume  time  and  space  harmonic 
solutions,  as  in  Section  3.3,  but  only  for  the  lowest  order  systems 
(6.7)  and  (6.9).  Harmonic  solutions  to  the  first  order  systems  (6.8) 
and  (6.10)  will  be  considered  later  and  boundary  conditions  will  be 
addressed  in  the  next  chapter. 

Thus,  we  let 


u(0)  =  X(c) 

exp{i(c+u>t) } 

(6.13a) 

0(0)  -  Y(C) 

exp{i(c+u)t) } 

(6.13b) 

«(0)  -  z(£) 

exp{i(£+wt) }  , 

(6.13c) 
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where  X,  Y  and  Z  are  the  (complex-valued)  non-dimensional  perturbation 
amplitudes  for  the  x-  and  y-velocities ,  and  the  interfacial  depth, 
respectively.  Note  that  the  time  dependence  assumed  in  Chapter  3 
differs  in  sign  from  the  one  assumed  here,  cu  is  the  dimensionless 
circular  frequency  of  the  disturbance,  which  is  related  to  the  dimen¬ 
sional  frequency,  0,  by 

(A  *■  fl/f  .  (6.14) 

Using  (6.13)  in  (6.7)  and  (6.9)  reduces  the  system  of  PDE's  in 
each  case  to  the  following  systems  of  ODE's  for  the  ^-dependent 
perturbation  amplitudes : 


iaX  -  Y  +  *  0 

(6.15a) 

ioY  +  X  +  ikXZ  =  0 

(6.15b) 

iaZ  +  ikX(l-Bv  E)Y  -  8v  X  +  Cl-Bv  E)XP  =  0 

00  00  00  V 

(6.15c) 

(PLANAR  FRONT,  ORDER  F°) 

b 

ig(C)X  -  Y  +  Z  *  0 

(6.16a) 

ig(C)Y  +  tX  +  ikXZ  =■  0 

(6.16b) 

ig(C)Z  +  ikXxY  +  (1-t)X  +•  tX_  »  0 

(6.16b) 

(EXPONENTIAL  FRONT,  ORDER  F°) 

b 
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in  which 


(6.17) 


a  *  10  +  kA(l-6)v 

00 

and 

g(5)  -  a)  +  kA{l+v  -t(£)}  .  (6.18) 

oo 

For  the  planar  front,  o  is  the  non-dimensional  Doppler  shifted  fre¬ 
quency,  and  g(£)  is  the  exponential  front  analogue  of  o.  The  subscript 
C  denotes  a  total  derivative  with  respect  to  5. 

6.5  ZEROTH  ORDER  ODE’s  FOR  THE  DEPTH  PERTURBATION  AMPLITUDE 

A  single  equation  for  the  depth  perturbation  amplitude,  Z(£),  is 
derived  from  (6.15)  for  the  planar  front  and  from  (6.16)  for  the  ex¬ 
ponential  front  by  eliminating  the  velocity  perturbation  amplitudes, 
X(£)  and  Y(£),  from  these  systems. 

The  planar  front  result,  which  follows  directly,  is 

(8voo5-l)(d2Z/d52)  +  gv^dZ/dO  +  A(£)Z(5)  =  0,  (6.19) 

where 

A(£)  *  1  -  o2  +  kA8v Jo  +  (kA) 2 (l-8v  O.  (6.20) 

(6.19)  is  linear,  homogeneous,  second  order  in  Z(£),  and  independent  of 
uo>  The  cross-stream  flow  does  not  en«-er  the  perturbation  dynamics  nor, 
as  we  will  see,  the  boundary  conditions  to  lowest  order  in  F^  for  the 
planar  front,  and  therefore  contributes  nothing  to  the  system's  sta¬ 
bility  characteristics.  If  (6.19)  is  redimensionalized,  the  equations 
developed  earlier  for  the  depth  perturbation  amplitude,  viz.,  (3.12) 
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As  for  the  planar  front,  the  depth  perturbation  amplitude  equation 
for  the  exponential  front  is  linear,  homogeneous  and  second  order.  We 
again  note  that  the  cross-stream  flow  does  not  enter  the  exponential 
front's  dynamics  to  lowest  order.  Since  we  are  primarily  interested 
in  the  influence  of  the  cross-stream  flow  on  the  frontal  zone's  sta¬ 
bility  characteristics,  it  is  now  clear  that  an  examination  of  the 
zeroth-order  dynamics  alone  is  inadequate.  The  cross-stream  flow  has, 
at  most,  a  first  order  effect  on  stability.  In  subsequent  chapters, 
therefore,  the  first  order  dynamics  will  be  developed,  and  attention 
will  be  focused  on  the  attendant  dispersion  relation. 


CHAPTER  7 


BOUNDARY  CONDITIONS 

7.1  INTRODUCTION 

In  this  chapter,  boundary  conditions  are  formulated  for  the 
lowest  order  perturbation  amplitudes  X(£),  Y(£)  and  Z(£),  and  for 
the  first  order  perturbations  u^}  v^  and  6^  for  both  the  expo¬ 
nential  and  planar  fronts.  Figures  7.1  and  7.2  show  the  simplified 
exponential  and  planar  frontal  zone  geometries,  respectively,  under 
the  influence  of  a  wavelike  disturbance  that  distorts  the  equilibrium 
boundary  separating  the  inner  and  outer  regions.  With  no  disturbance, 
the  surface  boundary  between  these  two  regions  is  contained  in  the 
plane  x-0. 

As  shown  in  Figures  7.1  and  7.2,  the  surface  displacement  of 

the  equilibrium  boundary  is  denoted  h(y,t)  in  dimensional  variables. 

For  small  perturbations,  |h|/X<<l.  The  interfacial  depth  in  the  plane 

x*0  is  5  D,  for  the  exponential  front,  and  D.(0)  for  the  planar  front. 

O  0  I 

The  boundary  condition  on  the  lowest  order  depth  perturbation  ampli¬ 
tude  Z  is  derived  by  invoking  the  kinematic  restriction  that  the  inter- 
facial  depth  beneath  the  displaced  inner-outer  region  boundary  remain 
at  its  equilibrium  value,  a  reasonable  restriction  as  long  as  the  per¬ 
turbation  amplitude  is  small.  In  addition  to  being  plausible,  this 
condition  may  be  formulated  rigorously  by  taking  only  the  lowest  order 
term  in  the  Taylor  series  expansion  for  the  time-varying  interfacial 
depth  about  its  basic  state  value.  The  boundary  condition  on  the 
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FRONTAL  ZONE  GEOMETRY 
FOR  THE  PLANAR  FRONT  WITH 
WAVE  DISTURBANCE 


FIGURE  7.2 


x-velocity  perturbation  is  derived  from  the  kinematic  condition 
that  the  normal  velocity  across  the  boundary  between  the  inner  and 
outer  regions  must  be  continuous.  We  also  require  no  velocity 
shear  across  the  boundary  in  the  long-stream  direction,  which  leads 
to  a  boundary  condition  on  the  y-velocity  perturbation. 

7.2  BOUNDARY  CONDITION  ON  THE  DEPTH  PERTURBATION 

The  kinematic  restriction  that  the  interfacial  depth  at  the 
inner-outer  boundary  remains  at  its  equilibrium  value  may  be  written 

6(h)  +  6{h(c,t),C,t}  =  6eq  ,  (7-l) 

where  6  *6  for  the  exponential  front  and  6  =  1  for  the  planar  front, 

eq  o  eq 

Note  that  (7.1)  is  written  in  dimensionless  variables,  that  is,  with 
h*-h/X  and  subsequently  dropping  the  asterisk.  6(h)  is  the  basic 
state  interfacial  depth  evaluated  at  5**h,  that  is,  at  the  displaced 
surface  boundary.  (7.1)  is  translated  into  boundary  conditions  on 
the  depth  perturbations  to  lowest  and  first  orders  in  F,  by  substitut- 

D 

ing  the  small  parameter  expansions  in  (6.6)  and  invoking  the  approxi¬ 
mations  Z(h)“Z(0)  and  e  ^=(l-h)  for  h«l.  To  lowest  order  in  F,  ,  we 

„  D 

obtain 


h  *  h  exp{i(5+ujt)  } 
o 


(7.2a) 


where 


h  -  Z(0) / (6  -1) 
o  o 


(7.2b) 
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for  the  exponential  front,  and 

h  -  Z(O)/0v  (7.2c) 

o  00 

for  the  planar  front.  At  first  order  In  F ,  we  obtain 

b 

5(1)(h(C,t),?,t}  -  0  (7.3) 

for  both  the  planar  and  exponential  fronts.  Condition  (7.3)  re¬ 
quires  that  the  first  order  depth  perturbation  amplitude  vanish  in 
both  cases  at  the  inner-outer  boundary. 

7.3  VELOCITY  BOUNDARY  CONDITIONS 

If  the  small  parameter  expansions  in  (6.6)  are  substituted 
into  the  perturbation  variable  decomposition  introduced  in  Section 
6.2,  the  lowest  order  flow  field  for  5>h  in  the  exponential  front 
case  becomes 


u(5,C,t)  -  l/S(0)(O  +  Q(0)a,H, t) 


(7.4a) 


v(£,C,t)  -  1  -  6<0)U)  +  Vco  +  v(0)(C,C,t)  (7.4b) 


To  first  order  in  F  ,  we  have 

D 
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(7.5) 


as  a  consequence  of  only  even  powers  of  F,  appearing  in  the  small 

b 

parameter  expansion  of  the  exponential  basic  state  interfacial  depth. 
Condition  (7.5)  is  found  to  apply  also  for  the  planar  front,  since 
there  the  equilibrium  flow  is  spatially  invariant.  At  first  order 
in  F,  ,  then,  the  velocity  perturbation  amplitudes  in  the  x-  and 

D 

y-directions  must  be  zero  for  both  the  planar  and  exponential  fronts 
At  zeroth  order,  of  course,  the  velocity  perturbations  are  not  zero, 
as  discussed  next. 

The  restriction  that  the  normal  velocity  at  the  boundary  be¬ 
tween  the  inner  and  outer  regions  be  continuous  may  be  written 


Lim  jjj;  (x-h)  =  G(0)(0)  (7.6) 

x-*-h 

in  dimensional  variables,  where  D/Dt  is  the  convective  derivative 
and  u(0)(0)  is  the  equilibrium  x-velocity  evaluated  at  the  boundary. 
Expanding  (7.6)  to  lowest  order  in  F.  yields 

D 

U(0)(h)  +  U(0)(h)  -  (v(0) (h)  +  v(0)(h)}(3h/3y) 

-  rah/at'*  - 


86 


in  which  all  velocities  on  the  left -hand-side  are  evaluated  at  x»h, 
i.e. ,  at  the  boundary.  The  basic  state  terms  involving  u  are  elimin¬ 
ated  by  the  approximation  u^  (h)au^  (0)  for  |h|/X<<l,  and  the  sec¬ 
ond  order  term  v^(h)3h/3y  may  be  neglected  altogether,  giving 

fl(0)(h)  -  v(0)(h)(3h/3y)  -  (3h/3t)  -  0  .  (7'7) 

(7.7)  consitutes  the  fundamental  zeroth  order  boundary  condition  on 
the  x-velocity  perturbation  amplitude.  Rewriting  (7.7)  in  non- 
dimensional  variables  in  accordance  with  (6.1),  substituting  the 
explicit  form  for  h(C,t)  in  (7.2a,b),  approximating  X(h)  by  X(0)  for 
|h|<<l,  and  neglecting  terms  of  order  h  compared  to  1  yields  the 
following  result: 


X(0)  -  ih  (u>  +  kX(v  +6  -1)} 

O  •  0 


(7.8) 


(7.8),  written  in  scaled  variables,  provides  the  connection  between 
the  depth  perturbation  amplitude  at  the  origin,  Z(0),  and  the  x- 
velocity  perturbation  amplitude  at  the  origin,  X(0) ,  through  the 
displacement  amplitude  hQ  appearing  in  (7.2b).  We  note  that,  for 
a  fixed  disturbance  frequency  and  wavelength,  X(0)  and  Z(0)  are  direct¬ 
ly  proportional,  but  in  phase  quadrature. 

For  the  planar  front,  we  follow  exactly  the  same  procedure 

and  arrive  at  the  following  boundary  condition  on  the  (dimensionless) 

x-velocity  perturbation  amplitude  to  lowest  order  in  F,  : 

b 


r 


i 
i 

n 


•  ^  Y*>  '•-Tf  ■•  ffc  •>  1  t' 
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X(0)  =  ih  {u>  +  kX(l-6)v  } 
o  00 


(7.9) 


Again  to  lowest  order  in  F,  ,  we  require  for  the  long-stream 

D 


velocity  that  no  velocity  shear  exist  across  the  inner-outer  region 
boundary,  although  a  long-stream  velocity  shear  can  exist  between 
the  frontal  zone  and  the  ambient  ocean.  This  restriction  is  physi¬ 
cally  realistic  for  small  amplitude  wave  disturbances,  and  may  be 
justified  rigorously  by  performing  a  Taylor  series  expansion  of  the 
time-varying  y-velocity  about  its  equilibrium  value  and  retaining 


only  the  leading  term.  Therefore,  as  a  boundary  condition  on  Y(5) 


at  order  F^,  we  require 


Y  (0)  -  0 


(7.10) 


for  both  the  planar  and  exponential  fronts. 
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CHAPTER  8 


ANALYSIS  OF  THE  EXPONENTIAL  FRONT  EQUATION  FOR  Z(H) 

TO  LOWEST  ORDER 

8.1  INTRODUCTION 

In  this  chapter,  the  exponential  front  depth  perturbation 
amplitude  equation  (6.23)  will  be  analyzed  in  some  detail.  Its 
singular  points  in  x-space  and  n-space  are  located,  the  asymp¬ 
totic  form  of  the  solution  far  from  the  surface  front  is  pre¬ 
sented,  and  the  condition  for  stable  solutions  at  infinity  is 
developed.  The  equation  for  Z(q)  is  then  integrated  numerically, 
and  several  of  the  resulting  mode  shapes  are  discussed.  Recall 
that  the  analysis  of  the  planar  front  DPA  equation,  (6.19),  was 
performed  in  Chapters  3  and  4. 

8.2  SINGULAR  POINTS  OF  THE  EQUATION  FOR  Z(n) 

Singularities  in  the  depth  perturbation  amplitude,  Z(n), 

occur  at  most  at  the  zeros  of  the  polynomial  A(n)  in  (6.24a).  Since 
A(n)  is  presented  in  factored  form,  its  zeros  are  readily  determined 
and  appear  in  Table  8.1.  Recall  that  the  domain  of  definition  of 
the  original  independent  variable  £  is  0<££“.  Transforming  to  the 
Independent  variable  t  maps  the  interval  0<5<®  into  the  interval 
60<t<1,  while  the  transformation  from  x  to  n  in  (6.23)  maps  the  do¬ 
main  of  definition  for  x  into  an  interval  in  n-space  with  variable 
endpoints  whose  locations  depend  on  the  values  of  ui,  k\  and  V^. 
Since,  in  general,  the  locations  of  the  singular  points  also  depend 
upon  the  values  of  w  ,  <\  and  Vm  ,  a  determination  of  whether  or  not 
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a  particular  singular  point  lies  in  the  domain  of  definition  of 
the  independent  variable  is  most  easily  made  when  the  domain  has 
fixed  endpoints.  For  this  reason.  Table  8.1  includes  the  locations 
of  the  singular  points  in  both  rv-space  and  x-space.  We  will  here¬ 
after  deal  almost  exclusively  with  the  t  variable,  since  its  domain 
of  definition  is  fixed. 

If  the  parameters  u,  kA  and  are  such  that  the  x^j*x. 

Table  8.1  for  ii* j  ,  that  is,  if  the  singular  points  are  distinct, 
then  the  ratios  B(n)/A(n)  and  C(n)/A(n)  in  (6.24)  possess  poles  of, 
at  most,  first  and  second  order,  respectively.  Each  singular  point 
is  therefore  regular  if  the  are  distinct.  Irregular  singularit¬ 
ies,  and  the  attendant  analytical  difficulties,  result  from  the  co¬ 
alescence  of  two  or  more  of  the  singular  points  in  Table  8.1,  which 
can  occur  for  certain  values  of  u,  k\  and  V  .  For  example,  if  these 
quantities  are  such  that 

l/(k*)2  *  -4<w/kA  +  vm  +  1)  , 


the  two  otherwise  distinct  regular  singular  points  t ^  and  x^  coalesce 
into  a  single  singular  point  which  is  Irregular.  Since  the  develop¬ 
ment  of  analytical  solutions  near  Irregular  singular  points  is,  in 
general,  very  difficult,  we  will  require  hereafter  that  all  five  of 
the  x^  are  distinct,  thereby  Insuring  regular  singular  points  for 
the  depth  perturbation  amplitude  ODE. 

Physically  significant  singular  behavior  in  the  solution  Z 

must,  of  course,  occur  in  the  interval  0<£<«,  or,  equivalently,  6  <x<l. 

o 

Singular  behavior  in  Z  for  values  of  the  Independent  variable  out- 
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side  this  range  Is  physically  meaningless  and  consequently  of  no 
concern.  The  fixed  singular  point  tq*0  Is  thus  irrelevant,  while 
the  fixed  singular  point  x^«l  Is  of  considerable  importance,  since 
it  corresponds  to  E*®.  It  thereby  determines  the  nature  of  the  sol¬ 
ution  at  infinity  in  the  original  x-coordinate  system.  The  ODE  for 
Z  thus  admits  the  possibility  of  a  solution  that  "blows  up"  at  infin¬ 
ity  for  some  values  of  the  parameters  o>,  k\  and  Va.  At  the  remain¬ 
ing  three  singular  points,  t^,  x^,  and  x^»  divergent  solutions  for 
Z  are  also  possible,  but  these  points  may  be  placed  outside  the  phys¬ 
ically  significant  interval  by  suitably  choosing  10,  < A  and  V^.  The 
locations  of  x^,  x^  and  x^  in  the  kA-w  plane  are  shown  in  Figure 

8.1.  The  nature  of  the  solution  near  x  »1  will  be  examined  in  de- 

4 

tail  in  the  next  section. 

8.3  BEHAVIOR  OF  THE  SOLUTION  NEAR  INFINITY 

Since  physically  acceptable  solutions  must  exhibit  regular 
behavior  at  Infinity,  we  wish  to  Investigate  the  nature  of  the  sol¬ 
ution  to  (6.23)  as  or,  equivalently,  y-*>l  and  n-*icA(u-l)  »  where 

u  *  u>/kA  +  v^  +  1  .  (8.1) 

We  also  require  ui*l,  since  when  u_l  the  two  singular  points  x^  and 
x^  coalsece  into  a  single  higher  order  singularity  at  infinity. 

With 

s  -  1  +  r\/k\  -  y  ,  (g.2) 

the  problem  reduces  to  an  examination  of  the  solution  to 
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s  Z  +  sZ  +  (d/a)Z(s)  -  0  (8. 

ss  s 

as  s->0  from  above,  where 

a  -  (l-u)2(l  -(kA)2(l-ii)2}  (8. 

and 

d  -  (k\)2(l-u)2  +  {(kX)2(l-u)2  -  1}2{1-(1-U)2}  -  1  .  (8- 

Under  the  transformation  s-e1",  we  obtain 

Z  +  (d/a)Z(t)  -  0  ,  (8. 

whose  solution  is  written  immediately  as 

Z  *  A  exp{ (-d/aj^t)  +  B  expf-C-d/aJ^t}  ,  (8. 

where  A  and  B  are  arbitrary  constants.  Rewriting  (8.7)  in  terms 
of  the  original  £ -coordinate  yields 
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ating  the  domain  of  allowable  (tcX.w)  values  from  that  in  which  di¬ 
vergent  solutions  exist.  Since  the  disturbance  wavelength  on  a  large 
scale  oceanic  front  is  generally  larger  than  the  Rossby  radius,  that 
is,  kA<1  typically,  (8.10)  and  Figure  8.2  show  that  the  disturbance 
periods  will  generally  be  longer  than  the  inertial  period,  or,  equiva¬ 
lently,  uj<1.  In  other  words,  we  anticipate  sub-inertial  frequencies. 

In  Table  8.2,  the  exact  dispersion  relation  for  the  exponential 
front  at  lowest  order  in  F.  ,  which  is  developed  analytically  in  Sec- 

D 

tion  10.1,  has  been  used  to  tabulate  the  disturbance  frequency  for 
various  combinations  cf  the  parameters  <5q  and  as  a  function  of  the 
normalized  wavenumber  kA.  ranges  from  0.1  to  0.5  in  steps  of  0.1, 
as  does  <A  ranges  from  0.1  to  1.9  in  steps  of  0.1.  Entries  carry¬ 

ing  a  superscript  "a"  correspond  to  values  of  u>,  <A  and  which  vio¬ 
late  the  stability  criterion  at  infinity  given  in  (8.9).  Entries 
carrying  a  superscript  "b"  correspond  to  parameter  values  that  place 
one  or  more  of  the  singular  points  of  the  DPA  equation  in  the  inter¬ 
val  6  <t<1,  as  discussed  in  Section  8.2.  An  examination  of  the  tabu- 
o 

lated  results  again  clearly  shows  that  the  system  is  sub-inertially 
dominated.  In  no  case  is  an  allowable  value  of  the  disturbance  freq¬ 


uency  greater  than  1. 
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Figure  0.2(b) 
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Figure  8.2(e) 


103 


Case  1:  V  »  0.1 
00 


kX 

0.1 

0.2 

0.3 

0.4 

0.5 

0.1 

0.98b 

0. 47b 

0.29 

0.20 

0.14 

0.3 

2. 94a 

1.413 

0.88b 

0. 60b 

0.42 

0.5 

4. 90a 

2.35a 

1.47a 

1.00 

0. 70b 

0.7 

6. 86a 

3.  29a 

2. 05a 

1. 40a 

0.98 

0.9 

8.82a 

4.23a 

2. 64a 

1. 80a 

1. 26a 

1.1 

10. 78a 

5. 17a 

3.23a 

2. 20a 

1.54a 

1.3 

12. 74a 

6.  lla 

3. 81a 

2. 60a 

1. 82a 

1.5 

14. 70a 

7.05a 

4. 40a 

3.00a 

2.10a 

1.7 

16. 66a 

7.99a 

4 . 99a 

3.40a 

2. 38a 

1.9 

18.62a 

8.93a 

5.57a 

3. 80a 

2. 66a 

Case  2: 

V  =  0.2 
00 

k\ 

0.1 

0.2 

0.3 

0.4 

0.5 

0.1 

0.97b 

0.46b 

0.28 

0.19 

0.13 

0.3 

2.91a 

1.38a 

0.85b 

0.57b 

0.39 

0.5 

4.85a 

2. 30a 

1. 42a 

0.95 

0.65b 

0.7 

6. 79a 

3. 22a 

1. 98a 

1.33a 

0.91 

0.9 

8.73a 

4.14a 

2. 55a 

1. 71a 

1. 17a 

1.1 

10. 67a 

5. 06a 

3. 12a 

2. 09a 

1.43a 

1.3 

12. 61a 

5. 98a 

3. 68a 

2.47a 

1. 69a 

1.5 

14.55a 

6. 90a 

4. 25a 

2. 85a 

1.9ba 

1.7 

16. 49a 

7.82a 

4.85a 

3.23a 

2.21a 

1.9 

18.433 

8. 74a 

5. 38a 

3. 61a 

2.27  a 

TABLE  8.2  Computed 

Values 

3 

4-1 

o 

for  Various  Values  of 

kX,  6 

0 

and  V 

00 

I 
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Case  3: 

V  -  0.3 

-J 

CO 

«  X 

0.1 

0.2 

0.3 

0.4 

0.5 

-f 

0.1 

0.96b 

0.45b 

0.27 

0.18 

0.12 

•A 

0.3 

2.88a 

1.35a 

0.82b 

0. 54b 

0.36 

0.5 

4.80a 

2.25a 

1.37® 

0.90 

0. 60b 

•  • 

0.7 

6. 72a 

3.15a 

1.91a 

1. 26a 

0.84 

0.9 

8. 64a 

4. 05a 

2.46a 

1. 62a 

1.08a 

mm 

1.1 

10.56a 

4.95a 

3.01a 

1.98a 

1. 32a 

1.3 

12.48® 

5.85a 

3.55a 

2.34a 

1.56a 

1.5 

14.40a 

6.75a 

4.10a 

2 . 70a 

1.80a 

1.7 

16. 32a 

7.65a 

4. 65a 

3. 06a 

2.04a 

1.9 

18. 24a 

8. 55a 

5. 19a 

3.42a 

2. 28a 

Case  4: 

V  =  0.4 

00 

x  X 

0.1 

0.2 

0.3 

0.4 

0.5 

— 

— 

■  "  1 

0.1 

0.95b 

0.44b 

0.26 

0.17 

0.11 

0.3 

2.85a 

1. 32a 

0.79b 

0. 51b 

0.33 

0.5 

4. 75a 

2.20a 

1. 32a 

0.85 

0.55b 

0.7 

6.65a 

3.08a 

1.84a 

1. 19a 

0.77 

0.9 

8.55a 

3.96a 

2.37a 

1.53a 

0.993 

1.1 

10.45a 

4.84a 

2. 90a 

1.873 

1. 21® 

1.3 

12.353 

5. 72a 

3.42a 

2. 21a 

1.43® 

1.5 

14.253 

6. 60a 

3. 95a 

2.55a 

1.65® 

1.7 

16. 15a 

7.48a 

4.48a 

2.89a 

1.87a 

1.9 

18.053 

8. 36a 

5.00a 

3.23a 

2.09® 
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Case  5: 

V  “  0.5 

00 

«  X 

0.1 

0.2 

0.3 

0.4 

0.5 

0.1 

0. 94b 

0.43b 

0.25 

0.16 

0.10 

0.3 

2.82a 

1. 29a 

0.76b 

0.48b 

0.30 

0.5 

4 . 70a 

2.15a 

1.27a 

0.80 

0.50 

0.7 

6.58a 

3.01a 

1. 77a 

1. 12a 

0.70 

0.9 

8.46a 

3.87a 

2.28a 

1.44a 

0.90' 

1.1 

10. 34a 

4.73a 

2.79a 

1.76a 

1.10' 

1.3 

12.223 

5. 59a 

3.29a 

2 . 08a 

1.30 

1.5 

14 . 10a 

6.45a 

3.80a 

2 . 40a 

1.50 

1.7 

15.98a 

7.31a 

4.313 

2 . 72a 

1.70 

1.9 

17 . 86a 

8. 17a 

4.8la 

3 . 04a 

1.90 

NOTES : 

a  -  the  values 

of  icl,  <5q  and 

V  are 

00 

such  that  the 

stability 

criterion 

at  infinity,  (8 

•9),  is 

violated 

b  -  parameter  values  are  such  that  one  or  more  singular 
points  are  contained  in  the  interval  <5Q<T  _<]_  : 

see  discussion  in  Section  8.2. 
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8.4  TYPICAL  MODE  SHAPES  FOR  THE  ZEROTH  ORDER  EXPONENTIAL 

FRONT  DPA  EQUATION 

Figure  8.3  Illustrates  typical  mode  shapes  for  the  depth  per¬ 
turbation  amplitude  Z(n).  The  curves  were  generated  by  numerically 
integrating  the  lowest  order  exponential  front  DPA  equation,  (6.23), 
using  a  finite  difference  method  that  reduces  the  ODE  to  an  alge¬ 
braic  system  which  is  then  solved  by  matrix  methods  (Conte  and  deBoor, 
1972,  Ch.  5).  The  algorithm's  accuracy  was  checked  against  a  sec¬ 
ond  order  ODE  with  known  analytic  solution.  The  curves  in  Figure  8.3 
were  generated  using  380  steps  in  the  numerical  integration,  and  the 
accuracy  is  expected  to  be  well  within  1%.  Values  of  the  input  par¬ 
ameters  and  «X  were  chosen  so  that  the  stability  criterion  at  in¬ 
finity  in  (8.9)  was  met  and,  additionally,  so  that  none  of  the  sin¬ 
gular  points  t. (i*l,3)  fell  in  the  interval  6  <t,<1.  These  restric- 
l  o  “  l” 

tlons  insure  regular  behavior  of  Z(q)  throughout  its  interval  of 
definition. 

The  value  of  Z  at  5*0  was  arbitrarily  set  equal  to  -0.05,  -0.07 
and  -0.06  in  Figures  8.3(a)  through  8.3(c),  respectively,  which  is 
consistent  with  the  assumption  of  small  amplitude  perturbations.  The 
points  nQ  and  in  these  figures  are  the  images  of  £*0  and  un¬ 
der  the  transformations  defined  in  (6.12)  and  (6.21d).  At  £»»  the 
boundary  condition  is  Z-0  as  demanded  by  the  asymptotic  behavior  of 
the  solution  for  Z(n)  developed  Sai  Section  8.3.  Since  the  DPA  equat¬ 
ion  (6.23)  involves  u  as  a  parameter,  specification  of  an  appropriate 
(kX,<i>)  pair,  that  is,  one  satisfying  the  lowest  order  exponential  front 
dispersion  relation,  is  required  in  order  to  perform  the  numerical 
integration.  This  was  done  by  using  the  results  of  Section  10.1,  in 


which  the  required  dispersion  relation  is  derived  explicitly.  The 

three  curves  in  Figure  8.3  are  representative  of  the  results  thus 

obtained,  and  serve  to  illustrate  typical  behavior  of  Z(n)  for 

values  of  V  ,  6  and  k\  less  than  1  which  meet  the  conditions  dis- 
00  o 

cussed  above.  The  solutions  are  smooth,  as  one  would  expect  from 


the  analysis  above. 


.500  V  =  0.300  u  =  0.840  k\  =  0.700  n  =  1.400 


CHAPTER  9 


DISPERSION  CHARACTERISTICS  OF 
THE  PLANAR  FRONT 

9.1  INTRODUCTION 

Duxbury's  (1963)  frontal  zone  geometry  is  shown  in  Figure 
9.1.  Unlike  the  model  in  this  paper,  his  frontal  zone  consists  of 
two  lnvlscid  regions  of  the  same  density  separated  by  a  vertical 
"vortex  sheet"  that  supports  a  shear  flow  parallel  to  the  boun¬ 
dary  between  these  regions.  In  Region  1,  which  is  analogous  to 
the  wedge  shaped  frontal  zone  modelled  in  this  paper,  a  uniform, 
geostrophic  equilibrium  flow  exists  whose  single  velocity  component 
is  parallel  to  the  vortex  sheet  boundary.  The  circulation  is  zero 
in  the  constant-depth  appendage  which  is  defined  for  y>0.  Duxbury 
develops  the  dispersion  relation  for  this  system,  which  functionally 
relates  the  disturbance  frequency  and  wavelength,  by  solving  the 
time-varying  perturbation  equations  in  Regions  1  and  2,  and  matching 
the  solutions  across  the  vortex  sheet  boundary.  The  result,  which 
appears  as  his  equation  (75),  thus  requires  a  priori  knowledge  of 
the  time-dependent  perturbation  equation  solutions,  since  these  then 
appear  explicitly  in  the  dispersion  relation. 

The  dispersion  characteristics  of  Duxbury's  system  rest  upon 
boundary  conditions  that  are  quite  different  from  those  in  this  model, 
which  consequently  leads  to  a  different  formalism  in  developing  the 
dispersion  relation.  It  will  be  seen  in  what  follows  that  the  dis- 
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persion  characteristics  of  both  the  exponential  and  planar  fronts 
are  available  directly  from  the  perturbation  equations  upon  appli¬ 
cation  of  the  boundary  conditions  developed  in  Chapter  7.  Since  no 
&  priori  knowledge  of  the  perturbation  solutions  is  required,  con¬ 
siderable  simplification  results. 

The  fundamental  nature  of  the  difference  between  Duxbury's 
model  and  this  one  is  readily  appreciated  by  more  closely  examining 
his  dispersion  relation  in  the  limiting  case  Dq-*-0.  Since  Dq  is  the 
(constant)  depth  in  his  Region  2,  letting  Dq-*0  removes  Region  2  from 
the  problem  and  presumably  recovers  the  planar  frontal  zone  geometry 
in  this  model  with  the  inner  region  width  set  to  zero.  We  expect  that 
this  is  a  physically  realistic  limit  which  should  correspond  to  well- 
behaved  dispersion  characteristics.  Taking  the  limit  Dq-K>  in  Duxbury's 
equations  (75)— (77) ,  however,  reveals  that  his  dispersion  relation 
exhibits  no  regular  behavior  in  this  limit  and,  in  fact,  becomes  un¬ 
defined.  This  peculiar  behavior  is  a  consequence  of  the  boundary  con¬ 
ditions  that  are  imposed,  which  require  a  coupling  of  the  dynamics  in 
Regions  1  and  2.  Removal  of  Region  2  therefore  leads  to  a  singular 
result . 

The  model  developed  in  this  paper,  however,  does  not  suffer 
this  shortcoming,  since  the  boundary  conditions  are  different.  The 
frontal  zone  geometry,  moreover,  is  formulated  more  realistically  at 
the  outset,  and  the  artifice  of  a  vortex  sheet  barrier  separating  two 
regions  of  the  same  density  is  unnecessary.  Garvine's  work  (1974, 
1979a,  1979b)  has  shown  clearly  that  the  steep  interface  between  the 
frontal  zone  and  the  ambient  ocean  is  maintained  by  local  dissipation, 
viz.,  turbulent  friction,  mixing  and  mass  entrainment  in  the  inner 
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dissipative  region,  which  ellnlnates  the  need  for  an  artificial  vert¬ 
ical  barrier,  as  in  Duxbury's  model.  In  the  following  two  sections, 
the  zeroth  and  first  order  dispersion  characteristics  are  derived  for 
the  planar  front. 

9.2  THE  ZEROTH  ORDER  DISPERSION  RELATION 

Derivation  of  the  planar  front  dispersion  relation  at  order 
F®  is  remarkably  simple  and  follows  Immediately  by  applying  the  boun- 

D 

dary  conditions  formulated  in  Chapter  7  to  the  lowest  order  pertur¬ 
bation  equations  in  Chapter  6.  Specifically,  we  substitute  (7.10) 
and  (7.9)  along  with  (7.2c)  into  the  y-momentum  equation  (6.15b)  eval¬ 
uated  at  e-0,  which  imposes  the  following  characteristic  condition  on 
the  scaled  frequency  o>  and  the  normalized  wavenumber  icA : 


w/kA  •  -V  (9.1) 

m 

Thus,  at  order  F^,  wavelike  disturbances  on  the  planar  front 

D 

in  the  moving  coordinate  system  are  non-dispersive  with  equal  phase 
and  group  velocities.  For  real  kA,  the  frequency  is  also  real,  and 
therefore  no  unstable  modes  exist  in  the  system.  Since  the  planar 
front  is  inherently  stable  and  non-dispersive  to  lowest  order,  it 
cannot  support  a  disturbance  whose  amplitude  grows  in  time.  In  the 
next  section,  we  will  therefore  consider  the  first  order  dispersion 
characteristics  for  the  planar  front  and  are  led  to  a  rather  surpris¬ 
ing  conclusion. 
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9.3  THE  FIRST  ORDER  DISPERSION  CHARACTERISTICS 

Recall  that  the  perturbation  variable  decomposition  was  intro¬ 
duced  in  Chapter  6,  and  that  the  zeroth  and  first  order  perturbation 
equations  were  derived.  The  lowest  order  systems  appear  in  (6.7)  and 
(6.9)  for  the  planar  and  exponential  fronts,  respectively,  while  the 
corresponding  first  order  systems  are  contained  in  (6.8)  and  (6.10). 
We  then  assumed  the  time-  and  space-harmonic  solutions  in  (6.13)  for 
the  zeroth  order  perturbations,  but  proceeded  no  further  with  the 
first  order  systems.  In  this  section,  we  will  introduce  time-  and 
space-harmonic  solutions  for  the  first  order  velocity  and  interfacial 
depth  perturbations  as  follows: 


a(1)(5,C,t)  -  X(1)(D  expUU-hot)} 


(9.2a) 


0(1)(5,C,t)  -  Y(1)(0  exp{i(c+u>t)} 


(9.2b) 


«(1)(5,C,t)  -  Z(1)(C)  exp{i(C-hi)t) } 


(9.2c) 


The  circular  frequency  u>  appearing  in  (9.2)  is  the  same  as  the  zeroth 
order  <i>  in  Chapter  6,  since  the  inhomogeneous  first  order  system  is 


driven  by  the  lowest  order  solutions. 

Substituting  (9.2)  into  (6.8)  leads  to  the  following  system 
of  ODE's  for  the  first  order  perturbation  amplitudes  X^,  Y^  and 
each  of  which  is  a  function  of  £  alone: 
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ioX(1)  -  Y(1)  +  Z*1*  -  -  uX? 


(9.3a) 


ioY(1)  +  X(1)  +  ikXZ(1*  -  -uY, 


(9.3b) 


ioZ(1)  +  Z5X(1)  +  (1-6VmC)X^1) 

+  ikX(l-8V  C)Y(1)  (9.3c) 

00 

-  -uZ_  -  (6u  /kX)CX_  -  Su  (l+iC)Y 

►  °°  r  oo 


where  a  is  defined  in  (6.17). 

Following  the  notation  introduced  in  Chapter  6,  X,  Y  and  Z 
with  no  superscript  are  lowest  order  quantities,  while  the  correspond¬ 
ing  first  order  quantities  are  marked  by  a  superscript  (1).  As 
before,  the  overbar  denotes  an  equilibrium  state  variable.  We  note 

that  the  equilibrium  state  and  the  solutions  to  order  F^  appear  both 

0 

as  driving  terms  and  as  non-constant  coefficients  in  the  first  order 
system.  The  first  order  dispersion  relation  may  now  be  derived  by 
applying  the  boundary  conditions  on  X^\  Y^  and  Z^  directly  to 
the  system  (9.3).  Doing  so,  however,  is  somewhat  involved  algebraic¬ 
ally  for  the  following  reason.  If  each  of  the  functions  in  (9.3)  is 
evaluated  at  5*0,  the  point  at  which  the  boundary  conditions  are 
applied,  we  see  that  (9.3)  reduces  to  a  system  of  three  algebraic 
equations  for  five  unknowns,  viz.,  X^  ,  Y^,  Z^\  X^  and  Z^ , 
each  evaluated  at  £*0.  The  lowest  order  solutions  and  their  deriva- 
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Cives  evaluated  at  £*0  are  presumed  known.  In  order  to  successfully 
develop  a  characteristic  relation  between  u>  and  <A ,  therefore,  we 
must  first  manipulate  the  system  in  (9.3)  to  generate  a  number  of 
equations  equal  to  the  number  of  unknowns.  This  procedure  is  unfor¬ 
tunately  tedious.  It  is  nevertheless  outlined  here,  since  it  will 
be  used  again  for  the  exponential  front  at  order  F* 

D 

We  begin  by  restating  the  known  zeroth  and  first  order  boundary 
conditions : 


Z(0)  -  BV  h  (7.2c) 

oo  O 

X(0)  -  iho{aH-kA(l-8)Voa}  (7.9) 

Y(0)  -  0  (7.10) 

Z(1)  =  0  (  with  Z(1)(h)“  Z(1)(0)  )  (7.3) 

X(1)  (0)  *  Y(1)(0)  =*  0  ,  (7.5) 


where  hQ  is  the  wave  disturbance  amplitude.  Since  the  lowest  order 

perturbations  are  homogeneous,  vanishing  with  h  ,  and  the  first  order 

o 

solutions  are  driven  by  the  lowest  order  solutions,  h  should  fac- 

o 

tor  out  of  the  dispersion  relation.  We  note  that  the  perturbation 
equations  at  order  appear  in  (6.15)  and  are  not  reproduced  here. 

Evaluating  (6.15a)  and  (6.15c)  at  OO  using  (7.2c),  (7.9)  and 
(7.10)  provides  explicit  representations  for  (0)  and  Z^(0).  (9.3b) 

evaluated  at  £-0  subject  to  (7.3)  and  (7.5)  gives  Y^(0)-0.  Thus, 
each  of  the  lowest  order  perturbation  amplitudes  and  its  first  der- 
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ivative  is  known  explicitly  at  the  origin.  Taking  d/d5  of  (6.15a) 
determines  Z^(0)  as 

Z^(0)  -  -ioX^CO)  ,  (9.4) 

where  X^(0)  is  known.  Differentiating  (6.15c)  once  with  respect  to 
5  and  (6.15b)  twice  provides  the  following  expressions  for  X^(0) 
and  Y^(0)  in  terms  of  previously  computed  quantities: 

X  (0)  -  -(2V  /kA)X.(O)  -  ioZ  (0)  (9.5) 

55  “5  5 

Y^(0)  *  (i/a){X55(0)  +  ikAZ^O)}  (9.6) 

All  relevant  lowest  order  quantities  are  now  known,  and  we  proceed 
to  close  the  first  order  system  in  a  similar  manner. 

From  (9.3a)  and  (9.3c)  with  5*0,  we  obtain,  respectively, 

Z^1}(0)  *  -Gx?(0)  (9.7) 

X^(0)  -  -uZ^0)  +  (uoB5/(kX)2}X5(0)  ,  (9.8) 

wherein  the  previously  derived  results  have  been  used.  Note  that 
u,  which  is  constant,  is  defined  in  (6.3a).  Finally,  differentiat¬ 
ing  in  order  (9.3b),  (9.3a)  and  (9.3c)  once  with  respect  to  5  and 
setting  5*0  closes  the  algebraic  system  by  providing  five  equations 
for  the  five  unknowns  X^X)  (0) ,  Y^1J  (0)  ,  Z^(0)t  X^  (0)  and  Z^  (0) . 

Note  that,  as  a  consequence  of  the  successive  differentiations,  these 
are  not  the  same  five  unknowns  obtained  by  simply  setting  £*0  in  the 
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original  system  (9.3). 

The  final  result  is  conveniently  written  in  matrix  form 
as  follows: 


■* 

Z5(0)+Ve8/kA 

ikA 

ia 

1 

0 

Y(1)(0) 

xc 

la 

-1 

0 

0 

1 

Y^1}(0) 

1 

io 

ikA 

0 

0 

Z^1}(0) 

1 

0 

0 

0 

0 

x^(0) 

0 

0 

1 

0 

0 

«) 

A 

B 

*  C 

D 

E 

where 

A  -  -GZ10)  +  {u c/(kX)2}X ff(0) 

SC  09  sc 

+  (u  /kX)(l+ic)Y  (0) 
00  SC 

B  -  -Gx^CO) 

C  -  -uYc5(0) 


(9.9) 


(9.10a) 

(9.10b) 

(9.10c) 
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D  -  -Gz^(O)  +  {u-C/(kX)2}X5(0)  (9.lOd) 

E  -  -uX^(O)  (9. lOe) 


With  a  bit  of  algebra,  the  Inverse  of  the  coefficient  matrix 
may  be  derived  and  the  system  In  (9.9)  solved  to  give 


where 

b14“1/a14’  b23“1/a13’  b24"_1/a13a14 

b25""a12/a13a14’  b4l"1/ai4’  b43“‘a12/a13a14 


(9.12) 


(9.9)  with  obvious  equivalences 


The  first  order  dispersion  relation  is  now  derived  directly 


from  (9.9)  and  (9.11).  From  (9.9),  we  have 


Substituting  (9.14)  into  (9.13)  results  in  the  characteristic  relation 


between  u  end  kA  demanded  by  internal  consistency  within  the  original 
system  of  dynamical  equations,  the  result  being 


where  the  a 
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the  are  elements  of  Its  Inverse  defined  In  (9.12).  Expanding 
the  coefficients  In  (9.15)  using  (9.12)  leads  to  the  rather  sur- 
prising  result  that  each  of  the  coefficients  of  C,  D  and  E  In  (9.15) 
vanishes  Identically  and  Independently.  Therefore,  no  dispersion 
relation  exists  to  first  order  for  the  planar  front,  and  the  original 
system  of  perturbation  equations  Is  Internally  self-consistent  with¬ 
out  the  constraint  of  any  characteristic  relation  between  u  and  <X. 

In  view  of  the  lowest  order  dispersion  relation,  which  Is  physically 
uninteresting,  the  first  order  result  developed  here  shows  a  grave 
deficiency  in  the  planar  front  model.  The  planar  front  model,  even 
with  appropriate  boundary  conditions  and  scaling,  is  incapable  of 
reproducing  the  dynamical  behavior  of  waves  on  oceanic  fronts.  In 
the  next  chapter,  however,  the  procedures  outlined  above  are  applied 
to  the  exponential  front  with  far  more  satisfying  results. 


CHAPTER  10 


DISPERSION  CHARACTERISTICS  OF  THE 
EXPONENTIAL  FRONT 

10.1  THE  ZEROTH  ORDER  DISPERSION  RELATION 

The  zeroth  and  first  order  dispersion  relations  for  the 
exponential  front  are  derived  in  the  same  manner  as  were  the  planar 
front  results  in  Chapter  9.  The  exponential  front  equations  for  the 
perturbation  amplitudes  X,  Y  and  Z  at  order  F^  appear  in  (6.16). 
Evaluating  the  y-momentum  equation  (6.16b)  at  5*0  then  gives 

ig(0)Y(0)  +  6qX(0)  +  ikAZ(O)  -  0  .  (10.1) 

Substituting  the  boundary  conditions  (7.10),  (7.8)  and  (7.2b)  into 
(10.1)  and  eliminating  a  common  factor  of  1Z(0)  imposes  the  follow¬ 
ing  characteristic  condition  on  m  and  kA  at  order  F®: 

0 

u»/kA  -  (l/5o)  -  «o-  Vw  .  (10. 2) 


The  phase  and  group  velocities  are  equal  and  constant,  as  they 


were  for  the  planar  front.  And,  since  id  Is  real  for  real  icA,  the 
exponential  front  is  dynamically  stable  at  order  F^.  A  plot  of 

D 


Cph(gr)+V»  aS  a  function  50  appears  in  Figure  10.1,  where  cph(gr) 

is  the  phase  (group)  velocity.  For  V >(1/6  -6  ),  C  .  ,  .<0;  other- 

“  oo  ptHgr; 
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vise,  Che  phase  and  group  velocities  are  positive  in  the  moving 
coordinate  system.  Realistic  values  of  are  not  likely  to  exceed 
approximately  0.5.  For  ^*0.5,  the  phase  (group)  velocity  will  be 
negative  only  if  6q  is  greater  than  about  0.8,  which  is  unrealist¬ 
ically  large.  The  range  of  V  and  <5  likely  to  be  encountered  in 

os  0 

an  oceanic  frontal  zone  is  therefore  such  that  the  phase  and  group 
velocities  in  the  moving  coordinate  system  will  be  positive.  Stable 
lowest  order  waves  on  the  exponential  front  are  expected  to  propa¬ 
gate  in  the  positive  y-direction  in  a  right-handed  coordinate  system 
with  the  x-axls  oriented  perpendicular  to  the  long-stream  direction 
and  into  the  frontal  zone.  This  result  is  consistent  with  observed 
meander  patterns  on  the  Gulf  Stream,  in  which  disturbances  propagate 
northward  and  eastward  along  the  Stream. 

10.2  FIRST  ORDER  DISPERSION  CHARACTERISTICS 

As  before,  time-  and  space-harmonic  first  order  solutions  of 
the  form  (9.2)  are  assumed  and  substituted  into  the  perturbation  var¬ 
iable  equations  (6.10),  which  leads  to  the  following  system  of  ODE's 
for  the  perturbation  amplitudes  X^  ,  Y^  and  Z^  : 

ig(C)X(1)  -  Y(1)  + 

-  ( (1-t)/t2}X  -  (1/t)X? 


(10.3a) 
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ig(C)*(1)  +  tX(1)  +  ikXZ(1) 

-  -(1/t)Y?  (10.3b) 


lg(OZ(1)  +  TX*1*  +  lkXtY(1)  +  (1-t)X(1) 
-  {(1-t)/t2}Z  -  (1/t)Zc 


(10.3c) 


where  the  quantities  T(?)  and  g(c)  are  defined  in  (6.12)  and  (6.18), 
respectively.  The  notation  Introduced  in  Section  9.3  is  also  employed 
here. 

The  boundary  conditions  for  the  exponential  front  are  repro¬ 
duced  here  from  Chapter  7: 


Z(0)  -  (6  -l)h 

0  0 

(7.2b) 

X(0)  -  ih  fui+kA(V  +6  -1)} 
o  °®  o 

(7.8) 

Y (0)  -  0 

(7.10) 

Z(1)(0)  -  0 

(7.3) 
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X(1)(0)  -  Y(1)(0)  -  0  (7.5) 

Once  again,  we  expect  the  disturbance  amplitude,  h  ,  to  factor  out 

o 

of  the  dispersion  relation,  since  the  perturbations  at  order 

b 

are  homogeneous.  Following  the  procedure  outlined  in  Section  9.3, 
the  system  (10.3)  is  closed  by  manipulating  it  to  produce  a  number 
of  equations  equal  to  the  number  of  unknowns. 

The  first  step  is  to  compute  the  relevant  lowest  order  quan¬ 
tities,  which  are  the  first,  second  and  third  derivatives  of  the 
perturbation  amplitudes  X,  Y  and  Z  evaluated  at  £-0.  This  is  done 
iteratively  using  the  zeroth  order  system  of  ODE's  in  (6.16).  Thus, 
from  (6.16a)  using  (7.10)  and  setting  £-0, 

Z?(0)  -  -ig(0)X(0)  .  (10.4) 

Evaluating  (6.16c)  at  £-0  with  (7.10)  gives 


X.(0)  -  -(1/6  ){ig(0)Z(0)+(l-6  )X(0) } 
^  o  o 


(10.5a) 


while  (10.3b)  with  (7.3)  and  (7.5)  at  £*0  yields 


Y?(0)  -  0 


(10.5b) 


Differentiating  in  order  (6.16a)  and  (6.16c)  once  with  respect  to 


5  ,  using  t  -t^^— t^»1-t,  and  aettin8  5*0  gives 

Zc_(0)  -  -i(gc(0)X(0)  +  g(0)Xc(0)}  (10.6) 

X^(0)  -  -(i/«o)(g5(0)Z(0)  +  g(0)Z5(0)> 

(10.7) 

+  {(6 -D/4  HX(0)  -  2x  (0)}  . 

0  0  £ 

Taking  d^/dC^  of  (6.16b)  with  5*0  determines  Y^(0)  as 

Y^(0)  -  (i/g(0) H (6q-1) (X(0)  -  2X^(0)) 

+  6  Xrr(0)  +  ikXZ  (0) } 
o  55  55 

(10.8) 


Computing  *-he  third  derivative  of  Z  from  (6.16a)  and  setting  5*0 


gives 

-  -i{g^(0)X(0)  +  2g^ (O)X^ (0)  + 
5*0 

g(0)X^(0)}  +  Y^(0) 


(10.9) 


2  2 

Taking  d  /d5  of  (6.16c)  with  5-0  and  the  previous  results  yields 


29 


— T  “  -(i/«  Hg„(0)Z(0)  +  2g  (0)Z  (0) 
dg  ^  0  «  ^  5 

4  g-0 


+  g(0)Z  (0)  +  kA6  Y__ (0) } 

44  o  44 


+  { («o-l)/6o>{X(0)  -  3X?(0)  +  3X  (0)} 


(10.10) 


And,  from  (6.16b),  computing  the  third  derivative  with  respect  to 
g  leads  to 


— J  -  (i/g(0)){3ig-(0)Y._(0) 

d£  *  ^ 

4  g-0 


+  (l-«o){X(0)  -  3X^(0)  +  3X  (0)} 


(10.11) 


+  «oXC5C(0)  +  ikAZc^(0)> 


All  relevant  lowest  order  quantities  at  £-0  thus  appear  in  (7.2b), 
(7.8),  (7.10)  and  (10.4)-(10.11) . 

The  system  is  now  closed  at  order  by  working  with  (10.3). 
Substituting  the  above  results  into  (10.3a,c)  and  setting  £-0  yields 


Z5(1)(0)  -  (l/«o)[((l-«o)/6o}X(0)  -  X?(0)J  (10.12) 


X5(1)(0)  -  (1/6*)  f{(l-«J/«JZ(0)  -  Z_ (0)1  . 
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The  quantities  that  remain  co  be  computed  are  the  first  derivative 
of  Y^,  and  the  second  derivatives  of  X^  ,  Y^  and  Z^,  each  eval¬ 
uated  at  the  origin.  This  is  accomplished  by  sequentially  differen¬ 
tiating  (10.3b),  (10.3a)  and  (10.3c)  once  with  respect  to  £  and  then 
(10.2b)  twice.  Putting  £*0  end  employing  the  previous  results  leads 
to  the  following  equations  which  determine,  respectively,  Y^(0), 

Z<X)  (0) ,  X<X)  (0)  and  Y^  (0) : 


“F 


ig(0)Y^1)(0)  +  6  X*X)(0)  +  ikXZ^X)(0) 

5  o  £  t 

(10.14) 

£J 

-  {(1-d  )/«J}Y.(0)  -  (l/«  )Y__(0) 

0  0  £  0  ££ 

ig(0)X<X)(0)  -  Y<X)(0)  +  Z^X)(0) 

-  |  {(1-t)/t2)X  -  (1/t)X? 

(10.15) 

11 

£«o 

ig(0)Z^X)(0)  +  2(1-5  )X^X)(0)  + 

<s  xfX)(0)  +  ikX5 
o  ££ 

Yp1^ (°) 

o  £ 

-  |((1-t)/t2>Z  -  (1/t)Z5 

' 

(10.16) 

£-0 

2ig-(0)Y^1) (0)  +  ig(0)Y^X) (0)  +  2(1-5  )x!X)(0) 

£  s  ££  o  ; 

-I 

+  6oX^X)  (0)  ♦  ikXZ^X)  (0) 

/ 

(10.17) 

.1 

£-0 

i 
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We  have  now  generated  a  closed  algebraic  system  of  six  equat¬ 


ions  for  the  six  unknowns  X^^(O),  Y^(0)t  Z^  (0)  ,  (0)  ,  Y^  (0) 


,(1) 


and  Z^7(0).  It  is  conveniently  written  in  matrix  form  as  follows: 


2(l-«o) 

2ig^(0) 

0 

<5 

0 

ig(0) 

ikA 

X^OJ) 

A 

2(l-40) 

ikA6 

0 

ig(0) 

6 

0 

0 

0 

Y^1} (0) 

B 

ig(0) 

-1 

0 

0 

0 

1 

zlx)(0) 

C 

« 

O 

ig(0) 

ikA 

0 

0 

0 

X 

'w' 

o 

* 

D 

0 

0 

1 

0 

0 

0 

Y^P(0) 

E 

6 

0 

0 

0 

0 

0 

0 

z^(0) 

F 

(10.18) 


where  the  six  constants  A  through  F  on  the  right-hand-side  of  (10.18) 
are  all  known  in  terms  of  the  zeroth  order  perturbation  amplitudes 
evaluated  at  £=0  as  follows : 


k  (<(1- 


t)/t  }Y  -  (1/t)Y 


) 


(10.19a) 


'«( 

'«( 


5=0 


{ (1-t) /x  }Z  -  (1/t)Z. 


) 


(10.19b) 


{ (1-t)/t  }X  -  (1/t)X 


) 


5=0 


(10.19c) 


5=0 


(10.19d) 
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D  -  {(1-6  )/6‘}Y,(0)  -  (1/6  )Yr(r(0) 
004  o  44 


E  -  {(1-6  )/6  ;X(0)  -  (1/6  )X_(0) 
oo  o  4 


(10. 19e) 


F  -  {(1-6  )/6‘}Z(0)  -  (1/ 6  )Z_ (0) 
oo  04 


(10. 19f ) 


a 


Writing  the  elements  of  the  coefficient  matrix  in  (10.18)  as 
aij '  w^ere  standard  row-column  notation  is  employed,  and  the  elements 
of  its  inverse  matrix  as  ,  we  have 


— 

”  - 

x<1}(0) 

0 

0 

0 

0 

0 

b16 

A 

Y^1} (0) 

0 

0 

0 

b24 

b25 

b26 

B 

z<1}(0) 

0 

0 

0 

0 

b35 

0 

C 

x^(0) 

0 

b42 

0 

b44 

b45 

b46 

D 

Y^}(0) 

b51 

b52 

b53 

b54 

b55 

b56 

E 

z^}  (0) 

0 

0 

b63 

b64 

b65 

b66_ 

F 

(10.20a) 
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where 


b5l“1/a15’  b42“1/a24’  b52=_a14/a24a15 


53“_al6/a15a36’ 

b63"1/a36’  b24=b51 

44"'a22/a24a15’ 

b64‘b63b51 

b54"('al6/a36+a22a14/a24'a12)/a15 
b25*b53*  b35*b63 

b45“(a22al6/a15'al5)/(a36a24) 
b55”al6(al6/a36“a22al4/a24+a12+a15a14/a16a24 
b65*b63b53’  b16*b42 ’  b26=b52 

b46’(a22al4/al5'a21)/a24 
b56*_(all/al5“a14a21/,al5a24~a16/a36+(a14/al5 
+a22al4/a24-a12))/a24 
b66"-(al5+al4/al5)/(a24a36) 


)/(a36aI5 


) 


(10.20b) 

)(-al6/a36 


The  dispersion  relation  is  developed  directly  from  (10.18) 
and  (10.20a).  Returning  to  (10.18),  we  have 


a24X£1)(°)  "  F  (10.21a) 

a36Z<1)(0)  -  E 


(10.21b) 
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(10.25a) 
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C2  "  “  d-60)2(1+2/60)  (10.25b) 

+  V4C3  *  V»> 

C3  “  1  +  v„  "  fi0  (10.25c) 

C.  -  1  -  «  +  cl  +  { (1-6  )/6  >C,  (10. 25d) 

4  o  J  o  o  j 

Ce  -  2V  -  C,  (10. 25e) 

5  00  3 

Given  values  for  6  and  V  ,  condition  (10.24),  which  is  im- 
0  00 

posed  by  the  requirement  of  self-consistency  within  the  set  of  pertur¬ 
bation  equations  at  order  F^,  provides  the  functional  relationship 

I 

between  the  disturbance  frequency  ui  and  the  normalized  wavenumber  icA, 
which  is  always  a  real  quantity.  Since  (10.24)  is  cubic  in  to  with 
real  coefficients,  fixing  tcA  leads  to  three  roots  for  to.  Thus,  they 
are  either  (i)  all  real  or  (ii)  one  is  real  and  the  other  two  a  com¬ 
plex  conjugate  pair.  Disturbances  whose  wavelength  is  such  that  the 
solutions  to  (10.24)  are  all  real  are,  of  course,  stable,  the  ampli¬ 
tude  neither  increasing  nor  decreasing  with  time.  Those  for  which 
lm(w)<0  (>0)  are  unstable  (evanescent),  where  Im(<o)  is  the  imaginary 
part  of  the  complex  frequency.  Since  complex  roots  of  (10.24)  must 
occur  in  conjugate  pairs,  each  unstable  mode,  i.e.,  one  whose  ampli¬ 
tude  increases  exponentially  with  time,  is  necessarily  accompanied  by 
an  evanescent  mode,  i.e.,  one  whose  amplitude  decays  exponentially 
with  time.  The  real  parts  of  to  for  these  two  modes  are  equal.  Thus, 
implicit  in  the  exponential  front's  perturbation  dynamics  at  order 
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CHAPTER  11 


ANALYSIS  OF  THE  FIRST  ORDER  DISPERSION 
CHARACTERISTICS  OF  THE  EXPONENTIAL  FRONT 


11.1  THE  LOW  FREQUENCY  LIMIT 

The  dispersion  relation  at  order  F*,  (10.24),  may  be  written 

w3  +  C,kAa>2  +  C„(kX)2u  +  kA(l-6  )2 

1  2  °  (11.1) 

+  (kA)3  |  {(6o-l)/«o}C2  +  C4C5  |  -0  , 

where  the  C^  are  defined  in  (10.25).  The  point  (0,0)  in  (kA,w)  space 
lies  on  the  curve  w(kX)  defined  by  (11.1),  and  we  wish  to  investigate 
now  the  asymptotic  behavior  of  (11.1)  in  the  limit  w*0  and  kA-*0,  that 
is,  the  low  frequency  limit. 

Since  all  quantities  in  (11.1)  are  presumably  order  one  under 
the  scaling  introduced  earlier,  the  asymptotic  form  of  <*>(k*)  as  k^*0 
is  determined  by  assuming 

u  -  <3(kA)n  ,  (11.2) 

where  u  is  order  one.  Substituting  (11.2)  into  (11.1)  gives 


a)3  +  C1ai2(kX)1'n  +  C2<i(kX)2(1"n) 

+  (1-6  )2(kX)1_3n  +  (kX)3(1"n)G  -  0 
o 


(11.3) 
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where  the  constant  G  is 

G  -  {(6o-l)/5o>C^  +  C4C5  .  (11.4) 

If  all  quantities  are  to  remain  finite  for  kX-K),  (11.3)  implies 
that  l-3n*0,  or  n™l/3.  This  Condition  leads  to  the  following  asymp— 
form 


w3  ^  -(i-«o)2kX  ,  (11.5) 

in  which  cne  original  variables  have  been  reintroduced. 

The  three  roots  of  (11.5)  for  kX>0  are 


-{(1-6  )2kA}1/3 

0 

(11.6a) 

*s(l+i/3)  {  (1*6  )2kX}1/3 

0 

(11.6b) 

ij(l-i/3){(l-6  )2kX)1/3 

0 

(11.6c) 

where  0<kX«l.  corresponds  to  a  stable  mode,  whereas  d>2  and  u>3 
correspond  to  evanescent  and  unstables  modes,  respectively.  We  note 
that  the  low  frequency  limit  of  the  first  order  dispersion  relation 
is  independent  of  and  depends  parametrically  only  on  &Q . 

The  phase  velocities  associated  with  these  three  frequencies, 
which  are  computed  from  Re(u  )/kX,  are: 
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Coh  "  -U-«  )2/3(kX)r2/3 
ph1  o 

(11.7a) 

CDh  *  Coh  “  ‘sd-S  )2/3<kX)"2/3 
ph2  ph3  o 

(11.7b) 

The  stable  wave  thus  travels  in  the  -y  direction  with  twice  the  phase 

velocity  of  the  evansecent  and  unstable  waves,  both  of  which  propa¬ 
gate  in  the  +y  direction. 

The  group  velocities  for  the  three  modes  are 

dRe(<*>  )/d(kX)  as 
n 

computed  from 

C  -  -(1/3) (1-6  )2/3(kX)"2/3 

8i  0 

(11.8a) 

C  -  (1/6) (1-6  )2/3(kX)~2/3 
g7  o 

(11.8b) 

The  group  velocity  of  a  particular  mode  is  thus  1/3  of  its  phase 

velocity,  and  both  the  group  and  phase  velocities  are  proprtional  to 
-2/3 

(kX)  for  all  disturbances.  As  kX-MD,  C  and  C  ,  become  Infinite. 

g  ph 
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11.2  REGIONS  OF  INSTABILITY  IN  kX-SPACE 

The  nature  of  the  roots  of  (11.1)  are  easily  determined  by 
employing  the  following  general  result  for  cubic  equations  (CRC  Hand¬ 
book  of  Mathematical  Tables,  1962,  p.  387).  The  equation 

3  2 

y  +py  +qy+r«0 

has  (A)  one  real  and  two  conjugate  complex  roots  for  M>0 

(B)  three  real  roots,  of  which  at  least  two  are  equal, 

for  M-0 

(C)  three  real  and  unequal  roots  for  M<0  , 

2  3 

where  M*b  /4+a  /27,  and  where  a  and  b  are  defined  in  terms  of  the 

2  3 

coefficients  p,  q  and  r  as  follows:  a*q-p  /3  ,  b«2p  /27-pq/3+r  . 

The  obvious  equivalences  between  the  quantities  p,  q  and  r  above 
and  the  coefficients  in  (11.1)  are  not  reproduced  here. 

For  several  values  of  the  parameters  6q  and  V^,  the  discrim¬ 
inant  M  is  plotted  for  (11.1)  as  a  function  of  kX  in  Figures  11.1a 
through  11. le.  Ranges  of  kX  for  which  the  discriminant  is  zero  or 
negative  correspond  to  three  stable  modes  on  the  exponential  front, 
while  values  of  kX  giving  a  positive  discriminant  result  in  one  stable, 
one  euansecent  and  one  unstable  mode  as  discussed  earlier.  Since 
M(kX)  is  an  even  function  of  kX,  only  the  domain  kX>0  is  plotted  in 
Figure  11.1. 
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Plots  of  M  vs  <\  for  Vw  *  0.100.  Value  of 

6  In  Parentheses  on  Each  Plot. 

0 


Figure  11.1(a) 


KtJ 


Plots  of  M  vs  kA  for  V  *  0.500.  Values  of 


6  in  Parentheses  on  Each  Plot. 
0 


Figure  11.1(e) 
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11.3  FURTHER  ANALYSIS  OF  THE  DISPERSION  RELATION 

(11.1)  may  be  further  analyzed  by  decomposing  the  complex 
frequency  into  real  and  imaginary  parts,  i.e., 

a)  »  +  i«  (11.9) 

The  stability  characteristics  of  the  frontal  zone  are  determined 
by  u>  ,  while  the  phase  and  group  'velocities  are  set  by  co  and  its 
first  derivative.  Substituting  (11.9)  into  (11.1)  leads  to  the 
following  equation  for  wR  alone: 


-8 <oR  -  80^^  -  2(cJ+C2)(kA)2u>R 

-(kA)3C1C2  +  (kA)(l-6Q)2  +  (kA)3G  -  0 
where  G  appears  in  (11.4). 


(11.10) 


Replacing  kA  by  -kA  in  (11.10)  and  adding  the  resulting  ex¬ 
pression  to  (11.10)  gives 


-8{u3(kA)  +  u)3(-kX)  } 

+  8CLkA  {(D2(-kA)  -  u>R(kA) } 

-2(cJ  +C2)(k»2  (aiR(kA)  +  uR(-kA)}  -  0  . 


(11.11) 
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This  last  expression  is  satisfied  only  if  the  coefficient  of  each 
power  of  kA  vanishes  independently,  i.e., 

UR(kX)  +  uR(-kX)  =  0 

u>R(-kX)  -  <oR(kX)  *  0  (11.12) 

<*JR(kX)  +  toR(-kX)  -  0 

These  three  conditions,  in  turn,  require 

uR(-kA)  -  -<oR(kX)  .  ,  (11.13) 

from  which  we  see  that  the  real  part  of  the  complex  frequency  is 
an  odd  function  of  kX. 


Proceeding  in  a  similar  manner,  we  can  quickly  generate  the 

following  result  that  expresses  to.  (kA)  in  terms  of  to_(kX): 

l  K 


wJ(kA)  -  io*(-kA)  -  3(ioR(kA)  -  toR(-kX) } 

+  2C  kX{to  (kX)  +  to  (-kA)} 

X  K  R 


(11.14) 


But,  in  view  of  (11.13),  the  right-hand-side  of  (11.14)  vanishes  iden¬ 
tically,  leaving 


to^(-kX)  -  +toi(kX) 


(11.15) 
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This  result  at  first  appears  unusual,  since  it  seems  to  imply  that 
w^kA)  simultaneously  exhibits  both  even  and  odd  symmetry.  This  appar¬ 
ent  inconsistency  is  resolved,  however,  by  realizing  that  complex 
roots  of  (11.1)  must  occur  as  a  conjugate  pair,  which  is  the  fact  re¬ 
flected  in  (11.15). 

The  real  and  imaginary  parts  of  (11.1)  are 


0,3  o,  o,2  a>2  ^ 2 . 

UR  *  +  C1<UR  "  V 


+  cJX,  +  {(1-6  )/kA}2  +  G  -  0 

L  R  0 


o,2  o,2  o, 

3uR  -  “i  +  2Ci“r  +  C2  -  0  , 


(11.16) 


(11.17) 


respectively,  where 

u>  ■w/kl  *  +  iw^  .  (11.18) 

'If 

Eliminating  between  (11.16)  and  (11.17)  generates  a  single  equation 
for  a)i(kA),  which  is  analogous  to  (11.10)  for  the  real  part  of  u.  The 
derivation,  however,  is  considerably  more  involved  algebraically  than 
the  one  leading  to  (11.10),  and  the  details  are  not  reproduced  here. 


The  final  form  is 
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+  e^u>i  +  +  eQ  -  0  (11.19) 

where 

e6  -  192/81  (11.20a) 

e4  -  (32/27) A(2A(2  +  3/6q)  -  3>  (11.20b) 

e2  -  (1/27) (36£2 (A)  -  24f(A)A2  +  4A4}  (11.20c) 

eQ  -  4A4(2A/27  -  1)  -  (4/27)f3(A) 

+(4/27) f (A)A2  (18A  -  A2  -  2f(A)} 

(11. 20d) 

+  4A{ A  -  A2/27  -  (l/3)*f(A)}(A/kX)2 

-  (A/kX)4 

f  (A)  -  A(1  -  A  (1+2/ <5q)  }  (11.20e) 

A  -  1  -  6 

0 


(11. 20f ) 
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This  is  a  remarkable  result,  in  that  the  coefficients  e^  through  e^  do 
not  depend  on  the  shear  velocity  V^.  Thus,  w^is  independent  of  V-t  or, 
equivalently,  the  frontal  zone's  stability  characteristics  are  indepen¬ 
dent  of  V  .  The  real  part  of  u,  on  the  other  hand,  does,  in  general, 

09 

depend  on  V  ,  so  that  the  phase  and  group  velocities  of  first  order 
frontal  disturbances  do  exhibit  a  parametric  dependence  on  V  .  As  we 
have  seen  in  Section  11.1,  only  in  the  low  frequency  limit  is  the  real 
part  of  <i>  independent  of  V^.  We  also  note  that  only  e^  contains  any  kA 
dependence;  each  of  the  other  coefficients  is  a  function  of  $Qalone. 
Thus,6Q  appears  to  be  the  dominant  parameter  in  determining  the  frontal 
zone's  stability. 


11.4  NUMERICAL  RESULTS 

(10.24)  has  been  solved  numerically  for  a  range  of  the  parameters 
&  and  V  .  For  fixed  kA  ,  the  three  roots  of  (10.24)  were  computed  it- 

0  oo 

eratively  using  the  IMSL^  algoritm  ZRPOLY,  which  provided  exceptional 

accuracy.  Substitution  of  the  calculated  roots  back  into  (10.24)  gave 

13 

equality  to  within  typically  one  part  in  10  . 

Table  11.1  compares  the  frequencies  calculated  from  the  exact  dis¬ 
persion  relation,  (10.24),  with  those  calculated  from  its  low  frequency 
asymptotic  form  in  (11.6)  for  the  case  6  *0.35  and  V  -0.20.  Since  the 

0  oo 

complex  frequencies  (roots)  must  occur  as  a  conjugate  pair,  only  one 
complex  root  is  tabulated  in  Table  11.1,  along  with  the  purely  real  root. 
Results  obtained  with  the  low  frequency  limit  are  in  parentheses.  It  is 
apparent  from  the  table  that  the  low  frequency  limit  is  a  good  approxi- 

I 

mat ion  only  for  large  disturbance  wavelengths,  or,  conversely,  small  k\. 
Agreeement  to  within  20Z  or  so  is  obtained  out  to  kA*0.10.  As  kA  increases, 
the  low  frequency  limit  becomes  a  progressively  poorer  approximation. 


MMMNMi 
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Parameter  Values:  6  *0.35, 

0 

V  =0.20 

00 

Complex  Root 

kX 

Real  Part 

Imag  Part 

Real  Root 

0.00 

0.000(0.000)* 

0.000(0.000) 

0.000(  0.000) 

0.01 

0.077(0.081) 

0.140(0.140) 

-0.166 (-0.162) 

0.02 

0.094(0.102) 

0.175(0.176) 

-0.213 (-0.204) 

0.03 

0.105(0.117) 

0.199(0.202) 

-0.248(-0.233) 

0.04 

0.114(0.128) 

0.218(0.222) 

-0. 278(-0. 257) 

0.05 

0.120(0.138) 

0.232(0.239) 

-0.303(-0.276) 

0.06 

0.126(0.147) 

0.245(0.254) 

-0. 327(-0. 294) 

0.07 

0.131(0.155) 

0.255(0.268) 

-0.349 (-0.309) 

0.08 

0.135(0.162) 

0.264(0.280) 

-0. 369(-0. 323) 

0.09 

0.138(0.168) 

0.271(0.291) 

-0.389 (-0.336) 

0.10 

0.141(0.174) 

0.278(0.302) 

-0. 408(-0. 348) 

0.15 

0.153(0.199) 

0.295(0.345) 

-0.493(-0.399) 

0.20 

0.160(0.219) 

0.291(0.380) 

-0.57 1 (-0.439) 

0.25 

0.165(0.236) 

0.269(0.409) 

-0.643(-0.473) 

0.30 

0.168(0.251) 

0.220(0.435) 

-0. 712(— 0. 502) 

0.35 

0.170(0.264) 

0.119(0.458) 

-0.777 (-0.529) 

0.40 

-0.006(0.276) 

0.000(0.479) 

-0.84K-0.553) 

0.45 

-0.126(0.288) 

0.000(0.498) 

-0.903(-0.575) 

0.50 

-0.226(0.298) 

0.000(0.516) 

-0.962(-0.596) 

*Results  in  parentheses 
limit 

computed  from 

the  low  frequency 

Table  11.1  Comparison  of  Low  Frequency  Limit 
and  Exact  Solution  to  the  Dispersion 
Relation 


The  exact  frequencies  show  a  band  of  instability  for  0<kX<0.40 
as  a  result  of  the  nonzero  imaginary  part  of  the  complex  root.  For 
kX20.40,  all  three  modes  are  stable,  since  the  imaginary  part  of  the 
complex  frequency  is  zero.  The  real  part  of  the  complex  frequency  in¬ 
creases  from  zero  at  kX=0  to  a  maximum  positive  value  near  <X*0.35, 
and  then  quickly  decreases  to  negative  values  for  kX>0.40.  The  behav¬ 
ior  of  the  real  and  imaginary  parts  of  the  complex  frequency  is  thus 
clearly  not  monotonic  in  kX.  The  low  frequency  limit,  however,  is 
monotonlc  in  kX,  which,  in  part,  accounts  for  its  poor  performance  as 
an  approximation  when  kX  is  large. 

Figures  11.2  through  11.4  show  plots  of  the  exact  dispersion  re¬ 
lation  (real  and  imaginary  parts  of  u  vs  kX)  at  order  F^  for  the  expo- 

b 

nentlal  front  computed  from  (10.24).  They  illustrate  the  qualitative 
behavior  of  the  dispersion  characteristics  as  the  parameters  6  and  V 

0  oo 

are  varied.  Figures  11.5  and  11.6  provide  typical  phase  and  group  vel¬ 
ocity  plots,  computed  from  Cp^»Re(«) /<X  and  Cgr=dRe(cj) /d(icX)  ,  respectively, 

for  various  values  of  6  and  V  . 

o  op 

Figures  11.2a  through  11. 2e  plot  the  imaginary  part  of  each  of 

the  three  complex  roots  (frequencies)  as  a  function  of  icX  for  <5  =0.5  to 

o 

0.1  in  increments  of  0.1,  respectively.  Since  w.  is  independent  of  V  , 

i 

the  results  in  this  set  of  curves  apply  for  any  value  of  V  .  The  real 

00 

part  of  the  complex  frequency  does  depend  on  V  ,  however,  and  its  value 

00 

therefore  appears  in  the  plot  annotation.  Each  of  the  three  roots  of 
the  dispersion  relation  is  marked  by  a  different  symbol  on  the  plots  in 
Figures  11.2  through  11.6  (circle,  diamond  and  triangle). 

(1)  International  Mathematical  and  Statistical  Libraries,  Houston,  TX  77036 
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In  Figure  11.2,  the  kA  axis  is  always  overplotted  because  the 
cubic  dispersion  relation  must  always  contain  at  least  one  purely  real 
root  with  a  corresponding  zero  imaginary  part.  Since  the  complex¬ 
valued  dispersion  relation  roots  must  occur  always  as  a  conjugate  pair, 
unstable  and  evanescent  modes,  corresponding  to  Im(co  )<0  and  lm(u)>0, 
respectively,  must  also  occur  as  a  pair.  Thus,  every  unstable  mode 
is  accompanied  by  an  evanescent  mode,  and  this  fact  is  reflected  in 
the  mirror  image  symmetry  about  the  ucA-axis  in  Figure  11.2.  In  add¬ 
ition,  Im(u)  is  also  symmetrical  about  the  m-axis,  as  discussed  in 
Section  11.3.  It  is  therefore  sufficient  to  consider  only  icAkO  in 
discussing  stability. 

Figure  11.2  reveals  two  symmetric  bands  of  instability,  that  is, 
two  regions  on  nonzero  Im(oj) .  The  one  referred  to  as  Band  1  is  cen¬ 
tered  about  the  icA  origin,  while  the  other.  Band  2,  is  bounded  away 
from  kA*0.  Band  1  is  characterized  by  the  maximum  normalized  wave- 

number  for  modes  in  this  band,  which  is  labeled  kA  in  Figure  11.2a, 

max 

and  by  the  (kA,u)  pair  that  gives  the  fastest  growing  unstable  mode  in 

this  band.  This  point  on  the  dispersion  curve  is  marked  by  the  largest 

value  of  |lm(u)|  ,  as  shown  in  Figure  11.2a.  Examination  of  Figure  11.2 

shows  that  |lm(tu)|  is  mono tonically  increasing  with  increasing  kA  for 

modes  in  Band  2,  so  that  modes  in  Band  2  exhibit  progressively  larger 
» 

imaginary  parts  of  the  disturbance  frequency  with  increasing  wavenumber. 
Therefore,  no  fastest  growing  mode  exists  in  Band  2.  The  only  parameter 
characterizing  Band  2  is  the  minimum  wavenumber  for  modes  contained  in 
it,  labeled  <Affl^n  in  Figure  11.2a. 

Each  unstable  mode  possesses  a  time  dependent  amplitude  which 
doubles  in  the  interval  t^-ln(2/F^)/ | Im(ui) | .  The  point  (0,0)  is  always 
contained  in  the  dispersion  curve.  Thus,  for  any  finite  values  of  the 
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parameters  6  and  F.  ,  as  the  wavelength  becomes  arbitrarily  large,  so 

O  D 

does  the  doubling  time.  The  doubling  time  also  becomes  infinite  for 
unstable  modes  with  nonzero  frequencies  as  F,-*0.  At  order  F*  for  the 

D  D 

exponential  front,  therefore,  unstable  modes  evolve  continuously  from 
stable  ones  with  increasing  cross-stream  flow.  For  zero  cross-stream 
velocity,  all  modes  are  stable,  but  the  doubling  time  for  a  given  un¬ 
stable  mode  decreases  from  Infinity  with  increasing  cross-stream  flow. 

The  existence  of  nonzero  cross-stream  velocity  is  therefore  of  paramount 
importance  in  determining  the  stability  characteristics  of  the  exponen¬ 
tial  front. 

The  limiting  normalized  wavenumbers  for  Bands  1  and  2,  kA  and 
°  max 

kA  .  as  discussed  above,  are  tabulated  as  functions  of  6  in  Table  11.2. 
min  o 

The  maximum  wavenumber  for  modes  in  Band  1  increases  with  increasing  6  , 

o 

as  does  the  minimum  wavenumber  for  Band  2  modes.  kA  is  nearly  equal 

max 

to  6  ,  and  this  observation  serves  as  a  useful  rule  of  thumb  for  setting 
0 

the  limits  of  Band  1  modes  in  the  range  0.1<6  <0.5.  Recall  that  the 

“  o 

results  in  Table  11.2  are  independent  of  V  . 

OO 

The  characteristics  of  the  fastest  growing  mode  in  Band  1  are  tab¬ 
ulated  as  functions  of  6  in  Table  11.3  for  V  *0.2.  Knowing  V  is  import- 

Q  00  00 

ant  in  this  case,  since  Re(u),  C  .  and  C  depend  upon  it.  <A,  Im(bi) 

pn  gr 

and  t^,  however,  do  not  depend  on  V^,  so  that  the  tabulated  results  for 

these  quantities  apply  generally  for  6  in  the  range  considered.  Obsetf- 

o 

vatlonal  data  tend  to  show  that  F^,  or,  equivalently,  the  cross-stream 
flow,  is  quite  small  in  large  scale  upper  ocean  density  fronts.  It  is 
unlikely,  therefore,  that  values  of  F,  much  larger  than  0.2  would  be 

D 

associated  with  the  frontal  structures  considered  in  this  paper.  Doubling 
times  in  Table  11.3  are  therefore  computed  for  F^«0.05  and  F^*0.20,  which 
are  taken  as  reasonable  representative  limits  for  this  parameter. 
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Table  11.3  shows  that  the  wavenumber  of  the  fastest  growing 

mode  is  nearly  equal  to  6  /2,  which  is  a  useful  rule  of  thumb  for 

0 

estimating  its  value.  As  6  increases  by  a  factor  of  5,  kX  increases 

o 

by  a  factor  of  about  4.  The  real  and  imaginary  parts  of  the  disturb¬ 
ance  frequency,  however,  are  far  less  sensitive  to  changes  in  6  .  The 

o 

real  part  of  a>  increases  very  slightly  as  6  increases  from  0.1  to  0.2, 

o 

and  thereafter  decreases  by  a  factor  of  less  than  1/3  as  6  continues 

o 

to  0.5.  The  imaginary  part  of  u>  is  even  less  sensitive  to  varying  6  , 

o 

and  remains  between  -0.263  and  -0.297  as  6  ranges  over  0.1-0. 5.  The 

o 

phase  and  group  velocities  both  decrease  with  increasing  6  and  show 

o 

considerable  variation.  The  shortest  doubling  time,  normalized  by  the 

Inertial  period,  occurs  for  6  *0.3.  It  is  not  particularly  sensitive 

o 

to  changes  in  F,  as  a  result  of  the  logarithmic  dependence  on  F.  .  As 

D  D 

F^  increases  by  a  factor  of  4,  the  doubling  time  decreases  by  somewhat 
less  than  50Z.  Note  that  all  entries  in  Table  11.3  are  dimensionless. 

In  order  to  compare  these  results  with  observational  data, 

Tables  11.4  through  11.6  provide  the  most  important  characteristics 
of  the  fastest  growing  mode  as  dimensional  quantities  for  the  data  in 
Table  11.3.  The  disturbance  period  in  days  is  computed  from  g* 
2ir/Re(u) *f *8.64x10^,  where  f  is  the  Coriolis  parameter.  For  convenient 
reference,  the  Coriolis  parameter  is  tabulated  in  Table  11.4.  The  wave 
length  is  computed  from  2irX/(«X).  Table  11.5  shows  the  disturbance  per 
iod  in  days  at  various  latitudes  as  a  function  of  Re(u>) ,  whose  values 
are  given  in  Table  11.3  for  the  fastest  growing  mode.  We  see  that  the 
disturbance  period  decreases  with  Increasing  latitude,  and  increases 
with  decreasing  Re(u>)<  For  latitudes  between  20°  and  50°,  the  fastest 
growing  mode  with  V  -0.2  exhibits  periods  between  3.63  days  and  12.50 

Oft 
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days  for  0.1<6  <0.5.  For  6  *0.3,  corresponding  to  Re(ui)>0.166,  the 
0  0 

periods  at  20° ,  30°  and  40  latitude  are  8.81,  6.02  and  4.69  days,  re¬ 
spectively.  These  values  are  in  fair  agreement  with  the  dominant  periods 
observed  using  satellite  imagery  of  the  Gulf  Stream  by  Maul,  et.  al  (1978). 
Their  spectral  analysis  of  Gulf  Stream  Meanders  revealed  dominant  per¬ 
iods  of  10,  5.8  and  4.8  days,  respectively,  at  latitudes  of  approxi¬ 
mately  27°N,  33°N  and  38°N.  The  variation  of  meander  period  with  lati¬ 
tude  seems  to  be  strongly  influenced  by  Coriolis  effects.  Using  the  fas¬ 
test  growing  mode  wavenumbers  in  Table  11.3,  the  corresponding  dimension¬ 
al  wavelengths  appear  in  Table  11.6  for  two  values  of  the  baroclinic  Rossby 
radius,  X.  Although  these  wavelengths  are  seen  to  be  quite  long,  they 
are  not  necessarily  inconsistent  with  observed  frontal  meander  patterns, 
which  result  from  the  superposition  of  a  continuum  of  wavelengths  ranging 
from  short  to  very  long. 

Table  11.7  displays  the  characteristics  of  the  fastest  growing 
mode  in  Band  1  for  V  >0.6  In  order  to  examine  the  sensitivity  to  changes 

OB 

in  V  .  Only  those  parameters  that  are  dependent  upon  Vn  are  Included  in 
this  table.  While  the  data  in  Table  11.3  show  that  the  real  part  of  u> 
does  not  change  substantially  with  Increasing  6  for  V  >0.2,  the  data  in 

o  OB 

Table  11.7  reveal  that  Re(iu)  is  far  more  sensitive  to  changes  in  5  for 

0 

V  >0.6.  As  5  changes  by  a  factor  of  5,  from  0.1  to  0.5,  Re(u>)  decreases 

OB  O 

by  a  factor  of  7,  from  0.155  to  0.023.  At  V  >0.2,  Re(o>)  was  monotonically  de- 

OB 

creasing  only  for  increasing  $  greater  than  0.2.  For  V  >0.6,  however, 

Q  OB 

Re(u)  decreases  with  increasing  6  between  0.1  and  0.5.  As  V  Increases 

O  OB 

from  0.2  to  0.6,  the  phase  velocity  at  6  -0.1  decreases  only  slightly, 

o 

while  at  6q-0.5  it  decreases  by  nearly  a  factor  of  5.  As  in  Table  11.3, 
is  monotonically  decreasing  with  Increasing  6q.  For  fixed  <Sq  the 
change  in  group  velocity  as  V  ncreases  from  0.2  to  0.6  is  larger  than 
the  corresponding  change  in  C^.  is  again  monotonically  decreasing 

with  Increasing  $oat  V  -0.6.  The  period  in  days  at  latitudes  of  20°,  30° 

40°  and  50°  corresponding  to  the  tabulated  values  of  Re(u)  constitute  the 
last  four  columns  in  Table  11.7.  Since  Re(u)  is  considerably  more  sensitive 
to  variations  in  6  at  V  -0.6  than  it  is  at  V  -0.2,  the  periods  in  Table 

Q  OB  08 

11.7  exhibit  far  more  variation  with  6  at  a  given  latitude  than  those  in 

o 

Table  11.5.  The  maximum  period  is  63.57  days  at  6  -0.5  and  50  latitude. 

0 

For  intermediate  values  of  6  ,  where  Re(u)  at  V  -0.6  is  not  drastically 

Q  OB 


different  from  its  corresponding  value  at  V  =0.2,  the  periods  in  Table 

00 

11.7  and  11.5  are  not  greatly  different.  Except  at  6^=0. 4  and  6q=0.5, 
even  at  V  =0.6,  the  disturbance  periods  are  thus  reasonably  consistent 
with  the  observational  data  of  Maul,  et  al, ,  (1978).  As  a  general  con¬ 
clusion,  with  increasing  ,  the  disturbance  period  is  progressively 
less  sensitive  to  the  change  in  V  as  6  decreases. 

oo  o 

Using  the  same  format  as  Table  11.7,  data  for  V  between  zero 

00 

and  2.00  are  summarized  in  Table  11.8. 
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TABLE  11.2  ! 
Maximum  and  Minimum  Normalized  . 
Wavenumber  for  Bands  1  and  2 


! 


Band  1 

Band  2 

6 

_ 0 

kX 

max 

kXmin 

ai 

0.1 

0.154 

*  0.77 

0.2 

0.248 

*  0.88 

-  » 

0.3 

0.321 

*  1.77 

0.4 

0.406 

>  2.00 

?J 

0.5 

0.490 

>  2.00 

a* 

9 


r 


-  - 


jrsw<7s»^r  — 
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TABLE  11.3 

Characteristics  of  the  Fastest 
Growing  Mode,  Band  1,  V*  -  0.2 


io_ 

K  A 

Re  (w) 

Im(w) 

C 

gr 

Fb=0.05 

Ld 

Fb=0.20 

0.1 

0.063 

0.180 

-0.263 

2.85 

1.17 

14.03 

8.76 

0.2 

0.111 

0.182 

-0.290 

1.65 

0.54 

12.72 

7.94 

0.3 

0.147 

0.166 

-0.297 

1.13 

0.26 

12.42 

7.75 

0.4 

0.189 

0.145 

-0.293 

0.77 

0.08 

12.59 

7.86 

0.5 

0.237 

0.117 

-0.281 

0.49 

-0.05 

13.13 

8.19 
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TABLE  11.4 

Values  of  the  Coriolis 
Parameter  as  a  Function  of  Latitude 


LAT  (Deg) 

f  (Sec 

-1) 

0 

0 

10 

2.526 

x  10-5 

20 

4.974 

x  10‘5 

30 

7.272 

x  10"5 

40 

9.349 

x  10~5 

50 

1.114 

x  10-4 

60 

1.260 

<3- 

1 

O 

»H 

X 

70 

1.376 

x  10"4 

80 

1.432 

x  10-4 

90 

1.454 

x  10~4 

1 

I 

1 

Jk 
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TABLE  11.5 

Period  at  Various  Latitudes 


Period  In  Days  at  Latitude  In  Degrees 


Re(w) 

20° 

30° 

40° 

50° 

0.180 

8.12 

5.56 

4.32 

3.63 

0.182 

8.03 

5.49 

4.27 

3.59 

0.166 

8.81 

6.02 

4.69 

3.93 

0.145 

10.08 

6.90 

5.36 

4.50 

0.117 

12.50 

8.55 

6.65 

5.58 
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kX 

0.063 

0.111 

0.147 

0.189 

0.237 


TABLE  11.6 

Wavelength  for  "Various  Values  of  kA 
and  Rossby  Radius 


X  «  50  km 
4987 
2830 
2137 
1662 
1326 


Wavelength  (km) 

X  =  100  km 
9973 
5661 
4274 
3324 
2651 
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CHAPTER  12 


concurs  iom 

This  work  has  examined  Che  stability  characteristics  of  a  shallow 
upper  ocean  density  front  with  son-zero  cross-stream  flow.  The  frontal 
zone  is  a  wedge  shaped  body  of  fluid  with  constant  density  that  floats 
on  an  ambient  ocean  of  slightly  greater  density,  the  dynamics  of  which 
are  unaffected  by  the  presence  of  the  frontal  zone.  The  model  is  hydro- 
dynamic  in  nature,  since  the  density  field  is  specified  a  priori.  Earth 
rotation  effects  are  incorporated  under  the  f-plane  approximation.  Time- 
varying  disturbances  of  the  frontal  zone  are  assumed  to  have  been  initiated 
by  some  driving  agent,  surface  wind  stress,  for  example,  which  is  no  longer 
active.  Thus,  only  free  oscillations  in  the  system  are  considered.  The 
ambient  ocean,  whose  flow  is  geostrophic,  and  the  frontal  zone  are  assumed 
to  be  in  isostatic  balance. 

The  frontal  zone  is  divided  into  two  regions,  an  inner  region  in 
which  dissipation  occurs  and  an  outer  region  in  which  the  dynamics  are 
strictly  inviscid.  Mass  entrainment  in  the  inner  region  produces  a  non¬ 
zero  mass  flux  at  the  boundary  between  it  and  the  outer  region,  so  that 
the  inviscid  outer  region  dynamics  respond  indirectly  to  dissipative  pro¬ 
cesses.  This  situation  corresponds  to  large  values  of  the  rotation  param- 
eter  introduced  in  Garvine's  (1979a)  frontal  model,  and  is  typical  of 
large  scale  oceanic  fronts  in  which  Coriolis  effects  are  dominant.  Time- 
and  space-harmonic  disturbances  are  assumed  to  exist  in  the  frontal  zone, 
and  are  treated  as  small  amplitude  time-varying  perturbations  to  a  time- 
invariant  equilibrium  flow.  Both  planar  and  exponential  frontal  zone 
geometries  are  considered,  with  respective  equilibrium  flows  that  are 
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geostrophic  and  aeostrophic. 

Duxbury's  (1963)  model  of  wave  disturbances  on  a  planar  front, 
which  was  formulated  with  zero  cross-stream  flow,  is  first  extended  to 
include  a  nonzero  cross-stream  velocity  component.  The  initial  analysis 
in  this  paper  is  performed  in  dimensional  variables,  as  was  Duxbury's. 

The  equilibrium  flow  for  the  planar  front  is  spatially  invariant,  and 
the  basic  state  is  consequently  geostrophic.  Constant  long-  and  cross¬ 
stream  velocity  components  are  driven  by  the  sloping  sea  surface  and 
frontal  zone  discontinuity  surface,  both  of  which  are  planes,  which  mot¬ 
ivates  the  planar  front  nomenclature. 

By  eliminating  the  velocity  perturbations  as  dependent  variables, 
a  general  seventh  order  ordinary  differential  equation  (ODE)  is  developed 
for  the  planar  front  depth  perturbation  amplitude  (DPA) .  An  asymptotic 
normal  form  solution  is  presented  which  is  valid  far  from  the  surface 
front,  and  an  integral  representation  of  the  solution  is  developed  which 

is  valid  in  any  region.  In  addition,  the  connection  between  these  solutions 

% 

is  established,  and  it  provides  an  important  insight  into  the  proper  decay 
scale  for  perturbations  to  the  planar  front  system.  With  these  results, 
the  problem  of  the  planar  front  with  cross-stream  flow  has  been  solved  in 
principle  by  reducing  its  solution  to  quadratures.  For  the  case  of  zero 
cross-stream  flow,  the  general  seventh  order  ODE  for  the  DPA  reduces  to 
a  second  order  equation,  whose  integral  representation  solution  is  evalu¬ 
ated  explicitly  to  recover  all  of  Duxbury's  earlier  results. 

However,  the  analysis  described  above  raises  fundamental  questions 
concerning  the  planar  front's  dynamics.  In  dimensional  variables,  the 
system  with  zero  cross-stream  flow  is  second  order,  while  the  system  with 
nonzero  cross-stream  flow,  however  small,  is  seventh  order.  This  result 


m 

raises  the  question  of  what  the  fundamental  order  of  the  system  is  and 
why  it  changes  so  drastically  with  vanishing  cross-stream  flow.  The  an¬ 
swers  lie  in  the  proper  scaling  of  the  problem,  in  which  spatial  variables 
are  scaled  to  the  invlscld  region  characteristic  length  (the  Rossby  radius) 
and  an  advective  characteristic  length  which,  in  fact,  is  so  small  as  to 
be  imbedded  in  the  narrow  dissipative  zone.  A  small  parameter  expansion 
is  then  introduced  which  involves  the  cross-stream  internal  Froude  number 
(a  measure  of  the  strength  of  the  cross-stream  flow  normalized  to  the  in¬ 
ternal  wave  phase  speed) .  At  zeroth  and  first  orders  in  the  Froude  num¬ 
ber,  the  system  normalised  by  the  invlscld  scale  length  is  second  order, 
while  the  corresponding  zeroth  crder  system  normalized  by  the  advective 
scale  length  is  fifth  order.  Since  the  advective  zone  is  imbedded  in  the 
dissipative  region  located  near  the  surfaoe  front,  however,  the  assumption 
of  invlscld  dynamics  is  no  longer  valid  and  the  model  with  the  fifth  order 
system  is  inapplicable.  On  the  Invlscld  scale,  which  is  the  only  region 
of  direct  applicability  of  this  model,  the  planar  front  dynamics  are  seen 
to  be  fundamentally  second  order. 

The  inviscid  scale  planar  front  with  nonzero  cross-stream  flow 
thus  reduces  essentially  to  Duxbury's  original  second  order  system,  which 
exhibited  no  instabilities  at  all.  This  result,  which  is  available  only 
with  an  appropriate  scaling  of  the  planar  front  problem,  motivates  con¬ 
sideration  of  the  exponential  front  with  cross-stream  flow.  The  entire 
analysis  in  this  latter  case  is  performed  in  scaled  variables.  For  com¬ 
pleteness,  the  planar  front  is  reexamined  using  essentially  the  same 
scaling  scheme  in  a  parallel  development  with  the  exponential  front.  The 
dispersion  relation  for  each  geometry  is  then  developed  at  zeroth  and  first 
orders  in  the  cross-stream  IntenuQ.  Froude  number,  which  is  used  to  intro¬ 
duce  a  small  parameter  expansion  for  each  dependent  variable  in  the  problem. 
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The  lowest  order  dispersion  characteristics  of  the  planar  front 
are  particularly  simple.  The  disturbance  frequency  and  wavenumber  are 
related  by  a  linear  equation  which  involves  only  a  single  parameter,  V  , 
the  shear  velocity  between  the  ambient  fluid  and  the  frontal  zone.  Wave 
disturbances  to  lowest  order  in  the  small  parameter  F,  ,  the  cross-stream 

D 

internal  Froude  number,  are  therefore  non-disperslve  with  equal  phase 
and  group  velocities.  In  addition,  real  values  of  the  wavenumber  corres¬ 
pond  to  purely  real  disturbance  frequencies.  The  planar  front  is  there¬ 
fore  dynamically  stable  under  all  conditions  at  lowest  order. 

At  first  order  in  F^,  however,  the  planar  front  dispersion  results 
are  somewhat  surprising.  The  system  of  differential  equations  satisfied 
by  the  first  order  perturbation  amplitudes  is  internally  self-consistent 
and  meets  the  boundary  conditions  without  imposing  any  characteristic 
relation  between  the  disturbance  frequency  and  wavenumber.  .Therefore,  no 
first  order  dispersion  relation  exists  for  the  planar  front. 

For  the  exponential  front,  whose  basic  state  is  formulated  in  Chap¬ 
ter  5,  the  equilibrium  flow  is  not  spatially  invariant,  and  it  is  con¬ 
sequently  ageostrophic.  Both  exact  and  approximate  solutions  for  the  equil¬ 
ibrium  interfacial  depth  are  developed,  with  the  approximate  solution  again 
involving  a  small  parameter  expansion  in  terms  of  F.  .  At  lowest  order  in 

D 

F^,  the  basic  state  interfacial  depth  profile  is  exponential,  which  leads 
to  the  exponential  front  nomenclature. 

A  single  second  order  ODE  for  the  exponential  front  DPA  at  order  F^ 
is  developed  by  eliminating  the  velocity  field  perturbation  amplitudes  as 
dependent  variables.  Analysis  of  this  equation  reveals  that  the  exponen¬ 
tial  front  is  sub-lnertlally  dominated.  Disturbances  leading  to  non-diver- 
gent  solutions  for  the  depth  perturbation  amplitude  far  from  the  surface 
front  must,  in  general,  have  periods  greater  than  the  inertial  period. 
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This  existence  condition,  moreover,  depends  only  on  a  single  parameter, 
the  shear  velocity  between  the  frontal  zone  and  the  ambient  fluid,  as 
did  the  lowest  order  dispersion  relation  for  the  planar  front. 

For  the  exponential  front  at  zeroth  order  in  F^,  the  dispersion 
relation  is  linear  in  u  and  kX,  as  it  was  for  the  planar  front  at  order 
F®.  It  depends  on  two  parameters,  however,  the  shear  velocity  V  and  the 

D  * 

interfacial  depth  SQ  at  the  coordinate  system  origin.  For  real  values 
of  wavenumber,  the  frequency  is  also  real,  and  the  system  is  dynamically 
stable.  Wave  disturbances  are  non-dispersive  with  equal  phase  and  group 
velocities.  With  realistic  values  for  the  shear  velocity  V-f  the  phase 
and  group  velocities  will  be  positive  along  the  +y  direction  in  a  right- 
handed  coordinate  system  oriented  such  that  the  -He  axis  is  in  the  cross¬ 
stream  direction  into  the  fromtal  zone  from  the  ambient  fluid.  This  re¬ 
sult  is  consistent  with  observed  oeaoder  patterns  on  the  Gulf  Stream, 
in  which  frontal  zone  waves  propagate  northward  and  eastward  along  the 
Stream. 

At  first  order  in  F^,  the  exponential  front  possesses  a  cubic  dis¬ 
persion  relation,  which,  for  certain  values  of  the  purely  real  wavenumber, 
can  lead  to  complex  frequencies  with  nonzero  imaginary  parts,  and,  there¬ 
fore,  to  instabilities.  Both  Vm  and  enter  as  parameters  in  determining 
the  real  part  of  the  complex  frequency,  which,  in  turn,  governs  the  phas« 
and  group  velocities.  On  the  othor  hand,  the  stability  characteristics 
of  the  first  order  system,  which  axo  determined  by  the  imaginary  part  of 
the  complex  frequency,  are  independont  of  the  shear  velocity  V  .  They 
depend  instead  only  on  the  interfacial  depth  at  the  inner/outer  zone  boun¬ 
dary,  dQ. 

In  the  low  frequency  limit  am  u  and  kX  go  to  zero,  an  unstable  and 
an  evanescent  mode  always  exist  am  a  pair  on  the  exponential  front  with 
equal  phase  velocities.  The  corresponding  group  velocities  are  also  equal. 
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and  equal  one-third  of  the  phase  velocity  for  each  mode.  The  phase  velocity 
of  the  stable  wave  in  the  low  frequency  limit  is  twice  that  of  the  unstable 
and  evanescent  waves  and  of  opposite  sign.  Its  group  velocity  is  also 
one-third  of  its  phase  velocity. 

Instabilities  of  the  exponential  frontal  zone  vanish  as  F.-K),  or, 

D 

equivalently,  as  the  equilibrium  state  cross-stream  velocity  vanishes. 

The  growth  time  for  unstable  waves  is  proportional  to  ln(l/F,  ),  and  thus 

o 

approaches  infinity  as  F^-K),  reflecting  a  continuous  transition  from  an 
Inherently  stable  system  at  order  F^  to  a  potentially  unstable  system  at 

D 

order  F^.  Computational  results  for  typical  values  of  between  0.1  and 
0.5  with  F^  between  0.05  and  0.2  show  doubling  times  on  the  order  of  1  to 
2  weeks,  which  seem  reasonable,  although  no  direct  observational  data  are 
available  for  comparison. 

Calculation  of  the  disturbance  periods  over  a  wide  range  of  and 
V  values  shows  that  they  are  also  typically  1  to  2  weeks,  which  is  in  gen¬ 
eral  agreeement  with  Gulf  Stream  observational  data.  In  addition,  the 
wave  periods  exhibit  a  strong  dependence  upon  the  Coriolis  parameter,  or, 
equivalently,  geographic  latitude,  which  is  also  evident  in  observations 
of  the  Gulf  Stream.  In  particular,  for  V^-0.2  and  6q=»0 . 3 ,  both  reason¬ 
able  values,  the  disturbance  periods  for  the  fastest  growing  mode  are  8.81, 
6.02  and  4.69  days  at  latitudes  of*  20°,  30°  and  40°,  respectively.  At  lati¬ 
tudes  of  approximately  27°N,  33°N  and  38°N,  Maul,  ££  al. ,  (1978)  observed 
dominant  periods  of  10,  5.8  and  4.8  days,  respectively,  by  performing  a 
spectral  analysis  of  available  satellite  infra-red  imagery. 

The  work  presented  here  shows  that  upper  ocean  density  fronts  can 
exhibit  barodinlc  instabilities,  while  all  previous  theoretical  efforts 
have  concluded  otherwise.  The  system's  stability  characteristics  are  seen 
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to  depend  critically  upon  the  existence  of  some  cross-stream  flow, 
which,  in  turn,  is  supported  by  dissipative  processes  in  a  narrow  zone 
near  the  surface  front.  The  quasi-steady  dynamics  of  this  region  have 
been  modelled  elsewhere  (Garvln«t  1979a,  1974).  A  natural  extension  of 
this  work  would  be  to  consider  earth  rotation  effects  under  a  beta-plane 
approximation  and  to  Include  coupling  between  the  frontal  zone  and  am¬ 
bient  ocean  dynamics.  Another  important  extension  would  be  the  develop¬ 
ment  of  an  energetics  analysis,  which  should  provide  insight  into  the 
physical  processes  at  work  in  the  frontal  zone  and  perhaps  display  the 
physical  origin  of  the  instability. 
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APPENDIX  1 


COEFFICIENTS  OF  £HZ  GENERAL  EQUATION 
FOR  Z(x)  FOR  THE  PLANAR  FRONT 

The  operators  through  L^  and  the  operator  polynomial  A(p) 

appearing  In  (3.11)  are  defined  In  (3.10).  These  operators  may  be 

expanded  and  substituted  Into  (3.11)  to  yield  the  final  form  of  the 

depth  perturbation  amplitude  ODE  In  (3.13),  with  the  constants  \>n 

and  Y  given  by: 
a 

2  3  2  2 

v0  -  la^'k w  (r-w  ) 

yQ  •  (f2-ta2){a2g,kw3  -f  i{g'k2<»>3(Df  (0)+a2y) 

-  a^g'k<i)2f  +  <o3(f2-u2)}> 

+  2a1g,k2a,2G(f2-3«2) 

-  a^g'k2(j2u(3f2-5«2)  (A.lc) 

Yl  -  g,k2u2u(3f2-3*i2)  (Df (0)+a2y) 

+  2a1g'ku5f(2«2-f2) 

+  la jg 'w(u2  (u»2-f2)  +4k2G2(6u2-f2)}  (A  ld) 

+  a2g'u.2f(f2-u»2) 

+  la2g,ku2u(5«2-3f2) 

2-  2  2  6  & 

-  u  u(10u  f  -3f  -7w  ) 


(A.  la) 


(A. lb) 


202 


203 


v2  -  iajg’wtk^UOw^f2)  +  u2(u»2-f2)> 

Y2  -  ig,u(Df(0)+a2yHk252(10a)2-3f2)  +  oi2(cu2-f2)} 
+  a1g’5{2k2u2(l8«2-£2)  +  w2(llu2-5f2) } 

+  ia1g,ku2f(f2-6u2) 

+  a2g'ku2u(10w2-3f2) 

+  2ia2g’u»u£(2a)2-£2) 

-2  4  224 

-  iuu  (21u  -20u  f  +3f  ) 

v  .-  2  2  2  4  2-2  2  2 

3  -  ajg'ufr  (£  -4»  )  +  5«  +k  u  (10b  -£*)} 

+  l*1g'k£u2(«2-£2) 

Y3  -  g’u(Df(0)+a2y){f2(f2-4a»2)+5tu4+k2G2(10co2-f2)} 
+  ig,kfu2(Df(0)+a2y)(a>2-f2) 

-  2ajg,kom2f 

+  iajg  'tou2  (9f  2-24k2u2-34w2) 

+  a2g’u2f (6u2-f2) 

+  ia2g'kG3(f2-10a)2) 

-  u3(35u>4+f4-20u)2f2) 

_3 

v4  •  2a^g'kfuu 

+  5ia1g'uu2(f2-2u2-k2u2) 


(A.le) 


(A.  If) 


(A.lg) 


(A.lh) 


(A.li) 
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APPENDIX  2 

PART  I  -  CORRESPONDENCE  WITH  DUXBURY ' S  NOTATION 

Th«  notation  employed  In  Duxbury's  (1963)  analysis,  which  Is 
a  special  case  of  the  problem  addressed  here,  is  sufficiently  diff¬ 
erent  that  the  following  table  is  included  in  order  to  conveniently 
summarize  the  notational  correspondence  between  the  two  models. 


Duxbury  Variable 
s 

x 

y 

o 

z 

It 

u 

a 

K 

a 

t 


Designation  in  This  Model 
al 

y 

-X 

V 

Z 

k 

-oi 

a 

(oi2-f2)/g'  +  fa^/u 
®d 

z  -  (2k/al)(Df(0)+alx) 
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appearing  in  (4.2)  are  given 


below  to  order  F' 


lo{o4(o  -l)+(kXF  )  (100-3) > 


5l(kXF  )  a(l-2o‘)  +  0(FJ) 
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}  -  (l-ff2){Fx<J2(2kA-o/Fy) 

+  (io3/Fy)  { (t/kA)  (l-o2)-kAt-F^-Fy/a} } 

,  -  kAF  o{  (a/7)  (5<j2-3)  (kAt+F  y)+2(2a2-l)} 
i  x  y  x 

+(Fx<72/Fy)  {o2-l+t(7o4+3-10o2)  } 

+  ia{a2(o2-l)+kAF2{(o2/F  ) (3-5o2)+4kA(2©2-l) >} 

*  y 

j  -  kAFx(o2(llo2-9)+(kAFxff/Fy)(3-10o2)} 
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♦  0<,*> 
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